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Abstract 


Nuclear  magnetic  resonance  (NMR)  imaging  is  a  sensitive,  non-invasive 
technique  for  monitoring  extravascular  lung  water  (EVLW),  in-vivo.  Previous 
studies  carried  out  in  this  facility  have  monitored  the  temporal  changes  in  EVLW 
content  during  the  development  of  increased  permeability  pulmonary  edema  (D. 
M.  Phillips,  M.Sc.  Thesis,  University  of  Alberta,  1987).  In  this  project  we 
continued  with  the  development  of  NMR  imaging  techniques  to  study  the 
clearance  of  lung  fluid  during  the  resolution  of  alveolar  edema. 

In  the  first  phase  of  this  study,  a  bolus  of  autologous  serum,  or  autologous 
serum  plus  platelet-activating  factor  (PAF),  was  instilled  into  one  lung  lobe  in  an 
anesthetized  cat  and  was  monitored  with  *H  NMR  imaging  for  4  hours.  PAF  is  an 
inflammatory  mediator  produced  by  the  lung  in  some  disease  states.  Images 
showed  that  serum  alone  was  cleared  slowly  from  the  lung  with  86%  of  that 
instilled  remaining  after  4  hours.  With  PAF  the  clearance  was  bi-phasic  with  a 
rapid  initial  clearance,  and  only  35%  of  the  instillate  remained  after  4  hours. 

To  investigate  whether  PAF  influenced  fluid  absorption  from  the  airspaces 
or  secretion  into  the  airspaces,  in  the  second  phase  of  this  study,  the  instillate  was 
a  50%  D2O  serum-like  solution  (SLS).  Fluid  absorption  was  monitored  with  2H 
NMR  imaging  while  the  combined  effects  of  fluid  absorption  and  secretion  were 
monitored,  as  before,  with  *H  NMR  imaging.  Without  PAF,  fluid  was  secreted 
into  the  airspaces  steadily  with  the  equivalent  of  30%  of  the  instilled  volume  over 
4  hours.  Fluid  absorption  from  the  airspaces  was  relatively  constant  when  SLS 
was  the  instillate  with  50%  absorbed  over  4  hours.  With  PAF,  however,  the 
absorption  was  markedly  bi-phasic  with  a  rapid  initial  absorption  that  slowed 
considerably  after  1  hour,  resulting  in  60%  of  that  instilled  being  absorbed  over  4 
hours. 


These  results  indicated  that  PAF  enhances  absorption  from  the 
alveolar  space  without  affecting  secretion  processes,  resulting  in  an  enhanced  net 
clearance.  Also,  these  results  are  the  First  in-vivo  observations  of  bi-directional 
fluid  movements  across  the  pulmonary  air-blood  barrier  and  demonstrate  its 
hitherto  unknown  and  unexpected  high  degree  of  activity. 
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1.0  Introduction  to  the  Thesis 


All  animals,  including  man,  depend  on  oxygen-based  metabolism  which 
provides  cells  with  the  energy  to  function.  For  most  terrestrial  animals  the  supply 
of  oxygen  originates  at  only  one  source;  the  lung,  the  prime  function  of  which  is  to 
match  ventilation  with  perfusion  (76).  To  this  end,  the  lung  has  a  series  of 
branching  airways,  matched  closely  by  the  pulmonary  arteries,  opening  into  300 
million  alveolar  sacs;  in  the  adult  human  the  alveolar  surface  area  of  70  m2  is  in 
contact  with  roughly  70  ml  of  blood  in  the  pulmonary  capillaries.  The  air-blood 
barrier  is  only  1 .5  Jim  thick,  on  average,  allowing  rapid  equilibration  (<  0.25  sec)  of 
O2  and  CO2  between  the  alveolar  air  and  the  blood. 

While  the  air-blood  barrier  of  the  lung  is  well  adapted  to  gas  exchange,  it 
also  provides  a  large  surface  area  for  fluid  movement  from  the  capillaries  into  the 
alveolar  space.  This  does  occur  when  the  alveolar-capillary  membrane  is  damaged, 
and  the  alveoli  become  flooded  with  fluid.  For  obvious  reasons,  the  balance  of 
fluid  movement  across  the  air-blood  barrier  of  the  lung,  and  the  mechanisms  which 
maintain  it,  have  long  been  a  topic  of  research.  With  the  advent  of  nuclear  magnetic 
resonance  (NMR)  imaging,  with  its  sensitivity  to  the  mobile  protons  in  water,  came 
an  exciting  means  of  monitoring  fluid  in  the  lung,  non-invasively.  Many  of  the 
earliest  attempts  to  observe  lung  water,  either  with  imaging  (18,  20,  35,  46,  65)  or 
magnetization  relaxation  time  measurements  (66,  67,  71, 83),  however,  were  only 
marginally  successful  because  of  respiratory  motions  of  the  lung,  in-vivo ,  and 
because  of  strong  local  magnetic  field  inhomogeneities  caused  by  magnetic 
susceptibility  differences  at  air-fluid  interfaces.  With  more  recent  developments  in 
NMR  hardware  and  measurement  techniques,  sensitive  measurements  of 
magnetization  relaxation  times  have  been  obtained,  in-vivo,  (68, 69,  70)  and  images 


have  provided  qualitative  observations  of  the  development  of  pulmonary  edema  (21, 
62). 
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The  original  purpose  of  this  project  was  to  improve  further  upon  available 
imaging  and  relaxation  time  measurement  techniques,  specifically  for  monitoring 
lung  water  in-vivo,  and  to  monitor  the  time  course  of  changes  of  lung  water  content 
after  a  solution  simulating  protein  rich  edema  fluid  was  introduced  into  the  alveolar 
space  of  the  lung.  This  study  explored  the  mechanisms  that  keep  the  alveoli  dry. 
We  raised  the  question  of  whether  fluid  was  only  being  absorbed  from  the  alveolar 
space,  or  if  there  was  actual  bi-directional  fluid  movement  across  the  air-blood 
barrier.  In  an  attempt  to  address  this  question,  a  second  phase  of  this  project  was 
undertaken  in  which  2H  NMR  imaging  was  developed  to  monitor  the  alveolar 
absorption  of  2H20  in  a  serum-like  solution.  Inter-leaved  with  the  2H  image 
acquisitions,  NMR  imaging  was  used,  as  in  the  first  phase  of  the  study,  to 
monitor  the  combined  effects  of  fluid  absorption  from  the  alveoli  and  secretion  of 
fluid  into  the  alveoli. 


1.1  The  Air-Blood  Barrier  of  the  Lung 

The  fluid  balance  of  the  lung  is  governed  primarily  by  the  air-blood  barrier 
which,  as  a  result,  will  be  the  focus  of  this  thesis.  For  details  of  the  structure  and 
function  of  the  lung  the  reader  is  referred  to  any  one  of  the  books  by  Staub  (76), 
Levitzky  (42),  West  (82)  and  Lim  (44). 

The  alveolar-capillary  membrane  forms  the  walls  of  the  alveoli  which  are 
polygonal  in  shape  with  a  diameter  of  approximately  300  }im,  and  have  a  total 
surface  area  of  approximately  100  m2  in  adult  humans.  The  alveolar  surface  is  lined 
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with  surfactant,  a  phospholipid  compound  which  serves  to  reduce  the  alveolar 
surface  tension  and  facilitates  expansion  of  the  alveoli  during  inhalation.  Surfactant 
also  serves  to  smooth  the  alveolar  surface  somewhat  by  filling  in  the  comers  created 
by  junctions  of  alveolar  walls  and  as  a  result  the  alveolar  surface  more  closely 
describes  a  sphere.  Alveoli  are  grouped  into  terminal  respiratory  units  which  are 
normally  supplied  with  air  via  a  single  terminal  bronchiole,  but  collateral  ventilation 
from  adjacent  terminal  respiratory  units  is  also  available  via  small  openings  in  the 
interalveolar  septae  known  as  the  pores  of  Kohn.  The  existence  of  two  or  more 
possible  routes  of  ventilation  can  greatly  influence  the  alveoli  when  there  is  regional 
edema,  particularly  if  sections  of  airways  are  filled  with  fluid  upstream  from 
partially  or  completely  air-filled  alveoli. 


Figure  1.1  A  schematic  representation  of  a  portion  of  the  pulmonary  air-blood 
barrier 


The  surface  for  gas  exchange  in  the  lung,  the  air-blood  barrier,  is  composed 
of  the  alveolar  epithelium  and  its  basement  membrane,  and  the  pulmonary  capillary 
endothelium  and  its  basement  membrane,  as  indicated  in  figure  1.1.  Two  types  of 
cells,  in  approximately  equal  numbers,  compose  the  alveolar  epithelium;  type  I  cells 
which  are  flat,  agranular,  lining  cells  covering  ~90%  of  the  alveolar  surface,  and 
type  II  cells  which  are  cuboidal  and  metabolically  active.  Type  II  epithelial  cells  are 
responsible  for  the  production  of  surfactant  and  for  maintaining  the  integrity  of  the 
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epithelium  by  differentiating  into  both  type  I  and  type  II  cells  when  required.  These 
epithelial  cells  are  joined  by  tight  junctions  and  normally  form  a  highly  impermeable 
barrier. 


The  basement  membrane  supporting  the  epithelium  is,  in  places,  fused  to  the 
basement  membrane  of  the  pulmonary  capillary  which  it  overlies,  and  otherwise 
there  is  an  interstitial  space  which  contains  interstitial  fluid,  fibroblasts,  collagen  and 
elastin  fibers,  and  interstitial  cells  (76).  It  is  generally  agreed  that  the  lymphatics, 
which  he  in  the  interstitial  space,  extend  only  to  the  level  of  the  terminal  bronchioles 
requiring  interstitial  fluid  to  be  moved  from  the  alveolar  septae  to  the  level  of  the 
lymphatics  by  a  pressure  gradient  (51, 75).  While  the  lymphatics  play  a  role  in  the 
alveolar  fluid  balance,  the  extent  of  their  role  in  the  clearance  of  edema  fluid  has 
been  shown  to  be  relatively  small  and  depends  somewhat  on  the  protein  and  solute 
concentrations  in  the  fluid  (52, 75). 

The  pulmonary  capillaries  which  underlie  70%  to  80%  of  the  surface  area  of 
the  alveolar  epithelium,  have  walls  formed  by  single  endothelial  cell  processes 
providing  a  minimal  barrier  to  the  diffusion  of  gases  into  the  blood.  Junctions 
between  these  endothelial  cells  are  not  as  tight  as  those  between  the  alveolar 
epithelial  cells  and  the  capillaries  are  likely  the  main  source  of  interstitial  fluid  and 
solutes  (51).  The  diameter  of  the  capillaries  is  such  that  red  blood  cells  must  pass 
through  in  single  file,  maximizing  the  exposure  of  each  red  blood  cell  to  the  air- 
blood  barrier.  The  supply  of  blood  to  the  alveoli  depends  both  on  gravity  and  the 
flow  resistance  of  the  supplying  arterioles  which  can  be  varied  by  changing  the  tone 
of  the  smooth  muscle  in  the  arteriole  walls.  The  pulmonary  blood  supply  is  reduced 
to  alveoli  which  are  poorly  ventilated,  resulting  in  an  increase  in  the  blood  supply  to 
better  ventilated  areas,  thereby  serving  the  lung's  prime  function  which  is  the 
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matching  of  ventilation  to  perfusion  (76).  There  are,  however,  a  number  of 
mechanisms  responsible  for  matching  ventilation  to  perfusion. 


1.2  The  Pulmonary  Fluid  Balance 

The  maintenance  of  "dry"  alveoli  for  gas  exchange  requires  a  delicate 
balance  between  the  mechanisms  affecting  fluid  pressures  on  each  side  of  the  air- 
blood  barrier.  Classically,  the  filtration  of  fluid  across  the  air-blood  barrier  has  been 
described  solely  in  terms  of  hydrostatic  and  osmotic  pressures.  This  model  was 
proposed  by  Starling  in  1896  for  fluid  filtration  across  a  semipermeable  membrane, 
in  what  has  become  known  as  the  Starling  equation: 

(1.1)  Qf=K[(p2-Pi)-c(n2-n1)] 

Here  Qf  represents  the  rate  of  fluid  filtration  across  the  membrane  and  is  positive 
when  the  flow  is  from  side  1  to  side  2  of  the  membrane,  K  is  the  membrane 
conductance,  Pi  and  P2  are  the  hydrostatic  pressures  on  sides  1  and  2  of  the 
membrane,  respectively.  The  reflection  coefficient  of  the  membrane  to  osmotically 
active  particles  is  G  which  has  values  of  1  for  an  impermeable  membrane  and  0  for  a 
permeable  membrane,  n  1  and  FI2  are  the  osmotic  pressures  on  sides  1  and  2  of  the 
membrane,  respectively.  In  applying  this  equation  to  the  capillary  endothelium  only 
protein  osmotic  pressures  need  to  be  considered  because  this  membrane  is  relatively 
permeable  to  solutes.  For  the  alveolar  epithelium,  however,  both  proteins  and 
solutes  must  be  taken  into  account  because  of  the  relative  impermeability  of  this 
membrane. 
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The  Starling  equation  has  proven  to  be  very  useful  in  describing  the 
development  of  hydrostatic  pulmonary  edema,  in  which  the  pulmonary  capillary 
pressure  is  greatly  elevated,  and  the  development  of  increased  permeability 
pulmonary  edema  in  which  the  capillary  endothelium  is  damaged  thereby  increasing 
the  membrane  conductance.  Both  of  these  situations  lead  to  an  increased  fluid  flux 
into  the  interstitial  space  resulting  in  interstitial  edema,  and  if  the  interstitial  pressure 
becomes  sufficiently  elevated,  or  the  alveolar  epithelium  is  also  damaged,  fluid  will 
move  into  the  alveolar  space  resulting  in  alveolar  edema  as  well.  The  pathogenesis 
of  pulmonary  edema  is  often  characterized  by  the  protein  concentration  of  the  edema 
fluid  relative  to  that  of  plasma  with  a  ratio  of  <0.65  indicating  that  hydrostatic 
pressure  is  the  primary  cause  and  a  ratio  of  >0.75  indicating  that  the  edema  is 
primarily  due  to  an  increase  in  membrane  permeability  (53).  What  the  Starling 
equation  fails  to  describe,  however,  is  how  alveolar  fluid  is  cleared  from  the  lung 
after  the  hydrostatic  pressure  is  reduced  to  normal  or  the  capillary  membrane  is 
repaired  and  once  again  becomes  impermeable. 

In  a  variety  of  studies  investigating  the  clearance  of  fluid  and  protein  from 
the  alveolar  space  it  was  observed  that  fluid  is  cleared  from  the  lung  much  more 
rapidly  than  albumin.  In  dogs,  sheep,  and  rabbits,  it  was  observed  that  blood  serum 
or  plasma  fluid  was  cleared  from  the  alveoli  with  clearance  half-times  (ti/2)  of  18 
hours,  9  hours  and  6  hours,  respectively,  and  that  albumin  was  cleared  at  a  rate  of 
only  ~1%  per  hour  (9,  50,  73).  The  result  of  these  vastly  different  fluid  and 
albumin  clearance  rates  is  that,  over  time,  albumin  is  concentrated  in  the  alveolar 
space  producing  a  large  osmotic  pressure  gradient  favoring  fluid  movement  into  the 
airspaces.  In  spite  of  this  pressure  gradient,  however,  the  clearance  of  fluid  from 
the  airspaces  proceeds,  and  this  led  Matthay  et  al.  (50)  to  propose  that  active 
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electrolyte  transport  might  be  an  important  mechanism  of  fluid  clearance  from  the 
alveoli. 


The  active  transports  of  Na+,  Cl'  and  K+  ions  have  been  observed  in  excised 
lungs  at  the  levels  of  the  trachea  and  the  subsegmental  bronchi  (14,  23,  59). 
Moreover,  angiotensin  II  was  observed  to  stimulate  active  Cl'  secretion  by  cultured 
tracheal  epithelial  cells  (79).  Alveolar  type  II  epithelial  cells  grown  in  culture  were 
observed  to  form  fluid  domes,  indicating  that  these  cells  are  actively  transporting 
electrolytes  with  water  following  passively  (31).  Inhibitors  of  Na+  transport  such 
as  amiloride  and  ouabain,  and  a  low  Na+  medium  were  each  observed  to  inhibit  the 
dome  formation,  while  Cl'  transport  inhibitors,  such  as  furosemide,  or  the  removal 
of  Cl"  or  K+  from  the  medium  had  no  effect  on  dome  formation  (32,  77)  indicating 
that  the  domes  were  formed  as  a  result  of  Na+  transport  from  the  apical  to  the 
basolateral  side  of  the  epithelial  cells.  Also,  a  variety  of  factors  such  as  (3- 
adrenergic  agonists,  cAMP  analogues  and  phosphodiesterase  inhibitors  were 
observed  to  increase  dome  formation,  indicating  that,  if  active  Na+  transport  is 
important  to  fluid  clearance  in-vivo ,  it  may  be  affected  by  nervous  stimulation  or  the 
supply  of  intracellular  cAMP  (33,  49). 

In  an  effort  to  determine  if  the  active  electrolyte  transport  observed  in 
cultured  epithelial  cells  is  a  significant  mechanism  for  clearing  fluid  from  the 
alveolar  space,  Basset  et  al.  carried  out  a  series  of  experiments  on  isolated  rat  lungs 
which  were  perfused  with  Ringer's  solution,  and  had  Ringer's  solution  instilled  into 
the  airspaces  (3,  4).  They  observed  both  passive  Na+  movements,  assumed  to  be 
paracellular,  and  active  Na+  transport  across  the  alveolar  epithelium  via  Na+ 
channels  and  a  Na+-D-glucose  symport  across  the  luminal  cell  membrane,  as  well 
as  a  Na+-K+  pump  across  the  abluminal  membrane.  Also,  they  observed  that  there 
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may  be  a  luminal  Na+-K+  pump  as  well,  suggesting  that  two  types  of  epithelial  cells 
actively  transport  electrolytes  as  it  is  unlikely  that  one  cell  would  have  both  luminal 
and  abluminal  Na+-K+  pumps  (19).  More  recently,  similar  results  were  observed 
in-vivo  when  Smedira  et  al.  (73)  instilled  plasma  into  the  airspaces  of  anesthetized 
rabbits  and  demonstrated  that  75%  of  the  fluid  clearance  could  be  inhibited  with  the 
addition  of  amiloride  to  the  instillate.  Also,  Berthiaume  stimulated  Na+  transport 
across  alveolar  epithelial  cells  in  anesthetized  sheep  by  raising  lung  tissue  cAMP 
levels,  and  as  a  result  observed  an  increased  lung  fluid  clearance  (10). 

The  clearance  of  protein  from  the  alveolar  space  presented  even  more  of  a 
puzzle  than  the  fluid  clearance  as  albumin  was  observed  to  be  cleared  from  the 
lungs  of  anesthetized  animals  at  a  rate  of  ~1  %/hour  and  yet  the  tight  epithelium  was 
intact  and  considered  impermeable  to  macromolecules  (9).  In  1947  Drinker  et  al 
(22)  proposed  that  proteins  must  be  broken  down  by  enzymes  into  fragments  small 
enough  to  cross  the  epithelium.  This  view  held  until  1967  when  Bensch  et  al.  (5) 
observed  that  the  lung  can  absorb  intact  proteins.  Bensch  et  al.  later  expanded  these 
findings  when  they  observed  that  macromolecules  can  be  transported  to  the 
capillary  blood  by  pinocytotic  vesicles  in  type  I  epithelial  cells  and  capillary 
endothelial  cells  and  that  the  lymphatics  play  only  a  minor  role  in  the  protein 
clearance  from  the  lung  (6). 

All  of  the  above  studies  indicate  that  in  addition  to  hydrostatic  and  osmotic 
forces,  the  pulmonary  fluid  balance  is  also  affected  by  the  active  transport  of  ions 
and  macromolecules  across  the  alveolar  epithelium.  Under  normal  circumstances 
the  lung  must  require  a  balance  between  all  of  these  mechanisms  to  maintain  a 
steady  state  of  "dry"  alveoli.  The  development  of  lung  edema  likely  depends,  for 
the  most  part,  on  hydrostatic  and  osmotic  pressures  as  it  has  been  shown  that  fluid 
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absorption  mechanisms  become  inoperative  if  the  normally  high  reflection 
coefficients  for  Na+  and  Cl'  are  lowered  (4).  The  clearance  of  edema  fluid  during 
the  recovery  from  a  pathological  state,  however,  is  likely  dominated  by  the  active 
transport  mechanisms  that  have  been  demonstrated  across  the  air-blood  barrier  of 
the  lung. 

Platelet-activating  factor  (PAF)  is  a  family  of  ether-lipids  whose  active 
component  is  l-0-Hexadecyl/octadecyl-2-acetyl-sn-glycero-3-phosphocholine  and 
is  produced  by  a  wide  variety  of  cells  in  the  body  including  alveolar  type  II 
epithelial  cells  (36,  37,  58, 63).  It  has  been  shown  that  intra-airway  PAF  can  cause 
bronchoconstriction,  inflammation  in  the  airways  and  capillary  spaces,  and 
pulmonary  neutrophil  sequestration  (11,  78,  80).  The  effects  of  PAF  on  the 
clearances  of  alveolar  fluid  and  protein  are  addressed  in  this  thesis  because  this 
agent  is  suspected  of  playing  roles  in  asthma  and  in  sarcoidosis,  and  the  role  of 
PAF  in  the  development  and/or  resolution  of  pulmonary  edema  has  not  been 
previously  evaluated  (2,  64). 


1.3  NMR  Studies  of  Lung  Water 

1.3.0  A  brief  historical  overview 

The  means  by  which  an  image  of  an  object  can  be  formed  from  an  NMR 
signal  was  first  described  in  1973  by  Lauterbur  (40)  and,  independently,  by 
Mansfield  and  Grannell  (47).  Three  years  later  Lauterbur  (41)  and  Frank  (27) 
concluded  that  it  may  be  feasible  to  use  NMR  imaging  to  detect  pulmonary  edema 
when  they  observed  a  linear  relationship  between  the  longitudinal  relaxation  rate 
and  the  dry  weight  to  wet  weight  ratio  of  lung  tissue.  Since  that  time  a  variety  of 
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studies  have  been  carried  out  showing  that  both  longitudinal  (Ti)  and  transverse 
(T2)  relaxation  times  increase  with  increasing  lung  water  (62,  66,  67,  71,  83)  and 
that  relaxation  times  can  be  used  to  differentiate  between  hydrostatic  and 
permeability  edema  types  (66).  The  NMR  signal  intensity  has  also  been  used  to 
detect  regional  edema  and  it  has  been  shown  to  be  well  correlated  with  the 
extravascular  lung  water  content  under  a  variety  of  conditions  (18, 20,  35,46  ,65). 

All  of  these  in-vivo  studies  suffered,  however,  from  the  typically  low  signal 
to  noise  ratio  obtained  from  the  lung  because  of  the  paucity  of  water  protons  in  the 
normal  lung,  short  transverse  relaxation  times,  plus  cardiac  and  respiratory  motions. 
The  means  to  reduce  the  detrimental  effects  of  these  properties  began  to  develop  as 
an  understanding  of  the  NMR  properties  of  lung  water,  and  more  sophisticated 
NMR  hardware,  were  developed.  For  example,  it  was  determined  that  air-water 
interfaces  in  the  lung  can  cause  local  magnetic  field  gradients  due  to  magnetic 
susceptibility  differences  across  the  interface,  and  that  these  local  gradients  enhance 
transverse  relaxation  effects.  However,  the  effects  of  local  magnetic  field  gradients 
can  under  most  circumstances  be  reversed  with  an  NMR  pulse  sequence,  which 
will  be  discussed  in  section  1.3.3.  Further  understanding  of  the  relaxation 
mechanisms  which  are  dominant  in  the  lung  was  provided  by  the  observations  of 
Kveder  et  al  (38).  Firstly,  they  observed  that  lung  water  underwent  fairly  restricted 
movement  and  diffused  much  more  slowly  than  free  water.  Secondly,  they 
observed  that  water  proton  relaxation  in  the  lung  is  influenced  by  the  rapid 
exchange  of  water  between  a  free  state,  and  a  state  in  which  it  is  bound  to  bio¬ 
polymer  segments  at  the  lung  tissue  surface.  As  a  result,  Kveder  et  al.  postulated 
that  longitudinal  relaxation  rates  of  lung  water  protons  should  be  proportional  to  the 
lung  tissue  surface  area  (39).  Moreover,  for  improving  the  quality  of  lung  images 
various  techniques  for  respiratory  gated  signal  acquisition  have  been  described  (24). 
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Also,  Phillips  et  al.  have  proposed  that  the  optimal  NMR  pulse  sequence  for  lung 
imaging  and  relaxation  time  measurements  is  a  spin-echo  sequence  which  includes 
respiratory  gated  signal  acquisition  with  regular  radio-frequency  excitation  pulses  to 
maintain  a  steady  state  (spin-conditioning),  and  excludes  cardiac  gated  signal 
acquisition  (61,  62). 

Making  use  of  these  developments  in  NMR  techniques,  Phillips  et  al.  were 
able  to  observe  the  changes  in  relaxation  times  that  occurred  as  permeability  edema 
developed  in-vivo  (62).  Moreover,  Shioya  et  al.  were  able  to  achieve  a  high  enough 
signal  to  noise  ratio  to  observe  the  bi-exponential  transverse  relaxation  of  water 
protons  in  normal  lung  tissue  in-vivo,  and  therefore  to  measure  the  transverse 
relaxation  times  of  two  distinct  fluid  components  in  the  lung  (68, 69, 70).  AJthough 
much  progress  has  been  made  in  the  use  of  NMR  for  the  study  of  lung  edema,  the 
bulk  of  the  studies  that  have  been  carried  out  were  aimed  at  demonstrating  that  some 
property  of  the  NMR  signal  or  some  technique  could  be  useful  for  detecting  and 
monitoring  edema.  Of  the  studies  that  had  been  done  prior  to  this  project,  only  four 
had  actually  involved  monitoring  extravascular  lung  water  or  one  of  its  properties 
for  a  period  of  time  (21, 62, 67, 70). 

The  capabilities  of  current  NMR  hardware  and  techniques,  for  both  clinical 
and  research  applications,  are  demonstrated  by  the  following  examples.  Firstly, 
recent  developments  of  a  projection  reconstruction  technique  for  reducing  the 
effects  of  motion  and  magnetic  susceptibility  differences,  have  enabled  the 
acquisition  of  images  which  demonstrate  the  large  blood  vessels  and  bronchi  in  the 
lung  (7).  Magnetic  resonance  imaging  has  been  shown  to  detect  experimentally 
induced  changes  in  the  lung  ventilation  and  perfusion  distributions  in  an 
anesthetized  animal,  with  the  administration  of  intra-airway  and  intravenous  contrast 
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agents  (8).  Clinically,  magnetic  resonance  imaging  has  been  shown  to  depict  the 
presence  or  absence  of  hilar  lymph  nodes,  and  to  demonstrate  differences  between 
patients  with  and  without  Pneumocystis  carinii  pneumonia  (54,  81).  Moreover, 
NMR  imaging  has  been  shown  to  be  useful  for  the  characterization  of  pleural  fluid 
and  is  able  to  distinguish  pleural  effusion  from  parenchymal  disease  and  from 
pleural  tumors  (56).  Relaxation  time  measurements,  and  relaxation  time  weighted 
imaging,  have  also  demonstrated  differences  in  pleural  effusions,  as  well  as 
demonstrating  the  extent  of  chest  wall  invasion  by  lung  tumors,  and  the  presence  of 
bronchogenic  carcinoma  and  distal  airway  collapse  (54,  56,  81).  While  the  ability  to 
characterize  regions  of  increased  signal  intensity  in,  or  adjacent  to,  the  lung  is  well 
demonstrated  by  these  studies,  we  have  found  no  clinical  NMR  studies  concerning 
the  quantification  of  extravascular  lung  water  discussed  in  the  literature. 

1.3.1  Macroscopic  Nuclear  Magnetization 

The  induction  of  nuclear  magnetization  in  a  spin  ensemble  by  means  of  an 
external  magnetic  field  requires  only  that  one  or  more  species  of  nuclei  comprising 
the  system  have  magnetic  moments.  The  nuclear  magnetic  moment,  |d,  is  related  to 
the  nuclear  spin,  I,  by  the  equation 

(1.2)  |d  =  yh  I 

where  y  is  the  gyromagnetic  ratio  and  is  unique  for  each  nuclear  species,  and  h  is 
Planck's  constant  divided  by  27T.  The  interaction  between  the  magnetic  moment  of  a 
nucleus,  Jd,  and  the  external  magnetic  field,  B0,  known  as  the  Zeeman  interaction, 
has  a  Hamiltonian,^,  given  by 
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(13)  tfz  =  -H-B0 

Taking  B0  to  be  along  the  z  axis,  and  substituting  for  JI  from  equation  1 .2,  results 
in  equation  1.4. 

(1-4)  #z  =  -yhB0  Iz 

where  Iz  denotes  the  z  component  of  the  nuclear  spin.  Substituting  this  Hamiltonian 
into  the  Schroedinger  equation  yields  a  solution  with  21+1  stationary  eigenstates 
with  energies,  E,  given  by 

(1.5)  E  =  -  y  h  B0  mz 

where  mz  represents  the  eigenvalues  of  the  operator  Iz,  and  has  values  -I,  -1+1, ... , 
1-1, 1.  For  a  spin  I  =  1/2  nucleus  there  are  two  energy  eigenstates,  identified  by 
their  orthonormal  eigenfunctions  I  a  >  and  I  (3  >,  corresponding  to  mz  values  of  1/2 
and  -1/2,  respectively.  In  the  general  case,  however,  the  spin  state  of  a  nucleus 
cannot  be  described  by  one  or  the  other  of  these  energy  eigenfunctions,  but  rather  is 
an  admixture  of  the  two. 

For  an  ensemble  of  non-interacting  nuclear  spins  in  thermal  equilibrium,  the 
relative  number  of  spins  in  each  energy  state  is  given  by  a  Boltzmann  distribution 
with  more  spins  in  the  lower  energy  states  than  in  the  higher  energy  states.  In  the 
case  of  spin  I  =  1/2  nuclei,  the  number  of  spins  in  state  I  P  >,  Np,  relative  to  the 
number  in  state  I  a  >,  Na,  is 

(Np/Na)  =  exp 


-AE 

kT 


(1.6) 
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where 

(1.7)  AE  =  Ep  -  Ea  =  y  h  B0 

In  equation  1.6,  k  is  Boltzmann's  constant  and  T  is  the  absolute  temperature.  The 
population  difference  between  these  two  energy  states  results  from  the  probability 
of  finding  a  spin  aligned  parallel  to  B0  being  greater  than  that  of  finding  it 
antiparallel  to  B0,  and  for  an  ensemble  of  N  spins  results  in  a  net  z-component 
magnetization  given  by 

(1.8)  Mz  =  yr,  £<iz> 

all  spins 

where  <  Iz  >  denotes  the  expectation  value  of  Iz  for  a  single  nuclear  spin.  For  an 
ensemble  of  N  spins  in  thermal  equilibrium,  the  net  z-component  magnetization  has 
magnitude  M0  (1) 

(1.9)  M0  =  Ny2h2l(I+l)Bo/3kT  =  XoB0 

In  this  expression  Xo  is  known  as  the  static  nuclear  susceptibility.  In  thermal 
equilibrium,  there  is  no  preferred  direction  for  the  transverse  component  of  an 
individual  spin  magnetic  moment,  and  all  transverse  spin  component  orientations 
can  be  found  with  equal  probability.  As  a  result,  the  transverse  magnetization 
components,  Mx  and  My,  are  zero  and  there  is  no  detectable  NMR  signal.  (The 
detection  of  the  NMR  signal  is  discussed  later  in  this  section.)  However,  in  the 
general  case 
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(i-io)  Mx  =  yii  £<ix> 

all  spins 

(1.11)  My  =  yh  £<ly> 

all  spins 

In  these  expressions  <  Ix  >  and  <  Iy  >  denote  the  expectation  values  of  Ix  and  Iy  for 
a  single  nuclear  spin.  Alternatively,  one  can  combine  equations  1.10  and  1.11,  and 
substitute  the  raising  operator,  I+,  for  Ix  +  i  Iy,  to  obtain  the  expression  for  the 
transverse  magnetization  of  the  spin  ensemble,  M+,  in  equation  1.12. 

(1.12)  M+  =  yh  X<I+> 

all  spins 

Because  a  nuclear  spin  state  can  generally  be  described  by  an  admixture  of  its 
energy  eigenfunctions,  the  spin  ensemble  expectation  values  of  Ix,  Iy,  and  I+,  can  be 
non-zero  when  the  spin  system  is  not  in  thermal  equilibrium.  For  example,  the 
single  spin  state,  \| /,  can  generally  be  written 

(1.13)  \|/  =  C+la>  +  C.lp> 
subject  to  the  condition  that 


(1.14) 


C+2  +  C_2  =  1 
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The  expectation  value  <  1+  >  can  then  be  computed  as  shown  in  equations  1.15  to 
1.19,  with  the  use  of  Dirac  notation  (43).  Firstly,  <  1+  >  is  given  by 

(1.15)  <I+>  =  <\j/II+l\j/> 

(43)  and  substituting  for  \J /  from  equation  1.13  gives 

(1.16)  <  1+  >  =  <  \J/  I  I+l  (  C+  I  a  >  +  C.l  (3  >  ) 

Applying  the  raising  operator  then  yields 

(1.17)  <  1+  >  =  <  \\f  I  C.l  a  > 

and  substituting  for  \|/  from  equation  1.13  results  in  the  expression 

(1.18)  <  4  >  =  (  C+*  <  a  I  +  C*  <  P  I )  C.l  a  > 

where  *  denotes  the  complex  conjugate.  Because  I  a  >  and  I  (3  >  are  orthogonal 
and  normalized,  equation  1.18  can  be  simplified  further  to  the  expression 

(1.19)  <I+>  =  C+*C. 

As  a  result,  the  quantum  mechanical  treatment  of  the  nuclear  spin  interaction  with  a 
magnetic  field  can  be  used  to  describe  the  net  transverse  magnetization,  M+ 
(equation  1.12),  of  a  spin  ensemble  and  is  not  limited  to  describing  magnetization 
which  is  parallel  or  antiparallel  to  a  static  B0  field. 
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To  describe  the  motion  of  the  net  magnetization,  M,  of  an  ensemble  of  spins 
in  any  external  magnetic  field,  Felix  Bloch  (12)  used  a  phenomenological  approach 
to  develop  a  set  of  equations  that  turn  out  to  be  valid  for  most  biomedical 
applications  of  NMR  because  of  the  highly  fluid  state  of  animal  tissue.  In  low 
viscosity  liquids  such  as  water,  extracellular  fluid,  or  plasma  for  example,  the 
magnetic  field  that  a  spin  experiences  due  to  adjacent  spin  magnetic  moments 
averages  to  zero  over  time  because  of  the  random  thermal  motions  of  the  spins,  and 
the  system  magnetization  approaches  its  equilibrium  value  exponentially  as  assumed 
by  Bloch.  The  basis  of  Bloch’s  approach  is  that  a  magnetic  moment,  JT,  interacting 
with  a  magnetic  field,  B,  will  experience  a  torque  given  by 


r  =  ji  x  b 


(1.20) 


Substituting  the  rate  of  change  of  the  spin  angular  momentum,  hi,  for  the  torque,  T, 
yields 


(1.21) 


,  dl 

hdt  =  ^xB 


With  the  additional  substitution  of  JLl/yh  for  I,  from  equation  1 .2,  gives 


(1-22) 


d-23) 


In  order  to  apply  this  same  treatment  to  a  spin  ensemble,  one  can  substitute  for  the 
magnetic  moment  jd,  its  expectation  value  <Jd>  (72).  Furthermore,  if  the  spins  in 
the  ensemble  can  be  thought  of  as  not  interacting  with  one  another,  equation  1 .23 
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holds  true  for  the  magnetic  moment  per  unit  volume,  i.e.  the  total  magnetization  of 
the  spin  system,  M  (72).  The  substitution  of  M  for  JI  in  equation  1.23  results  in 
the  expression 


(1.24) 


Without  a  loss  in  generality  one  can  define  the  z  axis  of  the  coordinate  system  to  be 
parallel  to  B.  The  solution  of  the  differential  equation  1.24  is  then  that  the  x  and  y 
components  of  M  precess  about  the  z  axis  at  a  frequency  C0o,  called  the  Larmor 
frequency,  which  is  equal  to  -yB. 

Because  the  spin  ensemble  magnetization  precesses  about  an  externally 
applied  magnetic  field,  the  interaction  between  M  and  B  is  much  easier  to  describe 
in  a  reference  frame  which  is  rotating  about  the  z  axis  at  the  frequency  C0o.  To 
allow  for  the  more  general  case,  however,  we  describe  this  interaction  in  a  reference 
frame  rotating  about  the  z  axis  at  a  frequency  Q,  and  then  investigate  the  special 
case  in  which  £2  =  C0o.  In  differentiating  M  with  respect  to  time  in  a  rotating 

frame,  it  is  necessary  to  allow  for  the  general  case  in  which  the  orthogonal  unit 
vectors  i,  j,  and  k  may  be  rotating  with  an  instantaneous  angular  velocity  Q  k  (72). 

In  this  case,  the  transformation  from  the  laboratory  frame  to  the  rotating  frame 
proceeds  as  follows.  Firstly, 


(1.25) 


and  in  a  rotating  frame  the  time  derivative  of  a  unit  vector  is  given  by 


19 


(1.26) 


di_Q 

dt  “i2 


x  1 


with  similar  expressions  for  the  derivatives  of  the  j  and  k  unit  vectors.  Substituting 
expressions  for  di/dt,  dj/dt  and  dk/dt  into  equation  1 .25  and  simplifying,  yields 


(1.27) 


dM  8M 


dt 


5t 


+  Q  x  M 


dM 


Now  by  substituting  for -jp-  from  equation  1 .24  results  in  the  expression 


(1.28) 


SM  w  „  ^ 

- =  yM  xBq-Q  x  M 

8t 


which  can  be  rearranged  to  give 


(1.29) 


8M 

St 


=  yM  x  (B0  +  Q/y) 


With  the  introduction  of  an  effective  magnetic  field  in  the  rotating  frame,  Beff, 
defined  by 

(1.30)  Beff  =  B0  +  ^/y 


equation  1 .29  is  reduced  to 


SM 

St 


=  yM  X  Beff 


(1.31) 
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Depending  on  the  magnitude  and  sense  of  f>,  Beff  can  be  larger  or  smaller  than  B0. 
In  the  on-resonance  case  the  reference  frame  rotates  at  the  angular  velocity  of  the 
precessing  magnetization,  M,  as  described  by  equation  1.32.  Thus,  in  the  on- 
resonance  rotating  frame  the  effects  of  the  B0  field  are  nullified  and  the 
magnetization  appears  to  be  static,  as  described  by  equation  1 .33. 


(1.32)  £2  =  C0o  =  -YBo 


(1.33) 


8M 

- =  0 

St 


If  one  superimposes  onto  the  static  B0  field  an  RF  field,  Bi,  which  is 
directed  along  the  Xp  axis  in  the  rotating  frame,  and  Q  is  chosen  to  be  the  same  as 
the  rotational  frequency  of  the  Bi  field,  the  resultant  effective  field  becomes 

(134)  Beff  =(B0  +  n/Y)kp  +  Biip 

Where  the  unit  vectors  ip, ip  and  kp  are  along  the  Xp,  yp  and  zp  axes,  respectively. 
Substituting  equation  1.34  into  equation  1.31  now  yields  the  new  equation  of 
motion  in  the  rotating  reference  frame: 

O  -E  K 

( 1 35)  —  =  yM  x  [(B0  +  f>/y)kp  +  B  i  ip] 

5t 

As  a  result,  the  magnetization,  M,  precesses  about  the  new  effective  field  in  the 
rotating  frame  at  a  frequency,  C0eff,  given  by  equation  1.36  and  as  shown 
schematically  in  figure  1.2. 
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Figure  1.2  Motion  of  a  system  of  non-interacting  spins  in  a  frame  of  reference 
rotating  at  angular  velocity  Q  in  the  presence  of  an  applied  Field,  Bj,  rotating  at 
angular  velocity  Q.  The  system  was  initially  at  thermal  equilibrium  with  a  net 
magnetization,  M,  along  the  z  axis. 


(1 .36)  COeff  -  y[(B0  +  F>/y)kp  +  Bi  ip]  -  yBeff 


In  the  special  case  in  which  the  B]  field  is  on-resonance,  Q  =  -C0o,  the  effective 
field  in  the  rotating  frame  is  equal  to  Bi  ip  and  the  magnetization  precesses  in  the 
ypZp  plane  as  shown  in  figure  1 .3. 
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Figure  1.3  Motion  of  a  system  of  non-interacting  spins  in  a  frame  of  reference 
rotating  at  angular  velocity  Cl  =  -yB0  in  the  presence  of  an  applied  field,  Bi, 
rotating  at  angular  velocity  f>.  The  system  was  initially  at  thermal  equilibrium  with 
a  net  magnetization,  M,  along  the  z  axis. 

As  illustrated  in  figures  1.2  and  1.3,  if  the  Bi  field  rotates  at  a  frequency  near  C0o 
and  is  applied  for  the  appropriate  duration,  the  magnetization  can  be  rotated  90° 
from  the  z  axis  into  the  Xpyp  plane.  In  practice,  only  brief  pulses  of  the  Bi  field,  on 
the  order  of  1  ps  to  hundreds  of  ps,  are  required  to  rotate  the  ensemble 
magnetization  through  90°  and  so  pulses  of  this  type  are  referred  to  as  90°  RF 
pulses.  Only  when  the  Bj  field  is  on-resonance  as  depicted  in  figure  1.3,  can  the 
magnetization  be  rotated  (with  a  180°  RF  pulse)  from  the  +z  axis  to  the  -z  axis.  The 
general  expressions  for  the  magnitude  of  the  magnetization  components  in  the  xpyp 
plane  and  along  the  z  axis  after  an  RF  pulse  of  duration  t  are  shown  in  equations 
1 .40  and  1 .41 .  The  derivation  of  these  expressions  is  included  in  Appendix  1 .1 . 


(1.37) 


Beff  =  [(Bo  +  C2/Y)2  +Bl2]1/2 
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(1.38) 


(1.39) 


(1.40) 


Mxy  =  MosinO  [cos26  (1  +  cos(yBeff  t))2  +  sin2(yBeff  t)] 1/2 


(1.41) 


Mz  =  Mq  [cos29  -  sin29  cos(yBeff  t)] 


These  equations  demonstrate  that  by  applying  an  on-resonance  RF  field,  Bi,  for  the 
appropriate  duration,  t,  the  magnetization  of  a  spin  system  can  be  rotated  through 
any  desired  angle  about  Bj.  The  means  to  manipulate  the  spin  system 
magnetization  in  this  way  is  employed  in  all  NMR  applications,  and  is  an  integral 
part  of  the  NMR  imaging  theory  discussed  in  section  1.3.3. 

Bloch  also  introduced  into  the  equation  of  motion  of  M,  two  terms  to 
describe  the  relaxation  of  the  magnetization  back  to  thermal  equilibrium,  subject  to 
the  conditions  that  inter-nuclear  interactions  are  weak  and  that  the  major  cause  of 
relaxation  is  thermal  agitation  of  the  spins  (12).  The  first  of  these  terms  describes 
magnetization  relaxation  in  the  longitudinal  direction,  parallel  to  B0,  whereas  the 
second  term  describes  relaxation  of  the  transverse  component  magnetization,  normal 
to  B0.  The  longitudinal  component  of  the  magnetization,  Mz,  was  assumed  by 
Bloch  to  approach  its  thermal  equilibrium  value  of  M0  exponentially  with  a 
characteristic  time,  Ti,  as  described  by  the  differential  equation 


dMz  Mz  -  M0 


dt  ’  Tj 


(1.42) 


The  transverse  magnetization  was  similarly  assumed  by  Bloch  to  return  to  its 
equilibrium  value  of  zero,  exponentially  with  a  characteristic  time,  T2,  so  that  the 
rate  of  change  of  Mx  and  of  My  are  described  by  the  equations 


24 


(1-43) 


dMx  _  Mx 


dt  -'T2 


(1.44) 


dMv  _  M 


dt  " '  T2 


Bloch  then  assumed  that  the  effects  of  nuclear  relaxation  could  be  superimposed  on 


the  motions  of  the  spins  due  to  externally  applied  magnetic  fields,  to  yield  a  general 


equation  of  motion  for  a  system  of  weakly  interacting  nuclear  magnetic  moments 
such  as  in  a  liquid  (equation  1 .45). 


(1-45) 


Equation  1.45  enables  one  to  describe  the  interaction  between  the  macroscopic 
nuclear  magnetization  and  externally  applied  magnetic  fields  in  most  biological 
systems.  Also,  this  equation  can  describe  the  basic  NMR  experiment  in  which  a 
90°  RE  pulse  is  applied  to  rotate  the  magnetization  of  a  spin  ensemble  from  thermal 
equilibrium  into  the  transverse  plane.  As  the  resultant  transverse  component  of  the 
magnetization  relaxes  back  to  thermal  equilibrium  it  also  precesses  about  the  static 
B0  field.  If  this  precessing  magnetization  is  linked  by  a  receiver  coil,  it  will  induce 
in  the  coil  an  oscillating  electromotive  force  (EMF),  which  in  turn  will  produce  a 
time  varying  potential  difference,  V(t),  across  a  load  at  the  output  of  the  coil.  It  is 
this  potential  difference,  V(t),  that  is  the  NMR  signal.  Descriptions  of  the  NMR 
hardware  and  methods  of  producing  static  and  RF  magnetic  fields,  and  detecting  the 
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NMR  signal  have  been  provided  in  great  detail  by  a  variety  of  authors  and  so  will 
not  be  repeated  here  (1, 28, 72). 

1.3.2  Magnetization  Relaxation  Processes 

The  longitudinal  relaxation  of  nuclear  magnetization  back  to  thermal 
equilibrium,  parallel  to  the  static  Bc  field,  is  also  known  as  "spin-lattice"  relaxation. 
The  term  "spin-lattice"  is  used  because  this  relaxation  is  due  to  a  transfer  of  energy 
from  the  excited  spins  to  the  thermal  energy  of  the  surrounding  medium,  termed  the 
"lattice".  Transverse  relaxation,  on  the  other  hand,  is  also  known  as  "spin-spin" 
relaxation.  In  this  case  "spin-spin"  refers  to  the  interaction  between  spin  magnetic 
moments  and  adjacent  spin  magnetic  moments  or  paramagnetic  ions.  Bloch  (12) 
however,  used  an  all  encompassing  description  of  transverse  relaxation  as  being 
due  to  inter-nuclear  interactions,  interactions  with  paramagnetic  ions,  and  with 
inhomogeneities  in  the  static  B0  field.  The  random  interactions  between  a  nuclear 
spin  and  another  spin  or  a  paramagnetic  ion,  cause  both  a  change  in  the  transverse 
components  of  the  spin  expectation  values  and  a  loss  of  coherence  of  their 
precession,  leading  to  an  irreversible  decay  of  the  transverse  magnetization  towards 
its  equilibrium  value  of  zero.  This  decay  is  characterized  by  the  time  T2.  Static  field 
inhomogeneities,  however,  result  in  a  dephasing  of  the  transverse  magnetization 
components  which  is  spatially  dependent  and  (in  the  absence  of  rapid  diffusion)  is 
reversible.  Because  the  transverse  relaxation  described  by  Bloch  includes  effects  of 
the  static  B0  field  inhomogeneities,  it  is  characterized  by  the  relaxation  time  T2* 
which  is  less  than  or  equal  to  the  natural  T2  of  the  spin  system. 

For  the  purposes  of  this  thesis  only  two  nuclear  species  are  of  interest, 
namely,  the  proton  and  the  deuteron,  and  so  the  following  discussions  of  nuclear 


26 


relaxation  will  be  limited  to  these  nuclei.  For  the  spin  I  =  1/2  proton,  relaxation  is 
dominated  by  the  magnetic  interaction  between  adjacent  dipole  moments,  known  as 
the  dipole-dipole  interaction.  The  deuteron,  on  the  other  hand,  has  a  spin  of  1  and 
so  has  an  electric  quadrupole  moment  which  can  interact  with  a  local  electric  Field 
gradient.  For  deuterons  in  a  liquid-like  system  it  is  this  quadrupolar  interaction 
which  dominates  the  nuclear  relaxation. 

Dipole-Dipole  Relaxation 

The  relaxation  process  which  has  been  termed  dipole-dipole  relaxation  is 
driven  by  the  interaction  between  a  spin  magnetic  dipole  moment  and  local  time 
dependent  magnetic  fields  which  are  due  to  adjacent  spin  magnetic  moments  or 
paramagnetic  ions.  The  time  dependent  nature  of  these  local  magnetic  fields  is 
provided  by  the  random  thermal  motions  of  the  spins,  and  by  the  natural  flip  rate  of 
the  ionic  spins.  It  is  the  magnetic  coupling  of  adjacent  spin  magnetic  moments  that 
provides  a  means  of  transferring  energy  between  spins,  and  from  spins  to  the 
thermal  lattice  so  that  the  spin  system  and  the  lattice  can  approach  a  common 
thermal  equilibrium. 

The  spin  wave  function  for  a  two  spin  system  in  a  static  magnetic  field, 
similar  to  that  for  the  single  spin  system  (equation  1.13),  can  be  described  as  a 
linear  combination  of  the  following  stationary  eigenstates  of  the  Zeeman 
Hamiltonian  (equation  1.3): 

(1.46)  \\f\  =  laa> 

Vo,s  =  I  a(3  >  +  I  pa  >  ] 


(1.47) 
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(1.48) 

Vo, A  =  1  ap  >  -  1  (5a  >  ] 

(1.49) 

V-i  =  1  PP  > 

Moreover,  transition  selection  rules  do  not  permit  the  system  to  undergo  transitions 
between  the  symmetric  spin  eigenstates,  \| f\,  \|/o,S  and  and  the  antisymmetric 

state  Xj/o^A  (25).  As  a  result,  we  can  express  a  general  two-spin  wave  function  *F  as 
a  linear  combination  of  energy  eigenfunctions  (similar  to  equation  1 .13),  as  follows 

( 1 .50)  *F  =  a  i  (t)\|/i  +  ao(t)\|/o  ,s  +  a.  i  (t)\| f.  \ 

l 

(1.51)  VT=San(t)Vu 

n=-l 

where  the  coefficients,  an(t),  are  given  by 

( 1 .52)  an(t)  =  Cn  exp(-i  En  t  /  h) 

In  this  expression,  the  energy  eigenvalue,  Eq,  corresponds  to  the  eigenstate  \J/n,  and 
the  coefficient  Cn  is  time  dependent,  as  will  be  demonstrated  below.  In  the  special 
case  in  which  the  two-spin  system  is  in  thermal  equilibrium  with  the  lattice,  it  is 
most  probable  that  the  system  will  be  in  the  lowest  energy  spin  eigenstate,  \\f[. 
Thus,  for  the  purposes  of  describing  nuclear  relaxation,  we  are  interested  in  the  rate 
at  which  a  spin  system  undergoes  transitions  from  the  states  \|/o ,s  and  to  the 
state  \j/p 
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The  Hamiltonian  of  the  dipole-dipole  interaction,  similar  to  that  of  the 
Zeeman  interaction,  depends  on  the  nuclear  magnetic  moment,  |I,  and  the  magnetic 
field,  b,  experienced  by  the  spin  as  a  result  of  an  adjacent  magnetic  dipole. 
Therefore,  to  derive  the  dipole-dipole  Hamiltonian  we  proceed  as  follows. 

The  magnetic  field,  b,  due  to  a  magnetic  dipole,  JLL,  is  described  in  spherical- 
polar  coordinates  in  equations  1.53  to  1.56,  and  the  coordinate  system  used  is 
illustrated  in  figure  1.4  (45). 


(1-53) 


b  =  br  r  +  b0  £  +  b^  £ 


(1.54) 


(1.55) 


(1.56) 


where  JIQ  is  the  permeability  of  free  space,  and  r,  9  and  <J>  are  axial  unit  vectors  in 
the  spherical-polar  coordinate  system.  Alternatively,  this  magnetic  field  can  be 
expressed  as  (45) 


(1.57) 
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x 

Figure  1.4  Spherical-polar  coordinate  system  used  to  describe  the  magnetic  Field 
due  to  a  magnetic  dipole  moment,  1 1. 

Similar  to  the  energy  of  the  Zeeman  interaction,  given  by  equation  1 .2,  the  energy  of 
the  interaction  between  a  magnetic  moment  jd2  and  the  field  b  due  to  the  magnetic 
moment  ]Ii  is  given  by: 


(1.58) 


where  r  is  now  the  distance  between  the  two  nuclei  and  r  is  a  unit  vector  pointing 
from  one  nucleus  to  the  other.  Because  of  the  relationship  between  the  nuclear  spin 
of  a  nucleus  and  its  magnetic  moment,  this  energy  can  be  expressed  in  terms  of  the 
nuclear  spins  as: 


(1-59) 
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Expanding  this  equation  further  in  terms  of  spin  operators  yields  the  Hamiltonian 
for  the  dipole-dipole  interaction,  as  shown  in  equation  1 .60  ( 13) 


(1.60) 


where: 


A  =  Izi  Iz2(1-3cos20i2) 

B  =  (- 1/ 4) [I_ !  42  +  l+i  L2]  (l-3cos2e12) 

C  =  (-3/2)[I+l  IZ2  +  1+2  Izi]-sin0i2  cos0i2  exp(-i  <j>i2) 
D  =  (-3/2) [L !  IZ2  +  1-2  Izi]*sin0i2  cos0i2  exp(i  <J>i2) 

E  =  (-3/4)(I+i  42 )  sin20i2  exp(-i  2(J)i2) 

F  =  (-3/4)(I_i  I.2)  sin20i2  exp(i  2§n) 


Expanding  the  dipole-dipole  Hamiltonian  into  this  form  helps  to  illustrate 
several  important  points  about  it.  Firstly,  the  A  and  B  terms,  called  the  secular 
terms  of  the  Hamiltonian,  do  not  involve  net  upward  or  downward  energy 
transitions  and  so  do  not  represent  longitudinal  relaxation  effects.  All  of  the  terms, 
however,  demonstrate  a  spatial  dependence  of  the  interaction  strength.  Term  B 
represents  an  exchange  of  energy  between  two  spins,  and  thus  together  they 
contribute  to  transverse  relaxation  effects.  The  C  and  D  terms,  on  the  other  hand, 
involve  energy  state  transitions  of  one  spin,  with  the  spins  exchanging  energy  at  the 
frequency  C0o  with  the  thermal  lattice.  Similarly,  the  E  and  F  terms  involve  energy 
state  transitions  of  both  spins  in  the  same  direction  and  energy  at  the  frequency  2C0o 
is  either  gained  or  lost  to  the  thermal  lattice  by  the  coupled-spin  pair.  As  a  result, 
spin  state  transitions,  and  therefore  longitudinal  relaxation,  occurs  only  as  a  result  of 
spins  absorbing  or  emitting  energy  at  specific  frequencies,  C0o  and  2(0o.  It  is 
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important  to  note,  however,  that  because  the  heat  capacity  of  the  thermal  lattice  is 
very  large  relative  to  that  of  the  spin  system,  the  lattice  is  essentially  always  at 
thermal  equilibrium  with  its  lower  Boltzmann  energy  states  more  populated  than  the 
higher  energy  states.  As  a  result,  the  spin  system  is  more  likely  to  donate  energy  to 
the  lattice,  than  it  is  to  absorb  energy  from  the  lattice. 

Now  that  the  form  of  the  spin  wave  function  and  of  the  dipole-dipole 
Hamiltonian  have  been  established,  we  can  investigate  the  probability  that  a 
transition  between  spin  states  and  a  change  in  the  energy  of  the  spin  system  will 
occur.  Initially,  we  can  assume  that  a  given  spin  is  in  a  state,  j.  The  probability  that 
this  spin  can  be  found  in  a  different  state,  k,  at  a  time  t  later,  is  given  by  Pjk  (25) 

(1.61)  Pjk  =  Ck*(t)  Ck(t) 

where  the  functions,  Ck(t),  are  as  described  by  equation  1 .52.  Thus,  to  describe  the 
spin  system  relaxation  toward  thermal  equilibrium  we  must  determine  the  functional 
form  of  the  coefficients,  Ck,  which  were  used  to  describe  the  spin  wave  function  'F. 

To  obtain  information  about  the  coefficients  of  the  spin  wave  function,  we 
can  make  use  of  the  time  dependent  Schroedinger  equation 

ST* 

(1.62)  y/tot  T' =  ili-gp 
where  y/tot  is  the  total  Hamiltonian  given  by 


(1.63) 


^tot  —  +  -^d-sec  +  ^d-ns 


and  where  yfd-sec  is  the  secular  part  of  the  dipole-dipole  Hamiltonian  which  contains 
only  the  A  and  B  terms,  and  y/h_ns  is  the  remaining  portion  of  the  dipole-dipole 
Hamiltonian  containing  the  C  to  F  terms.  It  is  important  to  note  that  the  spin 
eigenstates  described  by  equations  1 .46  to  1 .49  are  also  stationary  eigenstates  of  the 
truncated  Hamiltonian,  y/Trunc>  which  is  defined  by 

(1 .64)  -^Trunc  —  -^z  +  -^d-sec 

Because  the  magnetic  field  experienced  by  a  spin  due  to  a  nearby  spin  magnetic 
moment  is  much  smaller  than  the  static  B0  field  provided  by  the  NMR  hardware, 
we  can  make  the  assumption  that 

(1.65)  yfz+#d  -sec  »#d  -ns 

and  make  use  of  time-dependent  perturbation  theory  to  investigate  the  effects  of  the 
dipole-dipole  interaction  on  the  spin  wave  function  (25). 

Time-dependent  perturbation  theory  requires  the  basic  assumption  that  the 
perturbation  Hamiltonian,  y/d-ns,  is  much  smaller  than  the  non-perturbed 
Hamiltonian,  j^Trunc,  as  we  have  already  described,  so  that  the  spin  eigenstates  of 
y/Trunc  can  be  assumed  to  be  reasonable  approximations  for  the  spin  eigenstates  of 
y/tot  as  well.  Making  this  assumption,  and  expanding  yftot  as  described  in  equation 
1 .63,  the  time  dependent  Schroedinger  equation  (equation  1 .62)  becomes 
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The  time  derivative  of 'T  can  also  be  expanded  by  substituting  from  equation  1.51 
as  follows. 


(1.67) 


W  V  5an(t) 


~5T 


n=-l 


5t  Vn  + 


5^,, 

6t 


Because  \J/n  represents  the  stationary  eigenstates  of  'Mz  +  yfd.sec,  the  time  dependent 
Schroedinger  equation  gives  the  result 

( 1  *68)  (7fz  +  ^d-sec)  Xan(^  ¥n  =  i  h  /  van(t) 

n=-l  n=-T 


This  equality  allows  these  terms  to  be  subtracted  from  both  sides  of  equation  1.66 
to  give 

l  1 

(1-^9)  7/(j-ns  ^an(t)  \J/n  =  i  h  ^  ",  tyn 

n=-l  n=-l 

Now,  to  investigate  a  single  Cn,  say  Ck,  we  can  multiply  the  left  side  of  all  terms  in 
equation  1.69  by  ak*(t)\|/k*  and  integrate  over  all  space  to  make  use  of  the 
orthonormal  properties  of  \j/k.  Secondly,  we  can  make  a  simplifying  assumption 
about  the  starting  conditions  by  supposing  that  at  t  =  0,  the  system  is  definitely  in 
state  j,  so  that  Cj(0)  =  1,  and  all  other  Cn(0)  =  0.  The  result  of  these  manipulations 
is  that  equation  1 .69  becomes 


<  \|/k  I  ^cJ-ns  I  Vj  ->  exP 


-i  Ej  r 
h 


•  k  Sck 

1  h_5TexPl 


■i  Ek  r 


(1.70) 
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which,  when  rearranged,  gives  an  expression  for^^ 


(1.71) 


Integrating  this  expression  over  time  yields  an  expression  for  Cp(t)  as  well. 


Ck(t)  =  T-  f  expT -1  (Ei  -  Ek)1']  <  Vk  I  x  I  v  >  dt- 
h  J  h  J 


(1.72) 


u 


Equations  1.71  and  1.72  demonstrate  that  as  a  result  of  the  non-secular  terms  of  the 


dipole-dipole  Hamiltonian,  the  coefficients  used  to  describe  the  two-spin  wave 
function,  are  functions  of  time.  Furthermore,  the  important  energy  difference  in  the 
time  dependence  of  Ck  is  Ej  -  Ek,  the  energy  of  the  spin  state  transition. 

Returning  to  the  problem  of  computing  the  probability  of  a  spin  state 
transition  taking  place,  we  can  now  compute  the  transition  probability  per  unit  time, 
Wjk,  given  by 


(1.73) 


where  Pjk  is  defined  in  equation  1 .61 .  Because  Ck(t)  is  a  function  of  time  only,  we 
can  substitute  the  partial  derivative  of  Ck(t)  (equation  1.71)  for  dCk(t)/dt. 
Substituting  the  expression  given  in  equation  1.72  for  Ck(t)  then  yields  the 


expression 
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o 


(1-74) 


where  CDjk  has  been  used  to  express  (Ej  -  Ek)/h  and  T  =  (t1  -  t).  Equation  1 .74  can 
be  simplified  somewhat  by  recognizing  that  the  second  integral  is  the  complex 
conjugate  of  the  first,  and  that  the  complex  conjugate  in  this  expression  amounts  to 
a  change  in  the  sign  of  T.  Therefore,  we  can  integrate  over  the  range  of  -t  to  t  and 
combine  the  two  integrals  as  follows 


(1.75) 


Furthermore,  if  the  integration  interval,  2t,  is  much  larger  than  the  period  27t/COjk  we 
can  change  the  integral  limits  to  be  -oo  to  oo  to  yield  equation  1.76.  The  usefulness 
of  such  a  change  will  become  apparent  in  the  following  discussion. 


OO 


Wjk  —  _  T/j-ns  7/d-u.s  ^k^xp[-iCOjkT]dT 

n  -oo 


(1.76) 


-00 


The  non-secular  part  of  the  dipole-dipole  Hamiltonian,  #d-ns>  can  be 
expanded  into  its  terms,  C  to  F,  and  each  of  these  can  be  expressed  as  a  product  of  a 
space  term,  Z,  and  a  spin  term,  O. 


Tfd-ns  -  ZnOn 

n  =  Cto  F 


(1.77) 
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With  Tfd-ns  expanded  into  this  form,  the  terms  within  the  integral  in  equation  1.76 
become 

(1-78)  <\\fk I  7/d.ns  l\j/j>  =  V  Zn  <\|/kl  On  l\J/j> 

n  =  tT  to  F 

(1.79)  7/d-ns  =  ^  Zn  <\J/jl  On  l\J/k> 

n  =  C  to  F 


and  equation  1.74  becomes 

00 

(1 .80)  wjk  =  X  l<Vkl  IVj>l2  Jz(t)Z(t')  exp[-i(OjkT]dT 

n  n  =  C  to  F  -oo 

Because  the  relative  positions  and  orientations  of  the  two  spins  vary  randomly  in 
time,  the  space  parts  of  the  Hamiltonian,  Z(t),  vary  randomly  in  time.  As  a  result, 
Z(t)  is  not  an  analytical  function  so  we  cannot  exactly  predict  its  value  at  every  point 
in  time. 


In  order  to  deal  with  the  random  function,  Z(t),  we  can  compute  an 
autocorrelation  function,  G(T)  which  describes  how  well  Z(t)  correlates  with  itself 
over  a  time  interval  of  duration  T.  This  autocorrelation  function  is  defined  by 

(1.81)  G(T)=  Z(t)Z(t') 


where  denotes  the  ensemble  average,  and  again  T  =  (t1  - 1).  To  describe  the 

distribution  of  the  random  motions  of  the  interacting  spins  it  is  assumed  that  the 
auto-correlation  function  decays  exponentially  with  the  characteristic  time,  Tc, 
known  as  the  correlation  time,  so  that  G(T)  has  the  form 
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(1.82)  G(T)  =  G(0)  exp(-T/Tc) 

To  determine  the  frequency  distribution  of  the  random  spin  motions,  and  therefore 
to  determine  at  which  motional  frequencies  the  thermal  energy  of  the  spins  is  stored, 
we  compute  the  Fourier  transform  of  the  autocorrelation  function 

+  00 

(1.83)  J(CO)  =  J  G(T)  exp(-i  COT)  dT 

-oo 


The  function  J(CD)  is  known  as  the  spectral  density.  Substituting  the  exponential 
form  of  G(T)  given  in  equation  1.82  yields  the  Debye  spectral  density,  where 

+  00 

(1.84)  J(CO)  =  J  G(0)  exp(-T/Tc)  exp(-i  COT)  dT 

-00 

( 1 .85)  J(CO)  =  G(0) - ^ - - 

[l+(C0Xc)2] 

The  Debye  spectral  density  describes  the  magnitudes  of  the  frequency  components 
of  the  random  thermal  motions  of  nuclei  in  a  liquid-like  spin  system,  as  shown  in 
figure  1.5. 
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Figure  1.5  The  dependence  of  the  normalized  Debye  spectral  density  on  frequency, 
expressed  relative  to  the  correlation  time,  Tc. 

This  figure  demonstrates  that  the  random  motions  of  the  spins  may  have  many 
frequency  components,  two  of  which  correspond  to  the  frequencies  C0o  and  2C0o. 
Furthermore,  it  can  be  seen  that  by  computing  the  ensemble  average  of  Wjk  in 
equation  1.80,  the  integral  term  in  this  equation  can  be  substituted  by  the  Debye 
spectral  density  described  in  equation  1.83.  As  a  result 

(1 .86)  Wik  =  A  Z  kVkl  <i>n  IVj>'2  J(®jk) 

h2  n=TtoF 

Since  we  are  interested  in  relaxation  of  the  spin  system  toward  thermal  equilibrium, 
and  the  thermal  lattice  is  more  likely  to  receive  energy  from  the  spin  system  than  it 
is  to  give  up  energy  to  the  spin  system,  we  are  concerned  only  with  those  terms  of 


^d-ns  resulting  in  transitions  to  lower  energy  states,  i.e.  the  C  and  E  terms.  The 
longitudinal  relaxation  rate,  Rj,  is  therefore  given  by 
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(1.87) 


Substituting  for  the  form  of  Wj^  from  equation  1 .86  we  obtain  the  expression 


(1.88) 


Ri  <*  ki  J(C0o)  +  ki'  J(2C0o) 


where  ki  and  k i '  are  constants  which  allow  for  the  different  magnitudes  of  the 
space  parts  of  the  C  and  E  terms  of  ^fd-ns-  For  the  more  specific  case  of  a  system  of 
weakly  interacting  protons,  this  expression  is  (74) 


(1.89) 


where  ki"  is  a  constant  which  depends  on  the  dipole-dipole  interaction  strength,  the 
nature  of  the  motion,  and  nuclear  constants  such  that 


(1.90) 


This  equation  describes  the  spin-lattice  relaxation  rate  dependencies  on  both  the 
resonant  frequency  of  the  system  as  well  as  the  correlation  times  of  the  motions  of 
the  interacting  spins.  The  relaxation  mechanisms  are  most  efficient,  and  Ri  is  at  a 
maximum,  when  00oTc  =  0.64  (i.e.  approximately  unity)  as  demonstrated  in  figure 
1.6.  The  NMR  hardware  employed  for  the  studies  discussed  in  Chapters  2  and  3 
resulted  in  a  proton  resonant  frequency,  C0o,  of  100  MHz  (108  Hz)  and  for  water 
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molecules  in  solution  at  room  temperature  the  value  of  Tc  is  roughly  3  x  10' 12 
seconds  (30).  The  product  of  these  two  constants,  C0oTc,  has  a  value  of  10~4  and 
figure  1.6  therefore  demonstrates  that  the  relaxation  of  pure  water  protons  is  not 
very  efficient  at  a  proton  resonant  frequency  of  100  MHz.  A  macromolecule  of 
molecular  weight  20  000,  however,  may  have  a  rotational  correlation  time  of  10'8 
seconds  and  relaxation  rates  may  be  nearly  maximal  for  protons  in  such  molecules 
at  a  resonant  frequency  of  100  MHz  (30). 


Figure  1.6  The  dependence  of  the  quantity,  C0o[J(C0o)  +  4  J(2C0o)],  and  therefore 
the  dependence  of  Ri,  on  the  value  ot  Tc,  for  a  fixed  frequency,  C0o. 

The  time  dependent  nature  of  the  coefficients,  Ck,  arising  from  the  non- 
secular  terms  of  the  dipole-dipole  Hamiltonian,  results  in  a  time  dependent 
expectation  value  for  the  individual  transverse  spin  components  (equation  1.19).  In 
a  description  of  the  relaxation  of  the  transverse  component  of  the  magnetization  it  is 
therefore  necessary  to  include  the  non-secular  terms  of  the  dipole-dipole 
Hamiltonian  as  well  as  the  secular  (A  and  B)  terms.  As  a  result,  the  transverse 
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relaxation  rate  also  depends  on  those  terms  used  to  describe  the  longitudinal 
relaxation  rate,  namely  J(C0o)  and  J(2(0o).  In  addition,  the  A  and  B  terms  of  the 
dipole-dipole  Hamiltonian  involve  near-zero  frequency  components  of  the  spin 
motions  so  that  the  transverse  relaxation  rate,  R2,  is  given  by  (1) 

(1-91)  R2  k2  J(0)  +  k2'  J(C0o)  +  k2"  J(2C0o) 

where  k2,  k2'  and  k2"  are  constants,  similar  to  those  employed  in  equation  1.88. 
The  transverse  relaxation  rate  of  a  system  of  weakly  interacting  protons  is  given  by 

(1 .92)  R2  =  k2"  {  3  J(0)  +  5  J(C0o)  +  2  J(2C0o)} 


where  k2"  is  a  constant  which  depends  on  the  dipole-dipole  interaction  strength  and 
nuclear  constants  (1)  such  that 


(1.93) 


* = (id 


'y4  fo2‘ 

nr. 


3XC  + 


5Tr 


+ 


2Tr 


[l+(C0oTc)2]  [1  +4(C0oXc)2]J 


This  expression  for  the  transverse  relaxation  rate,  in  conjunction  with  figure  1 .6  and 
equation  1.90,  demonstrates  that  for  a  system  of  protons  having  a  resonant 
frequency  of  100  MHz  and  a  correlation  time  of  10" 12  seconds,  Ri  and  R2  are 
approximately  equal.  However,  at  the  same  proton  resonant  frequency,  for  a  spin 
system  having  a  correlation  time  of  10'8  seconds  the  value  of  R2  is  approximately 
twice  that  of  Ri.  Furthermore,  because  the  transverse  relaxation  rate  depends  on 
J(0)  it  does  not  reach  an  upper  limit  at  aXc  value  near  l/C0o  as  does  Ri,  but  rather 
continues  to  increase  with  increasing  Tc.  As  a  result,  it  is  possible  for  R2  to  be 
orders  of  magnitude  greater  than  R]. 
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Quadrupolar  Relaxation 

Quadrupolar  relaxation  occurs  only  when  a  nucleus  with  a  spin  I  >  1/2,  and 
therefore  possessing  an  electric  quadrupole  moment,  interacts  with  a  time  dependent 
local  electric  field  gradient.  The  source  of  this  electric  field  gradient  is  usually  the 
electronic  structure  of  the  molecule  which  contains  the  nucleus  of  interest,  but 
nearby  ions  or  molecules  can  also  contribute.  For  the  case  in  which  a  nucleus  with 
spin,  I,  and  electric  quadrupole  moment,  Q,  interacts  with  a  field  described  by  an 
electric  potential,  V(x,y,z),  the  Hamiltonian  of  the  interaction,  yfq,  is  given  by  (72) 

0-9 4)  .Tfq  —  4I(2L-1)  +  B'  +  O  +  D'  +  E'  ] 

where:  A'  =  Vzz  (3IZ2  - 12) 

B'  =  (Vzx  +  i  Vzy)  (LIZ  +  IZL) 

C'  =  (Vzx  -  i  Vzy)  (I+Iz  +  IZI+ ) 

D'  =  [j(Vxx  -  Vyy)  +  i  Vxy]  (I.)2 

E'  =  [j(Vxx  -  Vyy)  -  i  Vxy]  (I+)2 

The  quantity  e  is  the  electronic  charge,  and  Vxx,  Vxy,  Vyz,  ...  etc.,  represent  the 
second  partial  derivatives  of  V  with  respect  to  x,  y,  and  z.  It  is  important  to  keep  in 
mind  that  the  electric  potential,  V,  is  determined  by  the  electronic  structure  of  the 
molecule  and  that  the  derivatives  of  this  potential  can  vary  randomly  at  the  location 
of  the  nucleus  of  interest.  Thus,  with  7/q  expressed  in  this  form,  and  recognizing 
the  fact  that  each  of  the  terms  A1  to  E1  contain  both  a  space  term  and  a  spin  term,  the 
similarities  between  7fq  and  #d-d  in  equation  1.60  are  apparent.  As  a  result,  the 
method  that  was  used  to  describe  dipole-dipole  relaxation  from  the  starting  point  of 
the  dipole-dipole  Hamiltonian,  can  also  be  applied  to  the  quadrupolar  Hamiltonian 


to  describe  quadrupolar  relaxation.  Expressions  for  the  transverse  and  longitudinal 
quadrupolar  relaxation  rates,  of  an  ensemble  of  spin  1  =  1  nuclei  in  a  liquid  like 
system,  are  given  in  equations  1.95  and  1.96,  respectively  (1). 
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(1-95) 


Ri  =  (is)  5 


+ 


4Xr 


[  1  +  (COoTc)2]  [1  +4(C0oTc)2]. 


(1.96) 


R2_(ioo)  ^ 


9Tc  + 


1 5Tr 


+ 


6  Tc 


[1  +(C0oTc)2]  [1  +4(C0oTc)2]J 


where  the  correlation  time,  Tc,  characterizes  the  time  dependence  of  the  spatial 
derivatives  of  V  at  the  location  of  the  nucleus  of  interest.  The  coefficient,  £,, 
expressed  in  a  principal  axis  coordinate  system  (x',y',z',  such  that  Vjj  =  0  for  i  *  j)  is 
given  by 


(1.97) 


In  this  expression,  r|,  known  as  the  asymmetry  parameter,  is  a  measure  of  the 
electric  field  gradient  deviation  from  axial  symmetry,  and  is  given  by 


(1.98) 


T)  =' 


V 


x'x' 


V 


Vy'y1 


z'z' 


The  rates  of  nuclear  relaxation  resulting  from  the  dipole-dipole  and  quadrupolar 
interactions  have  similar  dependencies  on  C0o  and  Tc  because  we  have  assumed  the 
same  form  of  the  autocorrelation  function  to  characterize  the  random  motions  which 
mediate  each  of  these  interactions.  However,  the  strengths  of  these  two  interactions 
can  be  considerably  different.  This  is  demonstrated  by  the  fact  that  the  Ti  of  water 


protons  in  a  cat  brain  is  dominated  by  the  dipole-dipole  interaction  and  has  a  value 
of  1.5  seconds.  In  a  hypothetical  situation  in  which  we  exchanged  these  water 
protons  for  deuterons,  and  in  which  the  dipole-dipole  interaction  was  still  dominant, 
the  Ti  value  resulting  from  equation  1.90  would  be  approximately  2700  seconds 
because  of  the  lower  deuteron  gyromagnetic  ratio  (Yproton/Ydeuteron  =  6.5).  The  Ti 
of  deuterons  in  D2O  in  a  cat  brain,  however,  is  dominated  by  quadrupolar 
relaxation  and  has  a  value  of  roughly  250  msec  (26).  The  ratio  of  the  interaction 
strengths,  i.e.  ^/(Y4h2/r6),  in  this  hypothetical  situation,  is  8.6  x  104.  Thus,  in 
comparing  these  two  specific  examples  it  can  be  seen  that  the  strength  of  the 
quadrupolar  interaction  can  be  several  orders  of  magnitude  greater  than  that  of  the 
dipole-dipole  interaction. 

The  Two-Site  Rapid  Exchange  Model  of  Relaxation 

In  biological  systems,  such  as  the  lung,  the  environment  affecting 
magnetization  relaxation  is  often  spatially  heterogeneous.  In  cases  in  which  there 
are  two  relaxation  environments,  the  relaxation  rates  observed  depend  strongly  on 
the  rate  at  which  mobile  nuclei  can  exchange  between  the  two  environments  (29). 
The  description  of  the  relaxation  in  a  system  with  two  relaxation  environments  is 
generally  tractable  if  it  can  be  expressed  in  either  the  "fast-exchange"  case,  in  which 
the  mean  exchange  time  between  the  two  environments  is  much  shorter  than  the 
relaxation  time  in  either  environment,  or  the  "slow-exchange"  case,  in  which  the 
relaxation  times  are  much  shorter  than  the  exchange  time  between  the  environments 
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Hydration  Sphere  Water  Bulk  Water 

Figure  1.7  Schematic  representation  of  water  protons  in  two  rapidly  exchanging 
relaxation  environments;  bulk  water  and  water  in  a  hydration  sphere  around  a 
macromolecule. 

The  "fast-exchange"  case  generally  occurs  when  two  distinct  relaxation 
environments  are  either  in  contact,  or  are  separated  by  a  physical  boundary  which 
only  marginally  impedes  the  diffusion  of  the  water  molecules  from  one  environment 
to  the  other.  An  example  of  this  case  is  illustrated  in  figure  1 .7  in  which  protons  in 
bulk  water  can  diffuse  rapidly  from  the  slowly  relaxing  bulk  water  environment  to 
the  hydration  sphere  where  they  can  exchange  spin  energy  with  the  much  faster 
relaxing  protons  in  the  macromolecule  in  solution  (29).  Similarly,  this  situation 
may  arise  for  fluid  in  the  airspaces  of  the  lung,  which  is  also  in  contact  with 
epithelial  cell  membranes.  Macromolecules  which  are  incorporated  into  a  cell 
membrane  also  contain  relatively  fast  relaxing  protons  which  can  exchange  spin 
energy  with  water  protons  in  the  hydration  sphere.  Because  a  given  water  proton 
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can  spend  time  in  the  bulk  water  and  in  a  hydration  sphere  while  relaxing  to  thermal 
equilibrium  with  the  lattice,  we  can  observe  only  an  average  relaxation  rate,  Robs> 
for  the  spin  system.  The  relaxation  rate  observed  is  a  weighted  average  of  the 
relaxation  rates  in  each  of  these  two  environments,  as  shown  in  equation  1 .99. 

(1 .99)  Robs  —  Pfree  Rfree  +  Pbound  Rbound 

The  relaxation  rates  Rfree  and  Rbound  are  those  of  water  protons  in  the  bulk  water 
and  in  the  hydration  sphere  water,  respectively.  The  quantities  pfree  and  pbound  are 
the  relative  fractions  of  water  protons,  contributing  to  the  water  signal,  in  the  bulk 
water  and  in  the  hydration  sphere,  respectively. 

The  "slow-exchange"  case,  on  the  other  hand,  gives  rise  to  two  or  more 
relaxation  components,  and  generally  occurs  when  two  different  relaxation 
environments  are  separated  by  a  physical  boundary  which  impedes  the  diffusion  of 
water  molecules  between  the  two  environments.  An  example  of  this  case  is  lung 
alveolar  fluid  with  a  water  proton  relaxation  time  of  Tajv,  which  is  separated  from 
interstitial  fluid  with  a  water  proton  relaxation  time  of  Tjnt,  by  the  highly 
impermeable  alveolar  epithelium.  Because  the  epithelium  prevents  the  rapid 
diffusion  of  water  between  alveolar  and  interstitial  fluid,  a  given  water  proton  will 
relax  under  the  influence  of  only  one  of  the  two  environments.  The  NMR  signal 
from  such  a  spin  system  is  thus  a  sum  of  the  signals  arising  from  water  protons  in 
each  of  the  two  relaxation  environments,  and  demonstrates  a  bi-exponential  decay 
as  indicated  by  equation  1.100. 


(1.100) 


S(t)  —  S0  [Palv  exP("l/Ta]v)  +  pint  exp(-t/Tint)] 
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In  this  equation,  the  quantities  pajv  and  Pint  are  the  relative  fractions  of  water 
protons,  contributing  to  the  water  signal,  in  the  alveolar  fluid  and  in  the  interstitial 
fluid,  respectively. 

The  studies  discussed  in  the  following  chapters  include  relaxation  time 
measurements  for  lung  water  protons.  Because  of  the  heterogeneous  nature  of  the 
lung,  the  measured  lung  water  proton  relaxation  times  are  highly  dependent  on  the 
two-site  rapid  exchange  model  of  relaxation.  Moreover,  it  is  a  consequence  of  the 
bi-exponential  signal  decay  detected  in  the  "slow-exchange"  case  that  we  believe  we 
are  able  to  measure  relaxation  times  of  two  distinct  fluid  compartments  in  the  lung. 


1.3.3  NMR  Imaging  Theory 

In  order  to  obtain  a  two  dimensional  (2D)  NMR  image  of  a  given  plane  in 
space  (that  presumably  transects  an  object  of  interest),  a  nuclear  excitation  pulse 
sequence  must  satisfy  three  fundamental  requirements  which  are; 

1)  to  generate  a  detectable  NMR  signal, 

2)  to  provide  a  means  of  defining  the  desired  plane  with  spatially  selective 
pulses,  and 

3)  to  provide  a  means  of  spatially  encoding  the  signal  in  two  dimensions  in 
that  plane. 

The  following  sections  detail  the  production  of  spin-echoes  and  gradient-recalled 
echoes,  spatially  selective  pulses,  the  spatial  encoding  of  magnetization  represented 
in  k-space,  and  some  common  image  reconstruction  techniques. 
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Theory  of  Spin  Echoes  and  Gradient-Recalled  Echoes 

The  presence  of  local  inhomogeneities  in  a  static  B0  field,  whether  caused  by 
variations  in  magnetic  susceptibility  at  regional  boundaries  in  a  heterogeneous 
sample,  or  by  the  NMR  hardware,  causes  difficulties  in  the  observation  of  the 
natural  transverse  relaxation  time  of  a  spin  system,  T2,  by  reducing  the  transverse 
relaxation  time  observed  from  the  natural  T2,  i.e.  T2*  <  T2.  To  overcome  these 
difficulties  we  make  use  of  the  fact  that  dephasing  of  the  transverse  magnetization 
caused  by  static  field  inhomogeneities  is,  in  many  cases,  reversible.  Figure  1.8 
illustrates  a  situation  in  which  the  net  magnetization  of  a  spin  ensemble  has  been 
rotated  onto  the  yp  axis  by  means  of  a  90°  RF  pulse.  The  ensemble  magnetization 
is  then  left  to  evolve  under  the  influence  of  B0  for  a  time  period  T,  during  which 
static  field  inhomogeneities  cause  the  transverse  magnetization  components  to  lose 
phase  coherence. 


-T  /  T 

M  =  M0  e  2 


Figure  1.8  Dephasing  of  transverse  magnetization  components  due  to  static 
magnetic  field  inhomogeneities,  depicted  in  the  rotating  frame  of  reference. 
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180x  RF  Pulse 


Figure  1.9  The  time  course  of  changes  in  the  transverse  magnetization  of  an 
ensemble  of  spins  due  to  the  effects  of  static  B0  inhomogeneities  being  refocussed 
at  a  time  2T  after  the  initial  90x  RF  pulse  rotated  the  magnetization  into  the 
transverse  plane. 


If  at  the  time  t  =  T,  as  defined  in  figure  1.8,  a  180°  RF  pulse  is  applied  with  B]  in 
the  Xp  direction  (written  180x),  the  spins  are  flipped  about  the  Xp  axis  and  the  phase 
distribution  about  the  -yp  axis  becomes  the  opposite  of  that  phase  distribution  about 
the  +yp  axis  immediately  prior  to  the  application  of  the  RF  pulse.  In  other  words,  a 
spin  which  precesses  faster  than  the  rotating  reference  frame,  such  that  it  is  at  an 
angle  (j)f  from  the  yp  axis  at  the  time  t  =  T,  is  rotated  by  means  of  the  180x  pulse  to 
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be  at  an  angle  -(J)f  from  the  -yp  axis.  Similarly,  a  spin  which  precesses  slower  than 
the  reference  frame  and  therefore  lags  the  +yp  axis  by  an  angle  <J)S  at  the  time  t  =  X, 
is  rotated  by  means  of  the  180x  pulse  to  lead  the  -yp  axis  by  the  angle  <J)S  .  Thus, 
during  a  further  time  interval  of  duration  X,  all  spin  dephasing  due  to  static  B0  field 
inhomogeneities  is  reversed  and  the  net  magnetization  in  the  transverse  plane 
reaches  a  maximum  after  this  second  X  and  is  then  directed  along  the  -yp  axis.  The 
magnitude  of  the  refocussed  transverse  magnetization  is  reduced  from  its  initial 
value  at  t  =  0,  however,  by  the  irreversible  transverse  relaxation  as  illustrated  in 
figure  1.9.  After  the  transverse  magnetization  has  been  refocussed,  the  evolution  of 
the  relative  phases  of  the  spins  continues,  and  once  again  the  spins  fan-out  in  the 
transverse  plane.  This  sequence  of  RF  pulses  is  known  as  the  Hahn  spin-echo 
sequence  (34)  and  is  written  symbolically  as  [  90x  -  X  -  180x  -  X  ].  The  NMR 
signal  that  can  be  detected  from  an  ensemble  of  spins  undergoing  this  process  is 
illustrated  in  figure  1.10. 


inhomogeneities  being  refocussed  at  a  time  2X  after  the  initial  90°  RF  pulse  rotated 
the  magnetization  into  the  transverse  plane. 
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Note  in  figure  1.10  that  the  signal  has  been  drawn  as  if  the  receiver  phase  was 
shifted  180°  after  detecting  the  free-induction  decay  signal  (FID),  in  anticipation  of 
the  180°  change  in  phase  of  the  echo  signal,  and  the  receiver  reference  frequency  is 
slightly  offset  from  the  spin  precessional  frequency.  Although  the  FID  signal 
decays  with  the  characteristic  time  T2*,  the  amplitude  of  the  refocussed  "echo" 
signal  depends  only  on  the  natural  transverse  relaxation  of  the  spin  system  and 
decays  with  the  transverse  relaxation  time,T2-  Furthermore,  with  the  application  of 
a  second  180°  RF  pulse  at  the  time  3T,  a  second  echo  can  be  formed  at  the  time  4T 
and  so  on,  for  as  long  as  the  echoes  have  non-zero  amplitudes.  Such  an  expansion 
of  the  Hahn  sequence  is  known  as  the  Carr-Purcell  (CP)  sequence  (17)  and  is 
written  [  90x  -  T  -  (180x  -  T)n  ].  By  fitting  an  exponentially  decaying  curve  to  the 
echo  amplitude  maxima,  S(2nX),  where  n  is  an  integer,  the  values  of  T2  and  SQ  can 
be  determined  as  defined  by  equation  1 .101. 

(1.101)  S(2nT)  =  S0  exp(-2nT/T2) 

Because  the  Hahn  and  CP  spin-echo  amplitudes  are  highly  sensitive  to  spatial 
variations  or  setting  errors  in  the  180°  refocussing  pulses,  the  Carr-Purcell- 
Meiboom-Gill  (CPMG)  (55)  spin-echo  sequence  [  90x  -  T  -  (180y  -  T)n  ],  was 
developed.  Spin-echoes  produced  with  the  CPMG  sequence  are  in  phase  because 
the  magnetization  is  initially  rotated  onto  the  yp  axis  and  the  180°  refocussing 
pulses  rotate  the  magnetization  about  this  same  axis.  Also,  if  the  180°  pulses  are 
imperfect  and  actually  rotate  the  magnetization  through  slightly  more,  or  less,  than 
180°,  the  error  in  the  first  pulse  is  negated  by  the  error  in  the  second  pulse,  and  so 
on  through  the  sequence.  As  a  result,  the  even  numbered  echoes  produced  by  the 
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CPMG  sequence  may  provide  more  accurate  data  for  measurements  of  T2  values, 
than  those  provided  by  either  the  Hahn  or  the  CP  spin-echo  sequences. 

An  alternative  to  generating  an  echo  signal  by  means  of  a  180°  refocussing 
pulse  is  to  reverse  a  gradient  and  produce  what  is  called  a  gradient-recalled  echo. 
Such  an  echo  signal  can  be  formed  by  first  applying  a  magnetic  field  gradient  to  an 
ensemble  of  spins  which  are  initially  in  phase,  so  that  the  spins  are  forced  into  a 
linear  phase  distribution  along  the  gradient  direction,  and  the  net  transverse 
magnetization  reduces  to  zero.  The  phase  distribution  produced  by  this  gradient  can 
then  be  reversed  by  applying  an  equal  and  opposite  gradient,  the  net  transverse 
magnetization  again  approaches  a  maximum  as  an  echo  signal  is  formed  which  is 
similar  to  that  produced  by  a  spin-echo  sequence.  The  use  of  a  magnetic  field 
gradient  reversal  to  reverse  the  spin  phase  distribution  has  no  effect  on  the 
dephasing  caused  by  random  inter-nuclear  interactions  or  by  static  magnetic  field 
inhomogeneities.  The  maximum  amplitude  of  a  gradient-recalled  echo  therefore 
decays  exponentially  with  the  characteristic  time  T2*,  a  salient  difference  from  the 
Hahn  spin-echo.  Because  the  formation  of  a  gradient-recalled  echo  does  not  require 
the  transmission  of  tailored  RF  pulses,  most  NMR  systems,  particularly  if  they 
have  actively  shielded  gradient  coils,  are  able  to  produce  gradient-recalled  echoes 
with  a  shorter  inter-echo  interval  than  is  possible  for  spin  echoes.  It  is  the  ability  to 
produce  echo  signals  in  rapid  succession  that  makes  the  gradient-recalled  echo 
sequence  preferable  to  the  spin-echo  sequence  for  some  imaging  applications  (57). 

Spatially  Selective  RF  Pulses 

The  interaction  between  an  ensemble  of  nuclear  spins  in  a  static  B0  field 
with  a  rotating  RF  Bi  field,  as  discussed  in  section  1.3.1,  is  highly  dependent  on  the 
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frequency  of  the  Bi  field  relative  to  the  resonant  frequency  of  the  spins  in  the 
ensemble.  It  is  the  frequency  dependence  of  this  interaction  that  provides  us  with 
the  means  to  excite  spins  selectively  at  a  chosen  location  in  space.  This  spatial 
selectivity  requires  the  resonant  frequency  of  the  spins  to  be  spatially  dependent  and 
so,  for  most  NMR  imaging  applications,  a  linear  magnetic  field  gradient  is 
superimposed  onto  the  homogeneous  B0  field.  The  resultant  static  field  is  then  as 
described  by  equation  1.102  for  a  gradient  in  the  z  direction. 


B(z)  =  B0  k  +  (Gz  z)k 


(1.102) 


In  this  expression  the  quantity  Gz  =  8b/5z  and  b  =  b  k  is  the  gradient  field.  As  a 
result  of  this  magnetic  field  gradient,  the  resonant  frequency  varies  along  the  z 
coordinate  axis  as  depicted  in  equation  1.103 


C0(z)  =  C0o  +  Y  Gz  z 


(1.103) 


If  a  Bi  field  at  the  frequency  C0o  is  now  applied  along  the  xp  axis  in  the  rotating 
frame,  the  effective  field  in  the  rotating  frame,  Beff  (defined  by  equation  1 .34),  is  a 
function  of  position  along  the  z  axis,  and  is  given  by 


(1.104) 


Beff(z)  =  (B0  +  Gz  z  +  fi/y)  Kp  +  Bi  ip 


and  because  Q  =  C0o  =  -yB0  (equation  1.32) 


(1.105) 


BefKz)  =  Gz  z  kp  +  Bi  ip 
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If  the  Bi  field  is  applied  for  the  appropriate  duration,  spins  at  z  =  0  can  be  rotated 
through  90°  into  the  yp  axis  direction.  However,  the  spins  at  z  ^  0  will  be  rotated  to 
lie  at  an  angle  of  less  than  90°  from  the  z  axis  by  such  a  Bi  pulse,  and  their 
transverse  components  may  not  be  aligned  with  the  yp  axis,  as  illustrated  in  figure 
1 .2.  For  large  z,  such  that  z  »  Bi/Gz,  the  effective  field  in  the  rotating  frame  is 
essentially  unaffected  by  the  Bi  pulse,  and  the  spins  remain  aligned  with  B0  field. 
Figure  1.11  illustrates  the  transverse  magnetization,  Mxy,  resulting  from  the 
application  of  a  rectangular  90°  RF  pulse  of  duration  2T  and  frequency  C0o,  as 
computed  using  the  Bloch  equations.  The  corresponding  phase  distribution  of  the 
transverse  spin  components  is  illustrated  in  figure  1.12. 


Figure  1.11  The  z  direction  variation  of  the  transverse  magnetization  of  a  uniform 
spin  system,  after  the  application  of  a  rectangular  90°  RF  pulse  of  duration  2T  in  the 
presence  of  a  magnetic  field  gradient  in  the  z  direction.  The  transverse 
magnetization  amplitude  has  been  normalized  to  equal  1  at  z  =  0. 
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Figure  1.12  The  phase  distribution  of  magnetic  moments  in  the  transverse  plane, 
along  the  z  direction,  after  the  application  of  a  rectangular  90°  RF  pulse  of  duration 
2T  in  the  presence  of  a  magnetic  field  gradient  in  the  z  direction. 

The  transverse  magnetization  distribution  resulting  from  a  rectangular  90° 
RF  pulse,  as  demonstrated  in  figure  1.11,  is  obviously  unsuitable  for  imaging 
purposes  because  of  the  multiple  bands  of  excitation.  Moreover,  as  demonstrated 
by  figure  1.12,  the  spins  are  roughly  linearly  distributed  in  phase  so  that  the  net 
transverse  magnetization  of  the  spin  system  approximates  zero.  However, 
Mansfield  et  al.  (48)  have  shown  that  the  response  of  a  spin  system  to  an  RF  pulse 
in  the  presence  of  a  linear  magnetic  field  gradient  is  roughly  equal  to  the  Fourier 
transform  of  the  RF  pulse  envelope.  This  similarity  is  also  demonstrated  by  a 
comparison  between  the  magnetization  distribution  illustrated  in  figure  1.1 1,  and  the 
Fourier  transform  shown  in  figure  1.13,  of  a  rectangular  function  with  amplitude 
1/2  and  width  2T,  and  provides  a  clue  to  the  shape  of  RF  pulse  that  would  be 
required  to  excite  a  single  band  of  spins.  However,  because  the  response  of  a  spin 
system  to  an  RF  pulse  is  non-linear  (i.e.  Mxy  is  not  a  linear  function  of  Bj),  the 
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similarity  between  the  transverse  magnetization  distribution  and  the  Fourier 
transform  of  the  RF  pulse  envelope  is  greatest  for  pulses  which  rotate  the 
magnetization  through  a  small  tip  angle,  ({),  such  that  sintj)  «  ({),  (see  equation  1.40) 
and  gets  progressively  worse  as  the  tip  angle  increases. 


Figure  1.13  The  amplitude  of  the  Fourier  transform  of  a  rectangular  function  with 
amplitude  1/2  and  width  2T. 

The  similarity  between  the  transverse  magnetization  distribution  and  the 
Fourier  transform  of  the  RF  pulse  envelope  leads  one  to  expect  that  an  RF  pulse 
modulated  as  a  sine  function  (i.e.  sinc(x)  =  sin(x)/x,  the  form  of  the  Fourier 
transform  of  a  rectangular  pulse,  as  illustrated  in  figure  1.13)  would  result  in  an 
approximately  rectangular  transverse  magnetization  distribution.  However,  the  long 
duration  of  a  sine  modulated  RF  pulse  (~  30  msec  may  be  required)  makes  it 
unsuitable  for  use  with  the  spin-echo  imaging  of  lung  alveolar  fluid  discussed  in 
Chapters  2  and  3  because  the  faster  relaxing  component  of  lung  water  (T2  ~  25 
msec)  would  undergo  appreciable  transverse  relaxation  during  such  a  pulse. 
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Figure  1.14  The  z  direction  variation  of  the  transverse  magnetization  of  a  spin 
system,  normalized  to  1  at  z=0  (solid  line),  after  the  application  of  a  Gaussian  90° 
RF  pulse  in  the  presence  of  a  z  direction  magnetic  field  gradient.  Also,  shown  for 
comparison  is  a  Gaussian  function  having  the  same  maximum  amplitude,  and  width 
at  half  maximum,  as  that  of  the  transverse  magnetization  distribution  (dashed  line). 


Figure  1.15  The  phase  distribution  of  magnetic  moments  in  the  transverse  plane, 
along  the  z  direction,  after  the  application  of  a  Gaussian  90°  RF  pulse  in  the 
presence  of  a  magnetic  field  gradient  in  the  z  direction. 
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A  more  suitable  choice  is  the  much  shorter  duration  Gaussian  modulated  RF  pulse 
which  results  in  an  approximately  Gaussian  transverse  magnetization  distribution 
along  the  z  direction.  The  actual  response  of  the  spin  system  to  a  90°  Gaussian 
pulse,  calculated  with  the  Bloch  equations,  is  shown  in  figures  1.14  and  1.15  as  the 
transverse  magnetization  distribution,  and  the  phase  distribution  of  the  magnetic 
moments  in  the  transverse  plane,  respectively.  Figure  1.14  also  demonstrates  a 
comparison  between  the  transverse  magnetization  distribution  resulting  from  a  90° 
Gaussian  pulse,  and  a  Gaussian  distribution  having  the  same  maximum  amplitude 
and  width  at  half  maximum.  The  transverse  magnetization  distribution  along  the  z 
axis  now  has  a  single  band  of  excitation  as  desired,  with  an  approximately  Gaussian 
distribution.  The  phase  distribution  of  the  magnetization  is  again  roughly  linear  and 
so  can  be  reversed  with  the  application  of  a  gradient  pulse  of  the  appropriate 
amplitude  and  duration,  after  the  RF  pulse  has  been  applied.  Because  of  these 
desirable  properties  of  the  spin  system  response  to  Gaussian  RF  pulses,  spatially 
selective  90°  pulses  (7.7  msec  duration)  of  this  type  are  employed  in  the  studies 
described  in  chapters  2  and  3. 

A  k-Space  Representation  of  NMR  Imaging  Theory  and  Image 

Reconstruction 

Generating  an  image  with  the  NMR  signal  from  an  object  requires  that  the 
components  of  the  signal,  whether  as  part  of  an  FID  or  an  echo  signal,  have  some 
dependence  on  the  spatial  location  of  their  points  of  origin.  As  discussed  earlier  in 
this  section,  the  precessional  frequency  of  spins  can  be  made  spatially  dependent  by 
applying  a  magnetic  field  gradient.  Also,  the  relative  phase  of  the  spin  magnetic 
moments  can  be  made  spatially  dependent  by  applying  a  magnetic  field  gradient  for 
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a  fixed  duration  during  the  sequence.  In  order  to  describe  the  response  of  the  spin 
system  to  these  space  encoding  mechanisms  we  first  assume  that  the  magnetization 
has  been  rotated  from  its  thermal  equilibrium  direction  onto  the  Xp  axis  by  means  of 
a  slice  defining  90°  RF  pulse  in  the  presence  of  a  magnetic  field  gradient  in  the  z 
direction.  The  magnetization,  m,  of  an  isochromat,  a  small  volume  element  over 
which  it  can  be  assumed  that  all  spins  experience  that  same  magnetic  field,  is  then 
given  by  equation  1.106. 

(1.106)  m  =  mQ  Xp 

where  xp  is  a  unit  vector  along  the  x  axis  in  the  rotating  frame,  and  mo  is  the  initial 
amplitude  of  the  isochromat  magnetization.  At  a  time  t  from  the  onset  of  the 
application  of  a  magnetic  field  gradient,  Gx,  in  the  x  direction,  the  magnetization  of 
the  isochromat  is 

(1.107)  m  =  mo  [cos(  7Gxxt)  ?tp  +  sin(  yGxxt)  yp  ] 

However,  we  choose  to  express  this  isochromatic  magnetization  in  a  more  compact 
form  which  is  given  by 

(1.108)  m  =  mo  exp(  i  yGxxt) 

Furthermore,  by  applying  a  magnetic  field  gradient,  Gy,  in  the  y  direction,  from  time 
t  =  0  to  t  =  Tp,  at  which  time  Gy  is  removed  and  Gx  is  applied,  the  magnetization  of 
the  isochromat  becomes,  at  time  t  ^  Tp 


(1.109) 


m  =  mo  exp{  i  [  yGxx(FTp)  +  yGyyTp]} 
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and  demonstrates  both  x  and  y  spatial  dependencies.  This  expression  is  simplified 
by  the  introduction  of  the  coordinates  of  k-space  which  are  given  by 


(1.110) 

kx  —  YGx(t-Tp) 

(1.111) 

ky  =  7^yTp 

so  that  by  substitution  equation  1 .109  becomes 

(1.112)  m  =  nio  exp[  i  (  kxx  +  kyy)] 

Neglecting  relaxation  effects,  the  magnetization  is  proportional  to  the  spin 
density,  p,  at  the  isochromat.  For  any  isochromat  in  the  xy  plane,  the  magnetization 
m(x,y)  is  therefore  given  by 

(1.113)  m(x,y)  p(x,y)exp[  i  (  kxx  +  kyy)] 

Moreover,  the  NMR  signal,  s(x,y),  from  an  isochromat  is  proportional  to  the 
isochromat  magnetization  m(x,y),  and  can  be  expressed  as 

(1.114)  s(x,y)  «  p(x,y)exp[  i  (  kxx  +  kyy)] 

Two-Dimensional  Fourier  Transform  Image  Reconstruction 

The  desired  NMR  image  is  a  plot  of  the  spin  density  distribution  in  the 
selected  image  plane,  and  so  p(x,y)  is  the  infomiation  that  must  be  extracted  from 
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the  NMR  signal.  Toward  this  end,  the  NMR  signal,  S,  from  the  entire  volume  of 
interest  can  be  obtained  by  summing  over  all  isochromats  as  shown  in  equation 

1.115 


(1.115) 


S(kx,  ky) 


oc 


Jp(x,y)exp[  i  (  kxx  +  kyy)]  dx  dy 

y 


and  can  be  easily  recognized  as  the  two-dimensional  (2D)  inverse  Fourier  transform 
(IFT)  of  p(x,y)  at  the  position  (kx,ky)  in  k-space.  In  order  to  reconstruct  an  image 
from  this  signal,  however,  it  is  necessary  to  obtain  the  entire  2D  IFT  of  p(x,y),  and 
so  the  signal  must  be  sampled  at  all  values  of  kx  and  ky.  This  can  be  accomplished 
by  applying  the  gradient  Gx  during  the  signal  acquisition  period  over  which  t  -  Tp 
spans  from  -Taq/2  to  Taq/2,  where  the  origin  has  arbitrarily  been  chosen  to  be  at  the 
center  of  the  acquisition.  As  a  result  of  this  gradient,  kx  is  spanned  over  the  range 

(1.116)  kx  =  -7GxTaq/2  to  +7GxTaq/2 
where  the  separation  of  sampled  points  is 

(1.117)  Akx  =  yGxAt 

In  this  expression,  At  is  the  digitization  interval  used  when  recording  the  NMR 
signal  and  is  also  known  as  DW,  the  dwell  time.  If  between  successive  data 
acquisition  steps  Gy  is  incremented  by  a  constant  amount  AGy,  so  that  over  the 
entire  image  acquisition  period  Gy  is  spanned  from  -Gymax  to  Gyraax>  then  ky  is 
spanned  over  the  range 
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(1.118)  ky  —  -'Y'GymaxTp  to  -'Y'GyniaxTp 

and  in  this  direction  of  k-space,  the  separation  of  sampled  points  is 

(1.119)  Aky  =  yAGyTp 

Thus,  Akx  and  Aky  have  constant  values,  and  the  sampling  of  k-space  is  rectangular 
as  shown  in  figure  1.16.  The  RF  pulse  sequence  used  to  sample  k-space  in  this 
way  is  shown  schematically  in  figure  1.17. 


Figure  1.16  Representation  of  k-space  with  rectangular  sampling 
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Figure  1.17  Schematic  representation  of  a  pulse  sequence  for  rectangular  sampling 
of  the  data  in  k-space  which  corresponds  to  a  transverse  imaging  plane. 

The  k-space  data  represented  in  figure  1 .16  is,  in  practice,  obtained  from  the 
NMR  signal,  S(kx,ky).  Equation  1.115  demonstrates  that  the  spatial  dependence  of 
the  spin  density  distribution,  p(x,y),  can  be  obtained  from  this  k-space  data  by 
means  of  a  Fourier  transform,  as  follows 


(1.120)  p(x,y)oc  [  [  S(kx,ky)exp[ -i  (kxx  +  kyy)]  dkx  dky 

kx  ky 

Although  the  effects  of  relaxation  have  been  neglected  in  this  discussion,  in 
practice,  relaxation  effects  provide  contrast  between  image  regions  depicting 
different  relaxation  environments,  and  can  greatly  increase  the  amount  of  diagnostic 
information  contained  in  an  image.  As  a  result,  relaxation  effects  can  be  of 
considerable  importance  to  NMR  imaging  applications  and  so  should  not  be 


overlooked. 
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Virtually  any  method  of  sampling  k-space  is  suitable  if  the  sampling 
efficiency  is  not  a  concern,  but  in  all  cases  there  are  constraints  which  affect  the 
quality  of  a  reconstructed  image.  Because  k-space  is  sampled  at  discrete  points,  the 
result  of  the  Fourier  transform  used  to  reconstruct  an  image  is  cyclic.  As  a  result, 
the  highest  frequency  signal  which  is  sampled  must  be  below  the  Nyquist 
frequency  to  avoid  aliasing  in  the  reconstructed  image  (60).  If  the  desired  field  of 
view,  FOV,  has  dimensions  X  by  Y  then  kx  and  ky  must  be  constrained  such  that: 


(1.121) 


kxmaxX 

2 


K 

Akx 


(1.122)  s  — 

z  Aky 


where  Akx  and  Aky  are  the  separations  of  sampled  points  in  k-space.  In  practice 
this  is  usually  achieved  by  using  magnetic  field  gradients  with  a  low  enough 
amplitude  that  the  FOV  is  larger  than  the  object  of  interest.  This  can  present 
another  problem,  however,  because  Akx  and  Aky  are  also  related  to  the  image 

resolution  that  will  be  achieved  as  indicated  in  figure  1.18.  Nonetheless,  a  property 
of  the  k-space  data  that  can,  in  general,  be  utilized  to  improve  the  image  quality,  is 
that  the  higher  frequency  data  points  are  usually  very  nearly  zero.  Points  in  k-space 
can  therefore  be  packed  with  zeroes  at  frequencies  higher  than  those  sampled  to 
effectively  increase  the  number  of  sample  points,  N.  As  a  result,  the  reconstructed 
image  has  closer  data  points  but  the  size  of  the  FOV  does  not  change.  This  does 
not,  however,  improve  the  actual  resolution  of  the  image  because  this  process 
increases  the  number  of  image  data  points,  or  pixels,  by  ideally  interpolating  points 
between  the  original  data.  As  a  result,  the  choice  of  the  magnetic  field  gradient 
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magnitudes,  the  sampling  rates,  and  the  number  of  points  to  sample  to  generate  an 
image,  must  be  a  balance  of  the  desired  resolution,  the  FOV,  and  the  time  to  acquire 
the  image  data. 


Figure  1.18  The  relationship  between  the  number  of  sample  points  and  the 
sampling  interval  in  k-space  (A)  and  the  resolution  and  the  size  of  the  FOV  of  the 
image  (B)  which  results  from  applying  a  Fourier  transform  to  the  data  in  k-space. 

Backprojection  Reconstruction 

An  alternative  pulse  sequence  for  sampling  k-space  is  as  shown  in  figure 
1.19  in  which  magnetic  field  gradients  are  applied  at  successively  incremented 
angular  displacements  around  the  image  plane.  With  such  a  pulse  sequence  k-space 
can  be  sampled  in  radial  lines  as  shown  in  figure  1.20.  An  image  can  be 
reconstructed  from  this  radially  sampled  k-space  data  by  interpolating  between 
sampled  data  points  and  applying  a  Fourier  transform.  However,  the  interpolation 
process  can  be  time  consuming,  and  generally  requires  a  number  of  assumptions  to 
be  made  about  the  form  of  the  data  in  order  to  predict  k-space  data  values  at  points 
which  lie  between  the  sampled  points  in  k-space.  An  image  is  more  readily 


reconstructed  from  radially  sampled  k-space  data  by  means  of  a  reconstruction  from 
projections  technique,  as  follows  (16). 


RF 
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Figure  1.19  Schematic  representation  of  a  pulse  sequence  for  polar  sampling  of  the 
data  in  k-space  which  corresponds  to  a  transverse  imaging  plane.  Values  of  Gx  and 
Gy  are  expressed  for  the  n^  of  N  data  acquisition  steps. 


Each  line  of  the  data  represented  in  figure  1.20  is  the  inverse  Fourier 
transform  of  the  projection  of  the  signal  from  the  object,  onto  a  line  defined  by  the 
gradient  direction.  The  signal  from  such  a  projection  at  an  angle  (f)  from  the  x  axis 
is  given  by: 


(1.123) 


S(t,cf>)  J  J  p(r,(}))  exp(  i  y  Gr  r  t)  dr  ds 


where  the  r,  s  and  (j)  are  spatial  coordinates  as  defined  in  figure  1.21,  and  p(r,(J))  is 
the  spin  density  in  the  image  plane,  and  is  therefore  the  information  that  is  displayed 
in  an  image. 


Figure  1.20  Representation  of  k-space  with  radial  or  polar  sampling 

Equation  1.123  can  be  simplified  by  substituting  a  function,  g(r,(j)),  for  the 
integral  in  s,  where  g(r,(j))  is  the  projection  of  p(r ,(}))  onto  a  line  at  an  angle  (|)  from 
the  x  axis,  and  is  given  by 

(1.124)  g(r,<]))  =  Jp(r,()))  ds 

S 


This  substitution  yields  the  expression 


(1.125) 


S(t,(j))  «  Jg(r,c}))  exp(  i  y  Gr  r  t)  dr 
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Figure  1.21  Definition  of  the  coordinates  r,  s  and  (J)  for  a  given  signal  projection,  in 
terms  of  the  fixed  cartesian  coordinates,  x  and  y. 


As  with  equation  1.109,  the  substitution  to  k-space  coordinates  can  be  made  such 
that 


(1.126)  kr=yGrt 


and  by  substitution  equation  1.125  becomes 


(1.127) 


S(kr,(|))  a  jg(r,(j))  exp(  i  kr  r)  dr 
r 


which  is  the  inverse  Fourier  transform  of  the  projection  signal,  g(r,({)).  The  desired 
image  data,  however,  is  proportional  to  p(x,y)  and  is  given  by  the  Fourier 
transform  of  S(kx,ky) 
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(1.128)  p(x,y)  a  [  [  S(kXrky)  exp[-i(kxx  +  kyy)]  dx  dy 

kx  ky 

and  so  a  change  of  coordinates  is  required  in  order  to  make  use  of  the  signal 
projection  profiles.  The  variables  in  the  cartesian  coordinate  system,  x  and  y,  can  be 
related  to  those  in  the  polar  system,  r  and  (|),  as  follows 


(1.129) 

x  =  r  cost)) 

(1.130) 

y  =  r  sin(j) 

and  the  element  of  space  over  which  the  signal  is  to  be  integrated  can  be  related  to 
polar  coordinates  as  shown  in  figure  1 .22. 


Figure  1.22  The  spatial  element  of  integration  as  expressed  in  cartesian  and  polar 
coordinates. 


Equation  1.128  can  thus  be  converted  into  polar  coordinates  to  give 
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(1.131) 


P(r,(f>)  a  J  J  S(kr,(J))  exp(-i  kr r  )  |r|  dr  d(j) 


The  inner  integral  can  be  recognized  as  the  Fourier  transform  of  S(kr,  ({))  with 
respect  to  r,  with  the  exception  of  the  term  r  that  has  been  introduced  by  the 

change  of  coordinates.  Relating  equation  1.131  to  1.127  one  can  define  a  quantity 
g*(r,(j))  with  the  expression 


(1.132) 


g*(r,CJ» 


S(kr,(j))  exp(-i  kr  r )  dr 


where  g*(r,(|))  is  the  filtered  projection  of  the  signal  from  the  object  of  interest.  The 
image  data  can  then  be  obtained  from  the  filtered  projections  by  integrating  over  all 
of  the  projection  angles 


n 

(1.133)  p(r,<|>)  oc  Jg*(r,((»  d((, 


In  practice,  however,  the  signal  projections  can  be  acquired  at  only  a  finite  number 
of  projection  angles  and  so  the  reconstructed  image  is  approximated  as 

N 

(1.134)  I(r,<t>)  = 

n=0 


where  A({)  is  the  angle  between  successive  projections. 
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Figure  1,23  An  image  plane  containing  three  cylindrical  objects  and  the 
corresponding  projection  profiles  at  three  polar  angles. 

In  the  work  described  in  Chapter  3,  the  number  of  projections  used,  N,  was 
8,  whereas  the  first  images  published  by  Lauterbur  (40)  were  reconstructed  from 
only  4  projections,  and  Bowtell  et  al.  (15)  applied  back-projection  reconstruction  to 
a  fast-imaging  technique  with  64  projections.  For  an  example  of  this  reconstruction 
technique  we  look  at  an  image  plane  containing  three  cylindrical  objects  as  shown  in 
figure  1.23.  When  three  projections  of  the  signal  from  these  objects  are  projected 
back  onto  the  image  plane  and  summed  as  in  equation  1 .134,  an  image  is  created  as 
shown  in  figure  1.24.  In  practice  however,  many  more  than  three  projections  are 
usually  required  to  generate  an  image.  As  can  be  seen  in  figure  1 .24,  the  difference 
in  signal  intensity  between  the  background  and  the  object  images  increases  with 
more  profiles  being  back- projected,  as  does  the  resolution  of  the  image. 
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Figure  1.24  The  reconstruction  of  the  original  image  from  three  projection  profiles, 
back-projected  onto  the  image  plane  and  summed. 

In  order  to  estimate  the  resolution  that  may  be  achieved  in  an  image  which 
has  been  reconstructed  by  means  of  a  back-projection  algorithm,  we  look  at  an 
image  plane  containing  a  single  circular  object  of  radius  d,  and  compute  the  image 
pixel  intensity  at  a  point,  P,  a  distance  D  from  the  object  center,  as  illustrated  in 
figure  1 .25.  We  can  choose  this  circular  object  to  have  a  triangular  cross  section  in 
every  plane  normal  to  the  xy  plane,  so  that  its  intensity  projection  profile  at  any 
angle  from  the  x  axis  can  be  approximated  by  a  triangular  function,  with  amplitude 
I,  where 

(1.135)  I  =  ^lr-dl  for-d^r^d 

(1.136)  1  =  0 


for  I  r  I  >  d 
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and  where  IQ  is  the  peak  profile  intensity,  and  r  is  the  distance  from  the  profile 
center. 


Figure  1.25  A  circular  object  within  an  image  plane,  and  definiuon  of  the  reference 
point  P. 

Figure  1.26  illustrates  a  case  in  which  a  single  profile  has  been  projected  onto  the 
image  plane  at  an  angle  (})  from  the  x  axis.  In  this  situation,  the  point  P  is  at  a 
distance  r  =  Dsintj)  from  the  center  line  of  the  projected  profile,  and  the  intensity,  Ip, 
is  given  by 


IP  =  ^  (d  -  Dsim))) 


(1.137) 


for  D  ^  d/sintj) 
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Figure  1.26  A  single  intensity  profile  of  the  object  illustrated  in  figure  1.25,  back 
projected  onto  the  image  plane  at  an  angle  (})  from  the  x  axis.  The  dashed  lines 
indicate  the  center  line  and  edges  of  the  back-projected  profile. 


Because  the  point  P  has  been  chosen  to  lie  on  the  x  axis,  after  two  profiles  at  angles 
0  and  kI2  radians  from  the  x  axis  have  been  back-projected  on  the  image  plane,  the 
intensity  Ip  is  equal  to  IQ.  The  intensity  at  the  object  center  in  this  case,  is  2Io.  The 
intensity  at  point  P  after  N  profiles  have  been  projected  onto  the  image  plane, 
however,  depends  not  only  on  the  number  of  back-projected  profiles,  but  also  on 
the  distance  D  from  the  object  center,  and  the  diameter,  d,  of  the  object. 


I w '  I 
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Number  of  Projections,  N  Intensity  at  point  P 

2  Io 

1°  +  <f(d-D) 

D  >  d 

D  <;  d 

4  Io 

Io  +  ~g— (d-Dsin(7r/4)) 

D  >  d/sin(7i/4) 

d  <  D  ^  d/sin(7t/4) 

Io  +  jf(d-D) 

+  -^(d-Dsin(7r/4)) 

D^d 

8  Io 

I0  +  ^(d-Dsin(7r/8)) 

D  >  d/sin(7r/8) 

d/sin(7t/4)  <  D  and 

1°  +  -g^(d-Dsin(7r/4)) 

+  ^(d-Dsin(jr/8)) 

D  £  d/sin(7r/8) 

d/sin(37i/8)  <  D  and 

D  £  d/sin(7r/4) 

I0  +  ^(d-Dsin(3n/8)) 

+  ^(d-Dsin(7t/4)) 

+  ^(d-Dsin(jt/8)) 

d  <  D  ^  d/sin(37r/8) 

Io  +  |j(d-D) 

+  ^“(d-Dsin(37t:/8)) 

+  ^(d-Dsin(K/4)) 

+  ^(d-Dsin(jt/8)) 

D  <;  d 

Table  1.1  Computed  intensities  at  a  point  P,  a  distance  D  from  an  object  of 
diameter  d,  after  N  signal  intensity  profiles  have  been  projected  onto  the  image 
plane. 
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This  theoretical  intensity,  Ip,  at  the  point  P,  determined  by  summing  the  intensity 
contributions  (computed  from  equation  1.137)  from  each  of  2,  4,  and  8  profiles 
which  are  back-projected  onto  the  image  plane  at  equal  angular  separations,  is  listed 
in  table  1.1.  Figure  1.27  illustrates  the  computed  intensity  distribution  in  an  image 
constructed  from  each  of  2,  4,  8  and  16  back-projection  profiles,  as  well  as  the 
theoretical  limit  that  is  approached  for  a  very  large  number  of  projections.  The 
intensity  at  the  object  center  after  N  profiles  have  been  back  projected,  is  equal  to 
NI0. 


Relative  Distance  from  Object  Center  (D/d) 

Figure  1.27  Theoretical  intensity  profiles  for  an  image  of  the  circular  object 
illustrated  in  figure  1.25,  reconstructed  by  means  of  a  back-projection  algorithm. 
The  intensities  were  computed  with  the  equations  listed  in  table  1.1.  The  distance 
from  the  object  center  is  expressed  relative  to  the  object  diameter  and  the  intensities 
are  expressed  relative  to  the  intensity  at  the  object  center.  The  number  of  back- 
projections  represented  in  each  of  these  curves  is  indicated  on  the  right. 
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Figure  1 .27  demonstrates  the  improved  image  contrast,  and  improved  definition  of 
the  object  edge,  that  can  be  obtained  by  increasing  the  number  of  back -projected 
profiles.  Moreover,  because  the  image  intensity  near  the  edge  of  this  object 
depends  both  on  the  size  of  the  object  and  on  the  number  of  back-projected  profiles, 
the  ability  to  resolve  two  such  objects  that  are  in  close  proximity  will  depend  on 
these  factors  as  well.  The  minimum  distance  between  two  objects  which  can  be 
resolved,  however,  is  determined  by  the  separation  of  data  points  in  the  projections 
themselves,  which  for  the  deuteron  images  discussed  in  Chapter  3  was 
approximately  1 .6  mm.  However,  for  cases  in  which  the  image  resolution  is  greater 
than  the  separation  of  the  profile  data  points,  the  dimensions  of  the  objects  being 
resolved  in  an  image  should  be  specified  for  the  quoted  image  resolution  to  be 
meaningful. 

For  the  purposes  of  obtaining  deuteron  images  of  instilled  alveolar  HOD,  as 
discussed  in  Chapter  3,  a  back-projection  reconstruction  from  8  projections 
provided  sufficient  image  quality,  because  of  the  unique  nature  of  this  NMR 
application.  Firstly,  when  imaging  the  2H  distribution  in  a  body  there  is  essentially 
no  detectable  NMR  signal,  with  the  exception  of  2H  which  has  been  artificially 
placed  into  the  body.  Secondly,  the  HOD-based  solution  we  instilled  into  a  lung 
lobe  was  well  localized,  and  therefore  generally  provided  only  one  object  of  interest 
in  our  field  of  view.  In  this  particular  situation,  the  background  pixel  intensity  in 
the  reconstructed  image  was,  at  most,  only  1/8  that  of  the  pixel  intensity  of  the 
object  of  interest,  as  demonstrated  in  figure  1.27.  Moreover,  the  resolution 
achieved  in  an  image  reconstructed  from  8  projections  was  estimated  from  images 
of  1.3  cm  diameter  phantoms  to  be  roughly  3  mm,  suitable  for  the  purposes  of  the 
study  discussed  in  Chapter  3. 


78 


1.4  Application  of  NMR  to  Studies  of  Lung  Water 

The  NMR  applications  described  in  Chapters  2  and  3  include  proton  spin- 
echo  imaging,  deuteron  FID  imaging  and  deuteron  total  signal  intensity 
measurements,  for  monitoring  the  clearance  of  fluid  which  has  been  artificially 
introduced  into  one  lung  lobe  of  an  anesthetized  cat.  Chapter  2  discusses  a  study  in 
which  the  instilled  alveolar  fluid  is  composed  of  blood  serum  taken  from  the  cat 
itself,  and  Chapter  3  discusses  a  study  in  which  a  serum-like  solution  containing 
50%  D2O  is  instilled  into  the  airspaces.  In  order  to  monitor  the  net  clearance  of 
alveolar  fluid,  and  to  monitor  bi-directional  fluid  fluxes  across  the  pulmonary  air- 
blood  barrier,  a  combination  of  proton  and  deuteron  NMR  techniques  were 
employed. 

Proton  spin-echo  imaging  was  employed  both  for  first-echo  image  pixel 
intensity  measurements,  and  for  proton  transverse  relaxation  time  measurements,  to 
enable  monitoring  of  water  proton  environments  in  the  lung.  The  spin-echo 
sequence  employed  was  a  CPMG  sequence  in  which  the  selective  90°  excitation 
pulse  amplitude  was  modulated  by  a  Gaussian  function  in  order  to  define  a  1  cm 
thick  transverse  imaging  plane.  Magnetic  field  gradients  were  employed  to  enable 
rectangular  sampling  of  k-space.  Moreover,  the  signal  acquisition  was  gated  to  the 
cat's  respiratory  motion  so  that  the  NMR  signal  was  only  sampled  when  the  cat's 
lungs  were  at  functional  residual  capacity  (FRC).  However,  the  net  transverse 
magnetization  immediately  after  the  90°  RF  pulse  depends  on  the  time  that  has  been 
allowed  for  the  spin  system  magnetization  to  return  to  thermal  equilibrium  after  the 
previous  90°  RF  pulse.  As  a  result,  it  was  necessary  to  employ  "spin-conditioning" 
in  which  successive  pulse  sequences  were  applied  with  a  constant  repetition  time, 
Tr,  regardless  of  whether  the  cat's  respiratory  cycle  was  at  FRC  and  the  NMR 
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signal  was  to  be  sampled.  In  this  way,  the  transverse  magnetization,  M,  after  each 
90°  RF  pulse  did  not  vary  with  the  cat's  respiratory  rate  but  had  a  constant  value 
given  by 

(1-138)  M  =  M0[  1  -exp(-TR/Ti)] 

where  Mq  is  the  magnitude  of  the  spin  system  magnetization  at  thermal  equilibrium. 

The  deuteron  imaging  technique  employed  was  designed  to  minimize  the 
time  between  the  90°  excitation  pulse  and  sampling  of  the  NMR  signal,  thereby 
minimizing  the  signal  loss  due  to  the  rapid  deuteron  transverse  magnetization  decay. 
Deuteron  images  were  obtained  by  applying  non-selective  90°  RF  pulses  and 
acquiring  FIDs  in  the  presence  of  magnetic  field  gradients  which  enabled  radial 
sampling  of  k-space.  As  a  result,  deuteron  images  were  reconstructed  by  means  of 
a  back-projection  algorithm.  Deuteron  signal  intensity  measurements  were  obtained 
in  a  similar  manner,  with  the  exception  that  spatial  encoding  by  means  of  magnetic 
field  gradients  was  not  required.  The  deuteron  images  and  signal  intensity 
measurement  thus  acquired,  were  employed  for  monitoring  the  time  course  of 
changes  in  the  deuteron  content  of  the  lung,  after  instillation  of  a  serum- like  solution 
containing  D2O. 

As  will  be  discussed  in  the  following  chapters,  these  deuteron  and  proton 
NMR  techniques  have  enabled  us  to  obtain  time  course  observations  of  alveolar 
fluid  clearance,  and  of  fluid  absorption  and  secretion  across  the  pulmonary  air- 
blood  barrier.  Moreover,  these  observations  were  obtained  non-invasively,  in  a 
living  animal,  for  the  first  time. 
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2.0  Evaluation  of  the  Effects  of  PAF  on  Alveolar  Fluid  Clearance 
Using  NMR  Imaging.1 


2.1  Introduction 

Although  proton  nuclear  magnetic  resonance  (NMR)  imaging  is  particularly 
suited  to  studies  of  the  soft  tissues  in  the  body,  image  quality  is  greatly  limited  in  in- 
vivo  lung  studies  by  the  paucity  of  water  protons  and  by  susceptibility  problems,  as 
well  as  by  artifacts  caused  by  respiratory  motion,  cardiac  motion,  and  pulsatile 
blood  flow.  However,  the  NMR  signal  intensity,  depending  as  it  does  on  three 
parameters,  namely,  the  local  proton  density,  p,  the  longitudinal  relaxation  time,  Ti, 
and  the  transverse  relaxation  component,  T2  (1)  can  also  provide  important 
information  through  the  relaxation  times.  These  times  have  been  shown  to  be 
sensitive  to  the  water  content  of  tissue  (16,  29),  the  mobile  protein  content  of  fluid 
(10, 20)  and  air-fluid  interfaces  (15). 

Incorporating  techniques  to  reduce  artifacts  and  increase  sensitivity,  NMR 
signal  intensity  measurements,  both  in  excised  lungs  (7)  and  in-vivo  (8,  25),  have 
been  shown  to  be  well  correlated  with  gravimetric  lung  water  measurements. 
Relaxation  times  of  lung  water  protons  have  also  been  shown  to  be  well  correlated 
with  gravimetric  measurements  (29)  and  have  been  used  to  differentiate  between 
hydrostatic  and  permeability  edema  (25).  Moreover,  in  normal  rat  lungs  in-vivo, 
two  transverse  relaxation  components  have  been  demonstrated  (27).  In  only  a  few 


1  A  version  of  this  chapter  has  been  accepted  tor  publication.  P.  W.  Stroman,  P.  S.  Allen,  D.  C. 
Lien,  G.  Machin,  S.  F.  P.  Man.  Evaluation  of  the  Effects  of  PAF  on  Alveolar  Fluid  Clearance 
Using  NMR  Imaging.  Journal  of  Applied  Physiology. 
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cases,  though,  were  both  quality  images  and  relaxation  time  measurements  obtained 
(23,  25)  or  were  NMR  techniques  applied  to  monitor  serial  changes  in  lung  water 
(8,  23,  27).  One  purpose  of  this  study  was  to  improve  on  the  quality  of  both  the 
images  and  the  relaxation  time  measurements  of  alveolar  fluid  protons,  so  that 
modifications  to  the  time  course  of  alveolar  fluid  clearance  by  the  administration  of 
platelet-activating  factor  could  be  studied  in-vivo. 

Platelet-activating  factor  (PAF)  is  a  family  of  ether-lipids  whose  active 
component  is  1-O-Hexadecyl/octadecy  1-2-ace tyl-sn-glycero-3-phosphocholine  (24). 
PAFs  are  produced  by  a  wide  variety  of  cells  in  the  body  (21,  24)  including 
alveolar  type  II  epithelial  cells  (13,  14).  It  has  been  shown  that  intra-airway  PAF 
can  cause  bronchoconstriction  (3)  and  inflammation  in  the  airways  and  capillary 
spaces  (30).  Because  the  role  of  PAF  in  the  development  and/or  resolution  of 
pulmonary  edema  has  not  been  previously  evaluated,  we  studied  the  effects  of  PAF 
at  a  dosage  that  is  comparable  to  that  used  in  studies  of  the  effects  of  PAF  on  lung 
endothelial  and  epithelial  permeability  (5,9). 

Our  findings,  using  all  the  NMR  parameters  at  our  disposal,  were  that  the 
addition  of  PAF  to  instilled  alveolar  fluid  dramatically  increased  water  clearance 
from  the  lung  but  had  little  or  no  effect  on  the  clearance  of  albumin  from  the  lung. 
From  the  gravimetric  and  NMR  data  acquired  at  the  end  of  the  experimental  time 
course,  we  were  able  to  demonstrate  a  correlation  between  NMR  intensity 
measurements  and  gravimetric  measurements  of  the  instilled  alveolar  fluid,  both 
when  serum  was  instilled  (p  >  0.40)  and  when  serum  plus  PAF  was  instilled  (p  > 
0.15),  where  the  p  values  quoted  are  the  results  of  unpaired,  2  tailed,  t  tests.  This 
correlation  plot  had  a  slope  of  0.95  and  a  regression  coefficient,  r,  of  0.74  (n  =  6). 
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2.2  Methods 

2.2.1  Animal  Model  Of  Alveolar  Fluid  Clearance 

We  have  studied  alveolar  fluid  clearance  in  anesthetized  mongrel  cats  (3.0- 
6.0  kg  body  wt.)  using  a  model  similar  to  that  described  by  Matthay  et  al.  for  the 
sheep  (Matthay  1985).  The  cat  was  anesthetized  with  an  intramuscular  injection  of 
ketamine  (Ketalean  25  mg/kg),  xylazine  hydrochloride  (Rompun  0.5  mg/kg)  and 
acepromazine  (Atravet  1  drop/4  kg).  Anesthesia  was  maintained  with  the  same 
mixture  as  required  (~  25%  initial  dose/hour).  The  supine  cat  was  intubated,  but 
breathed  spontaneously.  Polyethylene  catheters  (PE-90,  o.d.  1.27  mm)  were 
placed  in  an  external  jugular  vein  and  a  femoral  artery  for  the  administration  of 
fluids  and  hemodynamic  monitoring. 

A  7-10  ml  sample  of  venous  blood  was  taken  from  the  animal  in  order  to 
prepare  4  -  5  ml  of  serum.  The  instillate  was  prepared  from  this  serum  by  adding  1 
mg  of  Evan's  Blue  dye  and  3  |iCi  of  125I-albumin  (Merck  Frosst,  Kirkland  Quebec) 
to  label  the  instilled  albumin.  For  experiments  which  included  PAF,  250  |il  of  10'5 
M  PAF  was  also  added  to  the  serum  just  prior  to  instillation;  a  dose  range  of  0.17 
to  0.31  pg  PAF/kg  body  weight  was  used  (mean  =  0.23  jig  PAF/kg).  The  volume 
of  solution  instilled  into  the  lung  was  0.7  ml/kg  body  weight  and  the  remaining 
aliquot  (1  ml)  was  used  for  125 1  and  water  content  determinations. 

Prior  to  fluid  instillation,  the  cat  was  ventilated  manually  6  -  7  times  with  a 
tidal  volume  of  10  ml/kg  body  weight.  A  PE-90  catheter  was  inserted  into  a  lower 
lung  lobe  (right  lung,  27  times  out  of  36  experiments)  and  the  fluid  was  injected 
with  a  syringe  over  15  sec,  followed  by  the  insufflation  of  several  ml  of  room  air. 
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The  precise  amount  of  fluid  instilled  was  determined  by  weight.  NMR  imaging 
and/or  gravimetric  analysis  of  the  lung  were  then  carried  out,  as  described  below, 
during  the  succeeding  four  hour  period. 

Five  minutes  prior  to  the  end  of  the  experiment  the  cat  was  given  5000  units 
of  heparin  intravenously.  Just  before  the  cat  was  sacrificed,  venous  blood  samples 
were  taken  for  determinations  of  hematocrit,  blood  density,  hemoglobin  content, 
and  125I  and  fractional  water  contents.  The  chest  was  opened  with  a  mid-line 
sternotomy  and  the  lungs  were  tied  off  at  the  hila  and  removed.  The  distribution  of 
the  instilled  fluid  in  the  experimental  lobe  was  identified  by  the  Evan's  Blue  dye. 

The  lung  lobes  were  processed  to  determine  the  clearance  rates  of  liquid  and 
125I-albumin  (17).  The  lobes  were  weighed,  added  to  a  known  amount  (~100  g)  of 
distilled  water,  homogenized,  and  half  of  the  homogenate  was  centrifuged  to  obtain 
a  clear  red  supernatant.  The  fractional  water  content  of  the  homogenate,  and 
125Iodide,  fractional  water  content,  density,  and  hemoglobin  content  of  the 
supernatant,  were  determined  for  each  lobe.  We  computed  the  total  weight  of  blood 
in  each  lobe  using  the  calculations  described  by  Pearce  et  al.  (22)  but  we  did  not 
include  a  correction  factor  for  pulmonary  hematocrit.  Within  these  calculations,  we 
computed  the  hemoglobin  concentration  in  the  homogenate  from  the  measured 
hemoglobin  concentration  in  the  supernatant  as  was  done  by  Selinger  et  al  (26). 

To  avoid  interference  from  Evan's  Blue  dye  the  cyanmethemoglobin 
concentration  in  the  supernatant  was  determined  from  the  light  absorbance  at  420 
nm  instead  of  the  usual  540  nm.  Using  this  absorbance,  Evan's  Blue  dye  at  a 
concentration  of  1  mg/4  ml  did  not  interfere  with  the  measurement  of  hemoglobin, 
and  was  accurate  to  within  3%  for  hemoglobin  concentrations  of  less  than  0.1 
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mg/dL.  The  concentration  of  dye  in  the  supernatant  in  these  studies  was  typically  in 
the  order  of  0.002  mg/ml. 

The  amounts  of  125Iodide  in  the  blood  samples,  instillate,  homogenate 
samples  and  supernatant  samples  were  determined  with  gamma  counts  of  samples 
with  uniform  geometry  using  a  LKB  Wallac,  model  1282  Compugamma.  The 
counts  were  corrected  for  background,  counter  dead  time  and  isotope  decay. 

Fractional  water  contents  of  samples  were  determined  by  weighing  the 
samples  before  and  after  drying  to  a  constant  weight. 

The  number  of  cats  included  in  each  experimental  group  is  listed  in  Table 
2.1.  Whereas  data  from  NMR  measurements  were  consecutive  45  minute  averages 
that  serially  spanned  a  4  hour  period,  the  gravimetric  measurements  provided  only 
discrete  data  points  and  so  were  obtained  at  either  1  hour  or  4  hours  after  fluid 
instillation  to  allow  comparison  to  the  NMR  data.  Cats  which  received  serum 
formed  the  control  study  group,  cats  which  received  serum  plus  PAF  formed  the 
first  experimental  study  group,  and  cats  which  received  serum  plus  lyso-PAF,  a 
biologically  inactive  form  of  PAF,  formed  the  second  experimental  study  group. 
Each  of  these  study  groups  is  further  subdivided  in  Table  2.1  into  a  1  hour  sacrifice 
and  a  4  hour  sacrifice  subset,  although  in  the  1  hour  subset  only  4  cats  were 
studied,  2  were  given  serum  and  2  were  given  serum  plus  PAF.  Moreover,  this  1 
hour  subset  was  only  evaluated  gravimetrically,  leaving  the  early  NMR  data  points 
to  be  obtained  from  the  4  hour  subset. 

Table  2.1  shows  how  the  aforementioned  study  groups  were  further 
subdivided  into  those  evaluated  gravimetrically,  those  evaluated  by  NMR,  and  those 
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evaluated  by  both  methods.  In  addition,  8  cats  were  assessed  for  histology  (2  for 
saline,  3  each  for  autologous  serum  plus  Evan's  Blue  dye  and  autologous  serum 
plus  Evan's  Blue  dye  plus  PAF)  so  that  morphological  changes  in  the  lung  caused 
by  the  experimental  procedure,  or  by  PAF,  might  be  detected.  The  excised  lungs 
were  fixed  in  2.5%  glutaraldehyde  and  1.5%  formaldehyde  in  a  sodium  cacodylate 
buffer,  osmolarity  295  mOsmol/kg  H2O  and  pH  of  7.4. 


Serum  + 
Lyso-PAF 

Data 

Serum  Instilled 

Serum+  PAF  Instilled 

Instilled 

Saline  Instilled 

Obtained 

1  hour  4  hours 

1  hour  4  hours 

4  hours 

1  hour  4  hours 

NMR  only 

0  2* 

0  1 

0 

0  0 

Gravimetric  _  _ 

only 

2  6 

2  5 

2 

0  0 

NMR  and  _ 

Gravimetric 

0  3 

0  4 

1 

0  0 

Histology 

1  2 

1  2 

0 

1  1 

*One  cat  had  2  lung  regions  containing  fluid  monitored  with  multi-slice  imaging, 
so  3  fluid  clearance  curves  were  obtained  for  this  group 


Table  2,1  The  number  of  cats  included  in  each  of  the  experimental  groups 


2.2.2  NMR  Imaging 

Imaging  was  carried  out  on  a  modified  Bruker  CXP  spectrometer  (Bruker 
Instruments  Inc.  Billerica,  MA)  with  a  40  cm  bore  magnet  at  a  static  field,  Bo,  of 
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2.35  T.  The  anesthetized  supine  cat  was  supported  by  a  curved  plexiglass  tray 
inside  the  magnet  bore.  Hydrogen  nuclei  in  a  1  cm  thick  transverse  imaging  slice 
were  excited  uniformly  with  selective  radio  frequency  (RF)  pulses  transmitted  by  a 
large  circumscribing  coil.  To  provide  high  receiver  sensitivity  in  the  dorsal  half  of 
the  thorax,  a  single  turn  rectangular  surface  coil  (4  cm  x  15  cm),  curved  around  the 
thorax,  served  as  the  receiver.  To  avoid  inductive  coupling  between  the  two  coils, 
active  switching  to  the  proton  resonant  frequency,  under  the  control  of  the  pulse 
program,  was  employed  between  transmission  and  data  acquisition.  Respiratory 
gating  was  brought  about  by  means  of  a  flexible  rubber  bellows  placed  around  the 
cat's  abdomen  and  linked  to  a  pressure  transducer  whose  analog  output  was 
converted  to  TTL  for  triggering  the  pulse  programmer.  The  pulse  sequence  was 
applied  at  regular  intervals,  regardless  of  the  cat's  breathing,  but  the  NMR  signal 
was  only  acquired  when  the  cat's  thorax  was  at  functional  residual  capacity  (FRC). 

Our  imaging  pulse  sequence  (Figure  2.1)  was  based  on  a  Carr-Purcell- 
Meiboom-Gill  (CPMG)  train  of  pulses  (19)  which  had  been  modified  for 
respiratory  gating  and  multi-slice  imaging.  The  selective  90°  RF  pulse  was 
gaussian  shaped  with  a  7.68  msec  duration  whereas  the  non-selective  180°RF 
pulses  were  square  pulses,  each  with  a  duration  of  350  |is.  The  inter-echo  interval, 
Te,  was  26.3  msec,  and  each  of  the  spin-echo  signals  was  acquired  to  form  a  series 
of  spin-echo  images.  Four  signal  acquisitions  with  a  repetition  time,Tp,  of  1.0  sec 
were  averaged  to  improve  the  signal  to  noise  ratio.  The  data  set  for  each  spin-echo 
image,  acquired  as  a  64  x  256  matrix,  was  zero-filled  to  create  a  128  x  256  matrix 
which  was  Fourier  transformed  to  create  a  128  x  256  pixel  image  (10  cm  x  20  cm 
field  of  view). 
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Sets  of  8  spin-echo  images  were  obtained  for  each  of  4  contiguous  slices. 
The  slices  were  interleaved  in  non-contiguous  pairs  thus  allowing  us  to  image  the 
entire  caudal  half  of  the  thorax  every  30  to  40  minutes. 


Multi-Slice  Loop  (x  2) 


Test  Respiratory  Gating  Signal 


Figure  2,1  Schematic  of  the  pulse  sequence  used  to  obtain  8  spin-echo  images  of 
each  of  2  non-contiguous  slices.  The  TTL  respiratory  gating  signal  indicates  when 
the  cat's  lungs  are  not  at  FRC.  In  this  case,  spin  conditioning  takes  place  and  no 
gradients  are  applied  and  the  NMR  signal  is  not  acquired. 


The  lungs  were  imaged  before  the  instillation  of  serum  solution  into  one 
lung  lobe,  and  then  continuously  over  the  subsequent  4  hour  period.  The  First  spin- 
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echo  image  of  each  set  was  used  for  monitoring  the  lung  water  content. 
Rectangular  regions  of  interest  were  defined  using  a  grid  as  shown  in  Figure  2.2, 
with  each  region  typically  containing  300  to  400  pixels  (1 .5  to  2.0  cc).  Most  of  the 
fluid-filled  lung  lobe  fell  into  one  of  these  regions  and  the  equivalent  part  of  the 
contralateral  control  lobe  was  also  arranged  to  be  a  region  of  interest.  Several 
regions  of  interest  were  analyzed,  namely,  the  experimental  lung  lobe,  the 
contralateral  control  lobe  and  several  adjacent  regions  of  the  chest  wall  and  back. 
Each  region  of  interest  was  then  characterized  by  the  sum  of  the  pixel  intensities  in 
that  region.  The  background  intensity,  Ibg,  was  calculated  from  the  set  of  control 
images  taken  prior  to  instillation,  and  subtracted  from  the  total  intensity  of  the  fluid 
filled  lung  lobe  region  to  determine  IinSt,  the  intensity  of  the  NMR  signal  produced 
by  the  instilled  fluid  in  the  lung  lobe.  The  calculated  IinSt  values  were  then  plotted 
as  a  function  of  the  time,  At,  between  instillation  of  the  fluid  and  the  mid-point  of 
the  image  acquisition  period.  Because  the  fluid  clearance  over  4  hours  was  non¬ 
linear  in  general,  the  value  of  IinSt(At=0)  was  approximated  by  linear  extrapolation 
from  the  two  values  of  IinSt(At)  closest  to  At=0.  The  estimated  value  of  Iinst(At=0) 
was  used  to  normalize  the  time  course  of  the  fluid  clearance  from  the  lung  lobe. 

The  same  regions  of  interest  were  also  used  for  transverse  relaxation  time 
measurements.  Using  the  total  pixel  intensity  of  the  chosen  region,  the  eight  spin- 
echo  images  provided  a  decay  curve  for  the  region-of-interest  magnetization,  from 
which  average  transverse  relaxation  time  components  were  derived  for  that  region 
using  a  non-linear  least  squares  algorithm  (4). 
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Before  Instillation 


Time  At  After  Instillation 


Intensity  =  I  Intensity  =  Ipp(At)  Intensity  =  I^q 

Figure  2.2  The  method  used  for  defining  regions  of  interest  in  transverse  images  of 
a  cat's  thorax.  Only  the  first  spin-echo  images  are  used  for  monitoring  the  time 
course  of  changes  in  the  lung  water  content.  The  intensity  of  the  acquired  signal 
which  is  due  to  the  instilled  fluid  is  given  by:  Iinst(At)  =  IpF(At)  -  Ibg- 


2.3  Results 

Examples  of  NMR  images  obtained  before  and  after  instilling  fluid  into  a 
cat's  lower  lung  lobe  are  shown  in  Figures  2.3a  to  2.3d.  The  time  dependence  of 
the  integrated  instillate  NMR  signal  intensity  data,  linst,  grouped  into  45  minute 
interval  averages,  is  shown  in  Figure  2.4.  For  comparison,  results  of  gravimetric 
measurements  obtained  at  1  and  4  hours  after  fluid  instillation,  are  also  shown.  All 
of  the  results  are  expressed  as  the  mean  of  n  experiments,  plus  or  minus  the 
standard  error  of  the  mean.  For  the  serum  instillate,  the  integrated  NMR  intensity 
data  showed  that  86%  ±  6%  (n=6)  remained  in  the  experimental  lung  lobe  between 
3.25  and  4  hours  after  instillation,  which  compared  favorably  to  gravimetric 
measurements  of  78%  ±  6%  (n=9)  at  4  hours  after  instillation.  With  PAF  added  to 
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Figure  2.3  Transverse  images  of  a  cat’s  thorax;  a)  before  instillation 
of  fluid  (upper  left),  b)  20  minutes  after  instillation  of  fluid 
(upper  right),  c)  120  minutes  after  instillation  of  fluid  (lower 
left),  and  d)  200  minutes  after  instillation  of  fluid  (lower  right). 
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the  instillate,  the  integrated  NMR  intensity  data  revealed  that  only  35%  ±  4%  (n=5) 
remained  in  the  lung  lobe  between  3.25  and  4  hours  after  instillation,  compared  to  a 
measurement  of  47%  ±  6%  (n=9)  after  4  hours  from  the  gravimetric  method.  That 
the  addition  of  PAF  to  the  instillate  dramatically  increased  fluid  clearance  over  a  4 
hour  period  was  demonstrated  by  both  measurement  techniques.  The  addition  of 
PAT  to  the  serum  instillate  introduced  a  bi-phasic  character  to  the  water  clearance 
which  could  be  readily  detected  from  the  sequential  NMR  images.  The  fast 
clearance  phase  had  a  ti/2  of  approximately  30  min,  whereas  the  slow  clearance  rate 
was  similar  to  the  clearance  rate  for  serum  alone  with  a  ti/2  of  approximately  670 
min. 


Figure  2,4  The  time  course  of  alveolar  fluid  clearance  as  determined  by  in-vivo 
NMR  imaging  and  by  gravimetric  measurements  of  excised  lungs.  Circles  indicate 
averaged  NMR  intensity  data  (serum  n  =  6,  serum  plus  PAF  n  =  5).  Triangles 
indicate  averaged  gravimetric  measurements  (each  group:  1  hour  n  =  2,  4  hours  n  = 
9).  The  dotted  line  shows  the  fitted  clearance  curve  without  PAF:  %  Remaining  = 
100exp(-t/966).  The  dashed  line  shows  the  fitted  clearance  curve  with  PAF:  % 
Remaining  =  55exp(-t/40)  +  45exp(-t/966). 
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When  lyso-PAF  was  added  to  the  instillate  it  did  not  appear  to  affect  fluid 
clearance;  gravimetrically,  71%  ±  19%  (n=3)  of  the  instilled  water  remained  after  4 
hours.  NMR  registered  a  corresponding  value  of  approximately  90%  (n= 1) . 

The  amount  of  125I-albumin  (in  the  gravimetric  groups)  remaining  in  the 
experimental  lung  lobe  after  4  hours  suggests  that  the  introduction  of  PAF 
decreased  the  residual  125I-albumin  from  91%  ±  3%  (n=9)  with  pure  serum 
instillate  to  82%  ±  5%  (n=9)  when  PAF  was  added.  Nonetheless,  the  significance 
of  this  discrimination  is  low  (p  =  0.06).  The  amount  of  radioactivity  we  were  able 
to  recover  was  ^  90%  (range  90%  to  100%)  of  that  instilled,  and  the  amount 
present  in  the  thyroid  (  as  an  index  of  free  radioiodine)  was  negligible  (<0.01%). 
We  were  therefore  confident  that  leaching  of  the  isotope  from  radiolabelled  albumin 
was  minimal. 

The  decay  of  the  proton  transverse  magnetization  from  the  region  of  interest 
in  the  fluid-instilled  lung  could  be  reproducibly  resolved  into  two  exponential 
components,  with  characteristic  times  T2F  (fast  relaxing)  and  T2S  (slow  relaxing). 
Measurements  of  T2F  and  T2S,  again  time  averaged  over  the  intervals  used  for  the 
integrated  signal  intensity  averaging  (45  minutes),  are  shown  as  a  function  of  time 
following  instillation  in  Figure  2.5.  It  is  clear  from  this  figure  that  any  differences 
in  T2F  due  to  PAF  addition  to  the  instillate  are  too  small  to  be  resolved  throughout 
the  4  hours  of  measurement.  Resolution  of  the  PAF  induced  differences  in  T2S 
appears  to  be  marginal,  however  the  uncertainties  were  too  large  to  claim  a 
significant  difference.  Throughout  the  4  hour  observation  period  the  more  slowly 
relaxing  T2S  component  accounted  for  approximately  40%  of  the  lung  proton  signal 
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in  the  region  of  interest,  irrespective  of  whether  serum  or  serum  plus  PAP  was 
instilled. 

Arterial  blood  gases,  pH  and  blood  pressure  were  monitored  in  experiments 
in  which  only  gravimetric  measurements  were  obtained.  All  of  these  monitored 
parameters  were  found  to  be  stable  and  without  systematic  change  throughout  the 
course  of  the  experiments. 


Figure  2.5  The  time  dependence  of  the  time  averaged  transverse  relaxation  times  of 
cat  lungs  containing  instilled  fluid.  The  filled  circles  represent  serum  instilled  into 
the  lungs.  The  open  circles  represent  serum  plus  PAF  instilled  into  the  lungs. 
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Figure  2.6a  Photomicrograph  of  cat  lung,  4  hours  after  instillation  of  autologous 
serum  plus  PAF. 
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Figure  2.6b  Photomicrograph  of  cat  lung,  4  hours  after  instillation  of  autologous 


serum. 
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Histological  sections  of  lungs,  examined  in  a  single  blind  fashion,  made  at  1 
and  4  hours  after  instillation  of  saline  plus  dye  showed  no  gross  morphological 
evidence  of  inflammation,  and  verified  that  instillate  was  delivered  to  the  alveoli. 
An  accumulation  of  neutrophils  in  the  vascular  spaces  was  observed  irrespective  of 
whether  the  instillate  was  either  serum  or  serum  plus  PAF,  but  no  accumulation  of 
neutrophils  in  the  airspaces  or  interstitium  could  be  detected.  Photomicrographs  of 
the  PAF  treated  lungs  in  Figure  2.6a  do  not  demonstrate  any  excess  of  fluid  in  the 
interstitium  or  vascular  spaces,  while  the  non-PAF  treated  lungs  in  Figure  2.6b 
show  distended  alveolar  capillaries  and  venules. 


2,4  Discussion 

The  NMR  images  we  have  obtained  have  enabled  us  to  observe,  non- 
invasively,  the  time  course  of  lung  alveolar  fluid  clearance  in  living  cats.  After  the 
instillation  of  either  serum,  or  serum  plus  lyso-PAF  into  the  airspaces,  we  observed 
a  mono-phasic  clearance  with  a  half  time  (ti/2)  of  approximately  11  hours.  In 
comparison,  as  summarized  by  Matthay  et  al.  (18),  alveolar  liquid  clearance  in  dogs, 
sheep  and  rabbits  proceeds  with  an  estimated  ti/2  of  18  hours,  9  hours  and  6  hours 
respectively.  With  the  addition  of  PAF  to  the  instilled  fluid,  however,  we  observed 
a  dramatic  change  to  the  nature  of  the  fluid  clearance  from  the  alveolar  space  of  the 
cat.  With  PAF  the  clearance  was  bi-phasic,  having  an  initial  rapid  clearance  with  a 
1 1/2  of  approximately  30  minutes  followed  by  a  clearance  which  was  similar  to  that 
observed  when  either  serum  or  serum  plus  lyso-PAF  was  instilled. 

The  clearance  of  protein  from  the  alveolar  space  was  seen  to  proceed  at  a 
much  slower  rate  than  that  of  the  alveolar  fluid  itself,  independent  of  the  species 
studied.  Protein  clearance  has  been  reported  to  be  1  %/hour  in  sheep  (17)  and 
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2%/hour  in  the  anesthetized  dog  (30).  In  cats  we  observed  an  alveolar  125I-albumin 
clearance  rate  of  2.3%/hour  when  serum  alone  was  instilled,  and  this  clearance  did 
not  appear  to  have  been  significantly  altered  by  PAF. 

The  influence  of  PAF  on  the  fluid  and  protein  clearances,  as  well  as  that  of 
lyso-PAF,  which  has  similar  biochemical  but  not  biological  properties,  can  be  better 
but  not  totally  understood  from  our  data.  With  PAF,  there  was  no  significant 
increase  in  albumin  clearance,  leading  to  the  postulate  that  the  increased  fluid 
clearance  occurred  not  because  of  a  change  in  the  sieve  properties  of  the  epithelium. 
Lyso-PAF  on  the  other  hand  did  not  show  an  effect  on  fluid  clearance  either.  We 
can  therefore  postulate  that  the  increase  in  fluid  clearance  caused  by  PAF  is  by  a 
receptor  mediated  mechanism.  Nonetheless,  we  cannot  be  certain  whether  this  is 
primarily  due  to  PAF,  or  secondarily  due  to  a  chain  of  events  that  have  been 
initiated  by  PAF  (3).  One  possible  explanation  arises  from  the  fact  that  PAF  has 
been  shown  to  increase  cytosolic  calcium  in  a  variety  of  cell  lines  (12, 14)  including 
canine  airway  epithelial  cells  in  primary  culture  (unpublished  data,  Man  et  al).  As 
cytosolic  Ca2+  is  implicated  in  intracellular  regulation  of  many  absorptive  or 
secretory  processes,  it  is  conceivable  that  this  agent  may  enhance  the  absorptive 
function  of  the  alveolar  and  airway  epithelial  cells.  Alveolar  protein  clearance,  on 
the  other  hand,  may  have  proceeded  via  pinocytotic  vesicles  in  alveolar  type  I  cells 
and  in  capillary  endothelial  cells  as  demonstrated  by  Bensch  et  al.  (2)  and  these 
vesicular  transport  processes  may  have  been  unaffected  by  PAF. 

Recently  it  has  been  shown  that  solute-coupled  fluid  transport  across  the 
alveolar  epithelium  plays  a  significant  role  in  the  removal  of  liquid  from  the  alveolar 
space  (6).  PAF  in  the  airspaces  appears  to  either  supplement  or  enhance  these 
active  transport  processes  in  some  manner,  but  only  transiently.  The  action  of  PAF 
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appears  to  be  relatively  short-lived  as  its  effect  was  only  observable  from  the  NMR 
data  within  ~60  minutes  of  the  introduction  of  serum  plus  PAP  to  the  airspaces. 
This  is  not  surprising  as  PAF  has  been  reported  to  be  metabolized  rapidly  by  type  II 
epithelial  cells  (14)  and  to  have  a  transient  effect  in  a  variety  of  situations  (5,9,  12). 
More  than  90  minutes  after  the  instillation  of  serum,  or  serum  plus  PAP,  the  rates 
of  alveolar  fluid  clearance  determined  from  integrated  NMR  signal  intensities  were 
approximately  equal  and  are  therefore  expected  to  have  been  driven  by  the  same 
processes  at  this  stage  of  the  clearance. 

The  ability  to  decompose  the  transverse  magnetization  decay,  for  all  cats 
with  fluid  instilled  into  their  lungs,  into  two  reproducible  components,  suggests  that 
we  were  detecting  NMR  signal  from  at  least  two  lung  water  compartments  between 
which  water  exchange  was  slow  with  respect  to  T2S-  Shioya  et  al  (28)  also 
observed  two  transverse  relaxation  times  for  lung  protons,  in-vivo.  In  normal  rat 
lungs  they  report  T2  values  of  9.5  msec  (~90%)  and  34  msec  (~10%),  where  the 
percentages  in  parentheses  denote  the  proportions  of  each  component.  By 
comparison,  our  overall  four  hour  mean  T2  values  from  serum  instilled  lungs  are  25 
msec  ±  1  msec  (n  =  24)  (59%  ±  2%)  and  177  msec  ±  7  msec  (n  =  24)  (41%  ±  2%). 
Again  the  parenthetic  percentages  denote  the  proportions.  Because  the  minimum 
inter-echo  interval  in  our  experiments  was  26.3  msec,  we  were  not  able  to  resolve 
the  faster  of  the  two  relaxation  components  observed  by  Shioya.  It  is  likely, 
therefore,  that  the  faster  relaxing  component  resolved  by  us  is  actually  a 
combination  of  the  two  components  observed  by  Shioya  in  the  normal  lung.  This 
component  is  then  assumed  to  be  due  to  tissue-bound  fluid  as  in  the  cells  and 
interstitial  space  of  the  lung.  The  more  slowly  relaxing  component  in  our  data, 
which  does  not  correspond  to  either  component  observed  by  Shioya,  could  then  be 
consistently  interpreted  as  the  instilled  fluid  residing  in  the  airspaces. 
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Additional  support  for  this  interpretation  was  observed  in  a  few  isolated 
cases  in  which  images  showed  that  some  of  the  instilled  fluid  remained  in  large 
airways.  Fluid  in  the  large  airways  was  found  to  have  T2  values  ranging  from  400 
msec  to  220  msec  as  compared  with  the  instillate  itself  which  had  a  T2  of 
approximately  480  msec.  In  contrast,  for  the  fluid  which  was  distributed  over  part 
of  a  lung  lobe,  T2  values  ranged  between  190  msec  and  130  msec,  as  discussed 
previously.  In  all  of  these  cases,  the  transverse  relaxation  in  the  lung  also  had  a 
faster  relaxing  component  with  a  T2  of  approximately  25  msec.  This  elevation  of 
the  slow  component  transverse  relaxation  time  for  fluid  in  large  airways  is 
consistent  with  an  expectation  that  relaxation  rates  (I/T2)  be  proportional  to  the 
tissue  surface  area  in  contact  with  the  airspace  fluid,  an  expectation  which  is  in  turn 
consistent  with  the  work  of  Kveder  et  al.  who  report  that  Ti  and  T2  relaxation 
within  the  lung  arises  from  a  rapid  exchange  between  bulk  water  and  water  protons 
which  are  tightly  bound  to  biopolymer  segments  (15). 

Although  there  is  little  or  no  signal  obtained  from  the  normal  lung,  as  can  be 
seen  in  figure  3a,  within  15  minutes  after  fluid  is  instilled  into  the  lung,  ~60%  of  the 
signal  corresponds  to  the  relaxation  component  that  we  have  attributed  to  tissue- 
bound  fluid  compartment.  Equally  surprising  is  the  fact  that  the  relative  signal 
proportions  corresponding  to  each  of  the  two  relaxation  components  changes  very 
little  during  the  observed  fluid  clearance,  even  with  as  much  as  60%  of  the  instilled 
fluid  leaving  the  lung  in  some  cases.  To  reconcile  these  data  one  would  have  to 
postulate  that  fluid  is  effectively  being  cleared  from  both  the  faster  relaxing 
compartment  (the  cells  and/or  interstitium),  and  the  more  slowly  relaxing 
compartment  (the  airspace  fluid),  by  means  of  a  flux  of  fluid  from  the  airspaces  into 
the  interstitial  space  of  the  lung  that  is  slow  on  a  T2  time  scale  but  fast  relative  to  the 
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overall  time  course  of  these  experiments.  This  postulate  is  supported  by  the 
observations  of  Gee  and  Staub  (11)  that  perivascular  fluid  cuffs  were  formed  within 
10  minutes  of  fluid  instillation  into  the  airspaces  of  a  dog  lung.  The 
photomicrographs  in  figures  2.6a  and  2.6b  also  lend  some  support  this  hypothesis 
since  it  is  clear  that  PAF  has  an  affect  not  only  the  fluid  clearance  from  the 
airspaces,  but  also  on  the  content  of  the  interstitial  and  vascular  spaces  of  the  lung 
as  well. 


2.5  Summary 


Results  obtained  from  in-vivo  NMR  images,  and  nuclear  relaxation  time 
measurements,  and  confirmed  gravimetrically  and  histologically  from  excised  lungs, 
demonstrate  the  unique  capabilities  of  NMR  techniques  for  studying  the  lung  fluid 
balance.  In  the  in-vivo  study,  we  were  able  to  obtain  both  quantitative  and 
qualitative  results  of  alveolar  fluid  clearance.  With  PAF  added  to  the  instillate, 
clearance  was  markedly  increased  over  a  short  30  minute  period,  but  was  bi-phasic 
overall.  Although  the  exact  mechanism  has  not  been  defined,  it  appeared  to  be  the 
result  of  a  pharmacological  action  of  PAF. 

The  resolution  of  transverse  relaxation  decay  curves  into  two  components 
has  allowed  us  to  speculate  that  tissue  bound  fluid  and  fluid  in  the  airspaces  are 
essentially  separated  on  the  short  time  scale  of  transverse  relaxation.  However,  the 
constancy  of  the  relative  proportions  of  these  two  compartments  over  the  4  hour 
observation  period  and  their  independence  of  PAF,  requires  the  additional 
hypothesis  that  exchange  of  airspace  fluid  does  exist  but  only  at  a  rate  sufficient  to 
equilibrate  populations  over  periods  of  many  minutes. 
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— In-vivo — NMR _ Assessment  of  Bi-directional  Alveolar  Fluid 

Fluxes.1 


3.1  Introduction 


The  clearance  of  excess  lung  alveolar  fluid  and  protein  has  been  studied 
under  a  variety  of  conditions  with  a  number  of  animal  species  and  in  every  case 
fluid  clearance  has  been  shown  to  proceed  more  rapidly  than  protein  clearance  (5, 

18,  24,  25).  Moreover,  alveolar  fluid  clearance  has  been  shown  to  proceed  against 
opposing  oncotic  pressure  gradients  which  rise  to  greater  than  50  cmH20  as 
proteins  become  concentrated  in  the  alveolar  space  (17).  Matthay  et  al.  proposed 
that  fluid  may  be  cleared  across  the  alveolar  epithelium  in  opposition  to  a  large 
pressure  gradient  by  means  of  active  cellular  transport  mechanisms  (17).  More 
recently,  this  proposal  has  been  supported  by  the  identification  of  several 
mechanisms  of  solute-coupled  transport  across  the  alveolar  epithelium  (3,  4,  6,  8, 

19,  24).  The  presence  of  these  active  cellular  transport  mechanisms  suggests  that 
the  alveolar  epithelium  may  be  much  more  active  than  was  previously  thought,  and 
that  there  may  be  a  certain  amount  of  fluid  transport  in  both  directions  across 
alveolar  epithelial  cells.  As  a  result,  the  goal  of  this  study  was  to  determine  if  fluid 
transport  in  both  directions  across  the  alveolar  epithelium  could  be  detected  in-vivo 
by  means  of  nuclear  magnetic  resonance  (NMR),  and  to  obtain  estimates  of  the 
relative  flows  in  each  direction. 


1  A  version  of  tills  chapter  is  to  be  submitted  for  publication.  P.  W.  Stroman,  S.  F.  P.  Man,  P . 
S.  Allen.  In-vivo  NMR  Assessment  of  Bi-directional  Alveolar  Fluid  Fluxes.  Journal  of  Applied 
Physiology. 
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Excess  alveolar  fluid  was  modeled  by  instilling  a  50%  deuterated  water 
based  serum-like  solution  into  one  lung  lobe  of  an  anesthetized  cat.  Because  there 
is  a  negligible  amount  of  2H  occurring  naturally  in  the  body,  this  model  provided  us 
with  a  fluid  tracer,  2H,  only  on  the  airspace  side  of  the  pulmonary  air-blood  barrier. 
In  order  to  monitor  the  fluid  and  tracer  clearances  from  the  airspaces,  we  applied 
proton  and  deuteron  NMR  techniques  consecutively  and  serially.  D2O  has 
previously  been  used  successfully  as  a  tracer  of  H2O  movements  (12).  We  have 
assumed  that  the  clearance  of  instilled  alveolar  2H  reflects  the  time  course  of  fluid 
absorption  across  the  pulmonary  air-blood  barrier  because  H2O,  D2O  and  HOD  are 
indistinguishable  to  the  alveolar  epithelium  and  because  in  liquid  water,  the  motion 
of  and  2H  is  governed  by  the  motion  of  intact  water  molecules  and  not  by 
chemical  exchange  (27).  The  time  course  of  the  changes  of  the  alveolar  !H  content 
is  influenced  both  by  the  absorption  of  HOD  from  the  airspaces  and  by  the 
secretion  of  H2O  into  the  airspaces,  whereas  the  clearance  of  the  tracer,  2H,  reflects 
only  the  rate  of  HOD  absorption  from  the  airspaces.  Thus,  the  combination  of 
these  two  NMR  techniques  enabled  us  to  quantify  the  rates  of  fluid  absorption  and 
secretion  across  the  pulmonary  air-blood  barrier. 

Because  we  have  demonstrated  previously  that  the  net  alveolar  fluid 
clearance  in  cats  can  be  enhanced  by  the  inclusion  of  PAF  in  the  instillate  (25),  we 
also  chose  to  assess  in  this  study  the  separate  effects  of  this  agent  on  airspace  fluid 
absorption  and  secretion.  Platelet-activating  factor  (PAF)  is  a  family  of  ether-lipids 
whose  active  component  is  l-0-Hexadecyl/octadecyl-2-acetyl-$n-glycero-3- 
phosphocholine  (22).  PAFs  are  produced  by  a  wide  variety  of  cells  in  the  body 
including  alveolar  type  II  epithelial  cells  (13,  14,  21,  22).  Moreover,  these 
molecules  may  have  physiological  and/or  pathological  significance,  as  it  has  been 
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shown  that  the  pulmonary  response  to  intra-airway  PAF  mimics  asthma,  and  PAP 
has  been  detected  in  the  lungs  of  humans  with  sarcoidosis  (2,  23). 


3.2  Methods 

3.2.1  Animal  Model  of  Alveolar  Fluid  Clearance 

Adult  mongrel  cats  (3.0  -  5.0  kg)  were  studied  after  being  anesthetized  with 
an  intramuscular  injection  of  ketamine  (Ketalean  25  mg/kg),  xylazine  hydrochloride 
(Rompun  0.5  mg/kg)  and  acepromazine  (Atravet  1  drop/4  kg),  and  anesthesia  was 
maintained  at  ~25%  of  the  initial  dose/hour.  Throughout  the  study  cats  were  in  a 
supine  position,  intubated,  and  were  ventilated  with  room  air  at  15  breaths/minute,  a 
tidal  volume  of  10  ml/kg,  and  2  cmH20  of  positive  end-expiratory  pressure 
(PEEP).  Arterial  blood  gas  tensions  were  monitored  to  assure  adequate  ventilation 
was  being  maintained. 

The  instillate  was  prepared  as  a  phosphate  buffered  physiological  saline 
solution  (PSS)  which  included  5.5  mM  D-glucose  and  equal  volumes  of  H2O  and 
D2O.  After  preparation  this  solution  was  filtered  through  a  0.22  jam  Millipore® 
filter  to  remove  any  possible  bacterial  contaminants.  Just  prior  to  instillation  into 
the  airspaces,  5  g/dL  bovine  serum  albumin  was  added  to  5  ml  of  PSS  to  make  a 
serum-like  solution  (SLS).  Also  prior  to  instillation,  1  mg  of  Evan's  Blue  dye  was 
added  to  mark  the  location  of  the  instillate  in  the  lung,  and  3  pCi  of  12^I-albumin 
was  added  to  label  the  instilled  albumin.  The  resulting  SLS  was  instilled  at  a  dose 
of  0.7  ml/kg  into  a  localized  region  of  one  lower  lung  lobe,  via  a  polyethylene 
catheter  (PE-90).  For  studies  which  included  PAF,  250  pi  of  10'5  M  PAF  was 
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also  added  to  the  SLS  immediately  prior  to  instillation;  a  dose  range  of  0.17  to  0.19 
|lg  PAF/kg  body  weight  was  used  (mean  =  0.18  fig  PAF/kg). 


Prior  to  fluid  instillation,  mechanical  ventilation  was  interrupted  and  a  PE- 90 
catheter  was  inserted  into  a  lower  lung  lobe  (left  lung,  7  times  out  of  10 
experiments).  The  SLS  was  then  injected  with  a  syringe  over  15  sec,  followed  by 
the  insufflation  of  several  ml  of  room  air.  Mechanical  ventilation  was  resumed 
immediately  after  fluid  instillation.  The  precise  amount  of  fluid  instilled  was 
determined  by  weighing  the  syringe  and  catheter  assembly  before  and  after 
instillation.  NMR  imaging  and/or  gravimetric  analysis  of  the  lung  were  then  carried 
out,  as  described  below,  during  the  succeeding  four  hour  period. 

At  the  end  of  the  4  hour  observation  period,  the  cat  was  euthanized  and  the 
lungs  were  removed.  The  distribution  of  the  instilled  fluid  in  the  experimental  lobe 
was  identified  by  the  Evan's  Blue  dye.  The  lung  lobes  were  processed  for 
gravimetric  measurements  of  the  excess  fluid,  and  for  gamma  counting  to  determine 
the  amount  of  instilled  125I-albumin  remaining  in  the  fluid-instilled  lobe,  as  in  our 
previous  study.  For  the  present  study,  however,  the  excess  alveolar  fluid  contained 
both  H2O  and  D2O. 

A  total  of  5  cats  were  studied  after  instillation  of  SLS  with  PAF  into  a  lower 
lung  lobe,  and  another  5  were  studied  without  PAF.  Gravimetric  measurements 
were  obtained  for  all  cats  but  the  changes  in  the  !H  and  2H  contents  of  the  instilled 
fluid  were  monitored  in  the  lungs  of  only  4  cats  in  each  of  these  two  groups. 
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3.2.2  !H  and  2H  NMR  Methods 

NMR  measurements  were  carried  out  on  a  modified  Bruker  CXP 
Spectrometer  (Bruker  Instruments  Inc.,  Billerica  MA)  with  a  40  cm  bore  magnet  at 
a  static  field  of  2.35  T.  Proton  data  were  acquired  using  a  surface  coil  attached  to 
the  plexiglass  tray  which  supported  the  cat  in  a  supine  position,  and  curved  around 
the  cat's  back  to  provide  maximum  sensitivity  to  the  lower  lung  lobes.  A  typical 
repetition  interval  for  proton  NMR  was  1  second.  A  second  smaller  surface  coil  (6 
turns,  8  cm  diameter)  attached  to  the  same  tray,  was  employed  for  both  transmission 
and  reception  at  the  deuteron  resonant  frequency  (15.4  MHz).  To  overcome  the 
rapid  deuteron  transverse  relaxation  in  the  lung  (T2*  ~  3  msec),  the  deuteron  NMR 
signal  was  acquired  in  the  form  of  FIDs  and  the  peak  signal  amplitude  was  sampled 
only  30  ps  after  the  RF  excitation  pulse  was  transmitted.  As  a  result,  the  2H  NMR 
signal  intensity  was  not  unduly  influenced  by  the  rapid  deuteron  transverse 
magnetization  decay  in  the  lung. 

Measurements  of  the  total  2H  NMR  signal  intensity,  without  spatial 
resolution,  were  made  with  a  technique  similar  to  that  used  by  Kim  et  al.  (12).  The 
total  signal  intensity  was  measured  from  1000  acquisitions  with  a  repetition  interval 
of  100  msec.  Our  assumption  that  the  2H  NMR  intensity  could  be  used  to  predict 
the  quantity  of  2H  within  our  region  of  interest  was  repeatedly  verified  by  a 
correlation  (r  ;>  0.993)  between  intensity  measurements  and  concentrations  of 
seven  100  ml  phantoms  ranging  from  1 .0%  to  0.05%  deuteration. 

Because  the  receiver  sensitivity  was  spatially  non-uniform,  a  change  of  the 
instillate  location  in  the  lung,  relative  to  the  receiver  coil,  could  have  a  considerable 
effect  on  the  detected  signal  intensity.  It  was  therefore  necessary  to  check  by 
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imaging  that  the  fluid  had  not  dispersed  during  clearance.  After  fluid  instillation, 
the  deuteron  signal  intensity  measurements  were  each  followed  by  a  deuteron  image 
and  a  proton  image,  and  the  sequence  repeated  throughout  the  4  hour  observation 
period.2 


Deuteron  NMR  images  were  obtained  by  applying  a  reconstruction  from 
projections  technique  to  eight  2H  signal  projections  from  a  transverse  slice.  Slice 
definition  was  provided  by  the  spatial  distribution  in  receiver  sensitivity,  a  method 
which  maximized  the  deuteron  signal  contribution  to  the  image.  The  projections 
were  generated  at  equal  angular  separations  between  0°  and  180°,  from  the 
corresponding  FID.  Employing  a  repetition  time  of  100  msec,  512  averages  were 
used  for  each  projection.  The  image  was  reconstructed  by  means  of  a  back- 
projection  algorithm  (15, 16). 

Transverse  proton  NMR  images  of  a  1  cm  thick  slice  of  the  cat's  thorax 
were  obtained  using  a  Carr-Purcell-Meiboom-Gill  (CPMG)  spin-echo  pulse 
sequence  having  an  inter-echo  interval  of  18.1  msec  and  data  acquisition  gated  to 
the  cat's  breathing  as  described  in  our  previous  study  (25).  Control  images  of 
several  slices  covering  the  caudal  half  of  the  cat's  thorax  were  obtained  prior  to  fluid 
instillation.  After  instillation,  the  region  of  the  cat's  thorax  containing  the  instilled 
fluid  was  imaged  alternately  with  2H  measurements  over  a  4  hour  period.  The  4 
hour  time  course  of  changes  in  the  content  of  the  fluid-instilled  lung  lobe  was 


2  Further  details  of  the  deuteron  NMR  techniques  employed  for  this  study  are  included  in 
Appendix  3.1  of  this  thesis. 


« 


determined  by  integrating  first-echo  image  pixel  intensities  over  the  lung  region  of 
interest.3 
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3.3  Results 


An  example  of  a  deuteron  image  (Figure  3.1a)  superimposed  on  the 
corresponding  transverse  proton  image  (Figure  3.1b)  obtained  after  instilling  fluid 
into  a  left  lower  lung  lobe  is  shown  in  Figure  3.1c.  The  time  dependence  of  the 
integrated  *H  signal  intensity  data,  Ih(0,  and  that  of  the  integrated  2H  signal 
intensity  data,  Id(0,  both  averaged  over  4  cats  and  with  data  grouped  into  40  minute 
intervals,  are  shown  in  Figure  3.2  for  the  case  in  which  the  instillate  was  SLS. 
Figure  3.3  shows  the  corresponding  data  for  an  instillate  of  SLS  plus  PAF.  The 
values  of  In(t)  and  Ig>(t)  are  expressed  relative  to  the  extrapolated  intensities 
immediately  after  instillation,  Ih(0)  and  Id(0). 


3  Transverse  relaxation  time  measurements  were  also  obtained  tor  the  same  lung  region  of  interest. 
These  relaxation  time  measurements  supplement  those  discussed  in  Chapter  2  ot  this  thesis  but  do 
not  contribute  to  the  current  discussion  of  bi-directional  fluid  movements,  and  so  are  discussed  in 


Appendix  3.2  of  this  diesis. 
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Figure  3 . 1  a)  A  NMR  image  of  a  1  cm  thick  transverse  slice 
of  a  cat’s  thorax,  30  minutes  after  the  instillation  of  SLS 
into  the  left  lower  lung  lobe  (upper  left),  b)  The  corresponding 
NMR  image  of  the  instilled  fluid  in  the  cat’s  left  lower 
lung  lobe,  as  well  as  two  of  the  projection  profiles  used  to 
construct  the  image,  acquired  20  minutes  after  fluid  instillation 
(upper  right),  c)  The  same  and  ~H  NMR  images  overlaid  for  a 
spatial  comparison  (bottom). 
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Figure  3.2  The  averaged  (n  =  4)  time  course  of  the  H2O  content  (black  circles)  and 
of  the  D2O  content  (black  squares)  of  a  lower  lung  lobe  after  instillation  of  SLS. 
The  data  points  correspond  to  averages  over  a  40  minute  data  acquisition  period 
centered  at  the  time  of  the  data  point.  The  error  bars  indicate  the  standard  error  of 
the  mean. 

Neither  2H  nor  NMR  images  demonstrated  any  change  in  the  location  of  the 
instilled  alveolar  fluid  over  the  time  course  of  an  experiment,  either  toward  or  away 
from  the  receiver  coil,  and  so  the  sensitivity  of  the  receiver  to  the  fluid  tracer  was 
assumed  to  remain  constant  during  an  experiment.  Because  the  original  instillate 
contained  equal  amounts  of  and  2H,  we  are  able  to  compute  a  net  fluid  clearance 
time  course,  lNet(0>  from  equation  3.1  using  both  sets  of  NMR  data 


lNet>  which  can  be  compared  with  the  corresponding  gravimetiic  results,  is 
illustrated  in  Figures  3.4  and  3.5  for  either  SLS  or  SLS  plus  PAF,  respectively. 
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Figure  3.3  The  averaged  (n  =  4)  time  course  of  the  *H  content  (white  circles)  and  of 
the  2H  content  (white  squares)  of  a  lower  lung  lobe  after  instillation  of  SLS  plus 
PAF.  The  data  points  correspond  to  averages  over  a  40  minute  data  acquisition 
period  centered  at  the  time  of  the  data  point.  The  error  bars  indicate  the  standard 
error  of  the  mean. 

For  the  SLS  instillate,  NMR  intensity  data  showed  that  the  proton 
component  of  lung  water  increased  to  121%  ±  19%  (n=4)  of  its  starting  value  in  the 
experimental  lung  lobe  over  a  period  of  3.3  to  4  hours  after  instillation,  whereas  the 
2H  NMR  intensity  data  showed  a  decrease  in  the  deuterated  component  to  45%  ± 
13%  (n=4)  of  its  starting  value  over  the  same  time  period  (Figure  3.2).  The  total 
amount  of  fluid  remaining  in  the  experimental  lung  lobe,  averaged  over  the  final 
data  acquisition  period  and  computed  by  combining  the  and  intensity  data,  is 
83%  ±  16%  (n  =  4).  This  just  agrees  within  the  range  of  experimental  uncertainty 
with  the  corresponding  gravimetric  measurement  of  56%  ±  13%  (n  =  5)  taken  at  4 
hours  after  instillation  (Figure  3.4).  All  of  the  results  are  expressed  as  the  mean  of 
n  experiments,  plus  or  minus  the  standard  error  of  the  mean. 
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Figure  3.4  The  net  clearance  time  course  of  instilled  alveolar  SLS,  computed  from 
in-vivo  NMR  measurements  (black  circles)  and  the  corresponding  gravimetric 
measurement  (black  triangle),  with  error  bars  indicating  the  standard  error  of  the 
mean.  The  NMR  measurements  are  averages  over  40  minute  data  acquisition 
periods  and  over  4  cats.  The  gravimetric  measurement  is  an  average  over  5  cats 
sacrificed  at  4  hours  post  instillation. 


With  PAF  added  to  the  SLS,  the  lH  NMR  intensity  data  showed  no 
increase  in  the  protonated  component  remaining  in  the  experimental  lung  lobe 
between  3.3  and  4  hours  after  instillation,  but  an  essentially  constant  value  of  101% 
±  8%  (n=4)  (Figure  3.3).  The  2H  NMR  intensity  data  however,  showed  a  greater 
decrease  in  the  deuterated  component  than  without  PAF,  since  only  39%  ±  3% 
(n=4)  could  be  detected  in  the  experimental  lung  lobe  between  3.3  and  4  hours  after 
instillation.  Together  these  data  indicate  that  the  total  averaged  amount  of  instilled 
alveolar  SLS  plus  PAF  remaining  in  the  lung  lobe  over  the  final  data  acquisition 
period  was  70%  ±  6%  (n  —  4)  of  the  amount  initially  instilled.  Again  the  NMR  data 
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point  just  overlaps  within  the  limits  of  error  with  the  corresponding  gravimetric 
measurement  of  49%  ±  15%  (n  =  5)  (Figure  3.5). 


Figure  3.5  The  time  course  of  the  net  clearance  of  instilled  alveolar  SLS  plus  PAF, 
computed  from  in-vivo  NMR  measurements  (white  circles)  and  the  corresponding 
gravimetric  measurement  (white  triangle),  with  error  bars  indicating  the  standard 
error  of  the  mean.  The  NMR  measurements  are  averages  over  40  minute  data 
acquisition  periods  and  over  4  cats.  The  gravimetric  measurement  is  an  average 
over  5  cats  sacrificed  at  4  hours  post  instillation. 


The  proportions  of  125I-albumin  remaining  in  the  lung  4  hours  after 
instillation  of  either  SLS  or  SLS  plus  PAF  were  82%  ±  6%  (n  =  5)  and  86%  ±  8% 
(n  =  5),  respectively,  showing  a  marked  consistency  within  the  limits  of  their  errors. 
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3.4  Discussion 


The  net  clearance  of  instilled  alveolar  fluid  is  the  result  of  considerable  fluid 
transport  in  both  directions  across  the  pulmonary  air-blood  barrier,  with  the  rate  of 
fluid  absorption  exceeding  that  of  fluid  secretion.  Serially  alternated  in-vivo  proton 
and  deuteron  intensity  measurements  of  the  clearance  of  a  partially  deuterated 
instillate  from  cat  lung  show  a  distinct  difference  between  the  proton  and  deuteron 
clearance  rates.  The  protonated  component  clears  more  slowly  (presumably  due  to 
H2O  secretion  into  the  airspaces)  thus  supporting  the  premise  that  bi-directional 
fluid  movement  across  the  air-blood  barrier  can  be  detected  using  a  combined 
!H/2H  monitoring  technique.  More  conclusive  support  for  this  premise  can  be 
found  in  the  agreement  of  (see  Fig  3.6)  the  combined  proton  and  deuteron  clearance 
of  the  present  study  with  the  net  clearance  of  H2O  observed  previously  for  a  fully 
protonated  autologous  serum  instillate  (25).  Matthay  et  al.  (17)  similarly  found  that 
the  clearance  of  instilled  Ringer's  lactate  with  5%  albumin  was  identical  to  that  of 
autologous  serum.  Accepting  the  premise,  the  individual  rates  of  absorption  and 
secretion  can  be  obtained  from  the  time  course  of  the  deuteron  and  proton  signals. 
For  example,  for  an  instillate  that  is  50%  deuterated,  the  time  dependence  of 
absorption,  A(t),  and  secretion,  S(t),  expressed  as  a  percentage  of  the  original 
instillate,  are4 

(3-2)  A(0=[l-K]*l°° 


4  The  derivation  of  equations  3.2  and  3.3  is  provided  in  Appendix  3.3  of  this  thesis. 
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and  they  are  illustrated  in  Figures  3.7  and  3.8  respectively. 


a> 

u 

c 

03 

Ih 

03 

<U 

0 

a 

E 

4> 

z 


0  60  120  180  240  300 

Time  Since  Instillation  (minutes) 

Figure  3.6  A  comparison  of  the  time  course  of  the  net  clearance  from  the  present 
experiments  and  from  the  previous  work  (25)  on  net  *H  clearance.  Squares  are 
used  to  indicate  values  computed  from  equation  3.1  from  the  present  combined 
and  2H  intensity  measurements,  whereas  the  circles  represent  previous  results  of 
intensity  measurements  (25).  The  instillates  are  represented  as  follows,  SLS 
(black  squares),  serum  (black  circles),  SLS  plus  PAF  (white  squares),  and  serum 
plus  PAF  (white  circles). 
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Figure  3.7  An  illustration  of  the  time  course  of  the  amount  of  fluid  absorbed  from 
the  lung,  relative  to  the  amount  instilled,  calculated  from  equation  3.2.  Absorption 
of  SLS  (black  circles)  and  of  SLS  plus  PAF  (white  circles)  are  shown. 
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Figure  3.8  An  illustration  of  the  time  course  ol  the  amount  ol  lluid  secreted  into  the 
lung,  relative  to  the  amount  instilled,  calculated  trom  equation  3.3.  Fluid  secretion 
after  the  instillation  of  SLS  (black  circles),  and  after  the  instillation  of  SLS  plus 
PAF  (white  circles),  are  shown. 
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Figures  3.6,  3.7  and  3.8,  indicate  that  over  the  4  hour  time  period  in  which 
~20%  of  instilled  SLS  is  cleared  from  the  lung,  -55%  of  the  instilled  fluid  volume 
is  absorbed  from  the  airspaces,  while  a  volume  approaching  35%  of  that  instilled  is 
secreted  back  into  the  airspaces.  This  result  illustrates  a  considerable  amount  of 
fluid  moving  in  both  directions  across  the  air-blood  barrier  of  the  lung  during  the 
net  clearance  of  instilled  alveolar  fluid.  It  is  consistent  with  a  model  in  which  the 
net  clearance  is  primarily  the  result  of  solute-coupled  ion  transport,  with  passive 
fluid  movements  being  of  secondary  importance,  as  has  been  suggested  by  Matthay 
and  by  Basset  et  al.  (4, 19).  The  efficacy  of  the  solute-coupled  Na+  transport  across 
the  alveolar  epithelium  has  also  been  clearly  demonstrated  by  Berthiaume  et  al.  and 
by  Smedira  et  al.  (6, 24).  Consequently,  although  there  may  be  a  certain  amount  of 
passive  fluid  movement  via  paracellular  pathways  contributing  to  the  observed  fluid 
secretion  into  the  airspaces  (17),  the  impermeable  nature  of  the  alveolar  epithelium, 
the  existence  of  a  Na+-  K+  pump  on  the  basolateral  side  of  type  II  epithelial  cells 
(4),  and  the  demonstrated  Cl'  secretion  by  tracheal  epithelial  cells  (26),  lead  us  to 
postulate  that  the  primary  source  of  secreted  alveolar  fluid  is  also  solute-coupled 
transport. 

The  effects  of  PAF  on  the  rates  of  net  fluid  clearance,  and  fluid  absorption 
and  secretion  across  the  pulmonary  air-blood  barrier,  are  also  illustrated  by  Figures 
3.6,  3.7  and  3.8.  For  example,  during  the  first  hour  after  fluid  instillation, 
approximately  25%  of  instilled  SLS  is  absorbed  from  the  airspaces,  whereas  with 
PAF  added  to  the  SLS,  the  amount  absorbed  is  increased  to  approximately  55%  of 
that  instilled.  The  net  clearance  is  increased  from  -10%  to  -30%  by  the  presence  of 
PAF,  and  over  the  whole  4  hour  period,  an  analysis  of  variance  (ANOVA)  shows 
that  the  whole  time  course  for  SLS  clearance  differs  significantly  fiom  that  of  SLS 
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plus  PAP  at  the  95%  confidence  level.  This  enhancement  of  the  net  clearance  by 
the  presence  of  PAF  in  the  instillate  was  also  found  previously  (25)  for  autologous 
serum  instillates.  However,  the  enhancement  was  more  marked  in  the  case  of 
autologous  serum,  probably  due  to  the  additional  components  in  blood  serum  (7, 
10) 


The  more  detailed  look  at  the  clearance,  which  is  afforded  only  by  the 
combined  1H/2H  technique,  enables  the  effect  of  PAF  on  the  components  of  the  bi¬ 
directional  flow  to  be  observed.  Figures  3.7  and  3.8  show  that  while  PAF 
enhances  clearance,  it  apparently  has  little  effect  upon  secretion.  One  possible 
mechanism  of  PAF  action  is  suggested  by  the  fact  that  this  agent  has  been  shown  to 
increase  cytosolic  calcium  in  a  variety  of  cell  lines  (1,  11,  14,  20)  including  canine 
airway  epithelial  cells  in  primary  culture  (Man  et  al,  unpublished  data).  Because 
cytosolic  Ca2+  is  implicated  in  intracellular  regulation  of  many  absorptive  or 
secretory  processes,  it  is  conceivable  that  this  agent  may  enhance  the  absorptive 
function  of  the  alveolar  and  airway  epithelial  cells.  An  alternative  possibility  for  the 
mechanism  of  PAF  action  is  suggested  by  the  fact  that  PAF  has  been  shown  to 
increase  intracellular  cAMP  (cyclic  3',5'-adenosine  monophosphate)  in  human 
lymphocytes  (28).  Agents  which  act  to  increase  intracellular  cAMP  have  been 
shown  to  result  in  increased  ion  transport  across  primary  cultured  monolayers  of 
alveolar  type  II  epithelial  cells  (9).  Moreover,  Berthiaume  has  demonstrated  that 
increasing  lung  tissue  cAMP  levels  in-vivo  results  in  an  increased  lung  fluid 
clearance  (6).  Thus,  a  plausible  explanation  of  the  observed  PAF-enhanced  alveolar 
fluid  absorption  will  be  provided  if  the  action  ol  PAF  can  be  shown  to  increase 
intracellular  cAMP  levels  in  pulmonary  epithelial  cells  as  well. 
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The  (8%  ±  6%)  clearance  of  protein  from  the  alveolar  space  over  a  4  hour 
period  was  indicative  of  a  much  slower  clearance  rate  than  that  of  the  alveolar  fluid 
itself.  Moreover,  protein  clearance  was  apparently  unaffected  by  PAP.  both  of 
these  findings  are  in  agreement  with  the  results  of  our  previous  study  on  the  net 
clearance  of  autologous  serum  (25)  showing  that  neither  SLS  nor  SLS  plus  PAP 
results  in  a  change  in  the  sieve  properties  of  the  epithelium  when  it  is  introduced 
into  the  airspaces  of  the  lung. 


3.5  Summary 


By  means  of  a  combination  of  *H  and  2H  NMR  techniques  we  have  been 
able  to  quantify  fluid  movements  in  opposite  directions  across  the  pulmonary  air- 
blood  barrier,  in  a  living  animal.  To  our  knowledge  this  is  the  first  such 
observation.  This  bi-directional  movement  of  fluid  is  thought  to  be  primarily  the 
result  of  solute-coupled  ion  transport,  with  passive  fluid  movements  being  of 
secondary  importance.  Furthermore,  by  obtaining  independent  measures  of 
absorption  and  secretion,  we  have  been  able  to  show  that  the  means  by  which  PAF 
enhances  alveolar  fluid  clearance  in  the  first  hour  after  fluid  instillation  is  by 
increasing  the  rate  of  fluid  absorption,  while  having  no  apparent  effect  on  the  rate  of 


fluid  secretion. 
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4.0  General  Discussion  and  Veracity  of  the  Conclusions 

Light  microscopy  of  the  lung  tissue  4  hours  after  the  instillation  of  either 
serum  or  serum  plus  PAF  showed  no  evidence  of  inflammation,  and  that  the 
alveolar  epithelium  and  capillary  endothelium  were  apparently  intact  at  that  time. 
We  can  postulate  that  PAF  did  not  cause  a  change  in  the  sieve  properties  of  the 
epithelium  because  the  clearance  of  albumin  from  the  alveolar  space  was 
observed  to  be  ~3%/hour  irrespective  of  whether  the  instillate  was  serum,  serum 
plus  PAF,  SLS,  or  SLS  plus  PAF.  Moreover,  we  can  be  fairly  certain  that  the 
cat's  overall  pulmonary  function  was  being  maintained  at  a  normal  level,  and  was 
not  disrupted  by  the  instillation  of  a  localized  bolus  of  fluid,  because  blood  gas 
tensions  remained  relatively  stable  in  these  anesthetized,  spontaneously  breathing 
cats.  Our  observations  indicate  that  fluid  clearances,  both  with  and  without  PAF, 
were  from  relatively  normal,  intact  lung  lobes  and  may  thus  have  direct 
physiological  relevance. 


4.1  NMR  and  Gravimetric  Measurements  of  Clearance  Rates 

The  NMR  studies  discussed  in  Chapters  2  and  3  have  demonstrated  that 
non-invasively ,  in  a  living  animal,  we  are  able  to  monitor  the  net  clearance  of 
instilled  alveolar  fluid,  and  to  quantify  the  absorption  and  secretion  components 
of  this  net  clearance.  The  net  clearance  (17%/4  hours)  of  a  serum  (or  serum-like) 
solution  was  observed  to  proceed  more  rapidly  than  that  of  albumin  (8%/4  hours), 
and  was  the  result  of  fluid  absorption  and  secretion  at  55  /olA  hours  and  38%/4 
hours,  respectively.  These  results  are  consistent  with  a  model  in  which  the  net 
clearance  of  alveolar  fluid  is  primarily  the  result  of  solute-coupled  ion  transport. 
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with  passive  fluid  movements  being  of  secondary  importance,  as  was  initially 
proposed  by  Matthay  (11). 

After  instilling  3.1  ml  of  SLS  into  a  3  kg  cat,  the  computed  rates  of 
absorption  and  secretion  were  approximately  80  nl/sec  and  55  nl/sec,  respectively. 
The  action  of  PAF  was  observed  to  increase  the  SLS  absorption  rate  to 
approximately  315  nl/sec,  but  only  over  the  first  hour  after  instillation,  and  had  no 
discernible  effect  on  the  fluid  secretion  rate  (see  Chapter  3).  In  combination  with 
an  autologous  serum  instillate  the  action  of  PAF  was  even  more  pronounced  than 
with  SLS,  as  demonstrated  by  the  experimentally  observed  clearance  rates  of 
65%/4  hours,  and  30%/4  hours,  respectively  (see  Chapter  2).  In  a  3  kg  cat,  these 
net  clearances  are  equivalent  to  average  4  hour  clearance  rates  of  95  nl/sec  and  44 
nl/sec,  respectively.  We  have  proposed  that  PAF  may  act  by  increasing  the 
cytosolic  Ca2+  in  pulmonary  epithelial  cells  to  enhance  the  absorptive  function  of 
the  epithelium,  or  by  increasing  the  levels  of  intracellular  cAMP,  thereby 
stimulating  epithelial  ion  transport.  Regardless  of  the  mechanism  of  PAF  action 
however,  the  observed  clearance  time  courses  of  serum  plus  PAF,  and  of  SLS  plus 
PAF,  indicate  that  there  are  one  or  more  components  of  serum,  not  included  in  the 
SLS  preparation,  that  act  to  enhance  the  effect  of  PAF  on  alveolar  fluid  clearance. 


4.2  Transverse  Relaxation  Time  Measurements  of  Water  Protons 

The  transverse  relaxation  time  measurements  discussed  in  Chapter  2  and 
Appendix  3.2  have  consistently  demonstrated  that  we  are  detecting  signal  from 
two  lung  water  compartments,  between  which,  fluid  exchange  is  slow  with  respect 
to  the  measured  transverse  relaxation  times.  In  Chapter  2  we  have  proposed  that 
the  two  relaxation  components,  which  are  identified  by  the  relaxation  times  T2F 
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and  T2S  (faster  relaxing  and  slower  relaxing,  respectively),  are  associated  with 
tissue-bound  fluid  in  the  interstitial  and  intracellular  spaces,  and  with  airspace 
fluid,  respectively.  The  measured  values  of  T2F  were  consistently  in  the  range  of 
25  msec  to  32  msec  and  demonstrated  no  systematic  variations  over  4  hours  of 
fluid  clearance  from  the  lung.  Similarly,  measured  values  of  T2S  did  not 
demonstrate  conclusive  systematic  variation  over  4  hours  of  fluid  clearance,  but 
the  4  hour  mean  values  of  T2S  varied  considerably  with  the  fluid  instilled  (serum- 
based  or  SLS-based)  and  with  the  mode  of  ventilation  employed  (spontaneous 
breathing  or  mechanical  ventilation  with  positive  end-expiratory  pressure). 
Averaged  values  of  T2S  ranged  from  150  msec  to  544  msec  and  strongly  support 
our  association  of  T2S  with  airspace  fluid.  Moreover,  the  dependence  of  T2S  on 
the  mode  of  ventilation  employed  supports  the  conclusions  of  Kveder  et  al.  (8), 
that  relaxation  rates  are  proportional  to  the  tissue  surface  area  in  contact  with  the 
airspace  fluid,  because  water  proton  relaxation  within  the  lung  arises  from  a  rapid 
exchange  between  bulk  water  and  water  protons  which  are  tightly  bound  to 
biopolymer  segments. 

One  aspect  of  the  observed  transverse  relaxation  components  that  has  not 
been  discussed  in  Chapter  2  or  Appendix  3.2,  is  the  relative  amount  of  signal,  pp 
or  ps,  contributing  to  either  the  faster  relaxing  or  the  slower  relaxing  component. 
Although  Figure  2.3a  demonstrates  that  there  is  little  or  no  proton  NMR  signal 
obtained  from  the  normal  lung,  within  15  minutes  after  fluid  instillation,  60%  to 
80%  of  the  NMR  signal  corresponded  to  the  relaxation  component  that  we  have 
associated  with  intracellular  and  interstitial  fluid.  To  investigate  this  surprising 
result  we  must  keep  in  mind  that  virtually  all  of  the  naturally  occurring  lung  fluid 
is  in  the  intracellular  and  interstitial  spaces,  and  also  conti ibuted  to  our  relaxation 
time  measurements.  Moreover,  the  amount  of  naturally  occurring  lung  fluid  in  a 
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cat  was  estimated  from  gravimetric  measurements  of  43  normal  lung  lobes  to  be 
4.0  g  H2O  per  kg  of  body  weight.  The  instillate  was  delivered  to  an  estimated 
14%  of  the  total  lung,  however,  and  this  portion  of  the  normal  lung  contains 
roughly  0.56  g  H20/kg  body  weight.  In  comparison,  the  instillate  itself  was 
delivered  at  a  dose  of  0.7  ml/kg,  a  dose  which  resulted  in  0.68  g  H20/kg  body 
weight  in  the  airspaces.  Based  on  these  values,  the  value  of  pp  immediately  after 
fluid  instillation  is  estimated  to  be 

PF  =  (0.68  g  H 20/kg  -T O^^g  H20/kg)  =045 

This  estimate  suggests  that  because  fluid  was  being  cleared  from  the  airspaces  and 
passing  through  the  interstitium,  the  minimum  proportion  of  the  NMR  signal  that 
could  be  expected  to  come  from  the  tissue-bound  fluid  was  approximately  45%. 
Support  for  our  association  of  T2F  with  the  tissue-bound  fluid  is  provided  by  this 
estimate,  but  it  does  not  explain  why  pp  was  consistently  observed  to  be  as  high 
as  60%  to  80%. 

As  surprising  as  the  seemingly  high  observed  values  of  pp,  is  the  fact  that 
the  values  of  pp  and  ps  changed  very  little  during  the  observed  fluid  clearance, 
even  with  as  much  as  60%  of  the  instilled  fluid  leaving  the  lung  in  some  cases. 
To  reconcile  both  of  these  features  of  our  data  at  once,  however,  one  can  postulate 
that  fluid  was  effectively  being  cleared  from  both  the  faster  relaxing  compartment 
(the  cells  and/or  interstitium),  and  the  more  slowly  relaxing  compartment  (the 
airspace  fluid),  by  means  of  a  flux  of  fluid  from  the  airspaces  into  the  interstitial 
space  of  the  lung.  Moreover,  this  fluid  flux  must  be  slow  on  a  T2  time  scale  but 
fast  relative  to  the  4  hour  time  course  of  these  experiments.  This  postulate  is 
supported  by  the  observations  of  Gee  and  Staub  (6)  that  perivascular  fluid  cuffs 
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were  formed  within  10  minutes  of  fluid  instillation  into  the  airspaces  of  a  dog 
lung  by  means  of  a  relatively  rapid  initial  fluid  flux  from  the  airspaces  to  the 
interstitium. 


4.3  A  Model  of  the  Lung  Fluid  Balance 

Using  all  of  our  NMR  observations,  we  are  now  able  to  construct  a 
detailed  postulate  of  the  mechanisms  responsible  for  the  clearance  of  instilled 
alveolar  fluid,  and,  by  inference,  for  maintaining  the  alveolar  fluid  balance.  We 
can  assume  that  the  observed  clearance  of  125I-albumin  from  the  airspaces  of  the 
lung  was  likely  via  pinocytotic  vesicles  as  demonstrated  by  Bensch  et  al  (3).  The 
fluid  clearance,  however,  is  governed  by  several  mechanisms  which  include  a 
passive  mechanism  (dominated  by  hydrostatic  and  osmotic  forces  with  fluid 
movements  via  paracellular  pathways),  and  active  cellular  transport.  Passive  fluid 
movements  are  responsible  for  the  formation  of  perivascular  fluid  cuffs  within 
~10  minutes  after  the  instillation  of  fluid  as  described  by  Gee  et  al  (6)  and  as 
indicated  by  transverse  relaxation  component  populations.  Also,  capillary  blood 
hydrostatic  pressures  may  or  may  not  have  been  modified  by  the  instillation  of 
fluid  into  the  airspaces  of  the  lung.  The  fact  that  alveolar  fluid  protein  osmotic 
pressures  are  highly  elevated  by  differential  alveolar  fluid  and  protein  clearances, 
on  the  other  hand,  has  been  well  demonstrated  (11).  Nonetheless,  active  fluid 
transport  mechanisms  are  expected  to  be  the  primary  sources  of  lung  fluid 
absorption  and  secretion,  as  well  as  being  the  mechanisms  which  are  stimulated  as 
a  result  of  PAF  action  to  enhance  airspace  fluid  absorption. 


In  order  to  estimate  the  relative  contributions  of  each  of  the 
aforementioned  mechanisms  to  the  maintenance  of  the  lung  fluid  balance,  and  to 
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the  clearance  of  instilled  alveolar  fluid,  we  have  employed  a  mathematical  model 
of  the  lung  fluid  balance.  The  basis  of  the  model  is  that  described  by  Blake  (4) 
which  describes  the  effects  of  hydrostatic  and  osmotic  forces  on  fluid  fluxes 
across  the  alveolar  epithelium  and  across  the  capillary  endothelium.  Moreover, 
the  parameters  used  in  the  model  are  to  simulate  one  lung  lobe  of  a  3  kg  cat,  in 
which  an  estimated  alveolar  surface  area  of  5000  cm^  is  in  contact  with  2.1  ml  of 
serum.  Additional  parameters  used  in  the  development  of  this  model  are  listed  in 
table  4.1  and  a  flow  chart  which  details  the  computer  algorithm  is  included  in 
Appendix  4.1. 

The  parameters  used  to  describe  the  cat's  lung  fluid  balance  were  first 
examined  with  a  model  of  the  normal  lung  which  included  only  hydrostatic  and 
osmotic  forces.  This  model  demonstrated  fluid  movements  only  from  the  blood 
into  the  interstitial  space,  and  the  rate  at  which  lymphatics  remove  fluid  from  the 
interstitium  was  adjusted  to  maintain  the  lung  in  a  steady  state  with  the  interstitial 
fluid  volume  listed  in  table  4.1.  The  number  of  endothelial  pores  per  unit  area 
was  also  adjusted  to  maintain  the  lymph/plasma  protein  concentration  ratio  near 
the  expected  value  of  0.7  (12).  Fluid  movements  across  the  alveolar  epithelium 
were  negligible  in  this  model  of  the  normal  lung.  Moreover,  altering  the  model  to 
include  serum  in  the  airspaces  had  a  negligible  effect  on  computed  trans-epithelial 
fluid  fluxes.  Thus,  this  initial  model  of  the  lung  fluid  balance  was  unable  to 
describe  the  lung  fluid  clearances  that  we  have  observed. 


144 


Epithelial  pore  diameter 

0.8  nm  (12) 

Number  of  epithelial  pores 

2.0e+10  pores/cm2 

Endothelial  pore  sizes 

2  nm  and  12.5  nm  (4) 

Ratio  of  2  nm  pores/ 12. 5  nm  pores 

1000/1  (4) 

Total  number  of  endothelial  pores 

Alveolar  diameter 

2.0e+10  pores/cm2 

100  pm  (14) 

Distance  across  cell  membranes 

0.5  pm  (13) 

Diameter  of  albumin  molecule 

3.4  nm  (4) 

Diameter  of  water  molecule 

0.15  nm  (4) 

Temperature 

311.7K 

Lymph  flow  rate 

1 .44e-03  ml/sec 

Before  fluid  instillation: 

Alveolar  fluid  volume 

0.100  ml 

Interstitial  fluid  volume 

0.25  ml 

Relative  airspace  pressure 

0  mmHg 

Relative  interstitial  pressure 

-4  mmHg  (9) 

Relative  capillary  pressure 

7  mmHg  (10) 

Interstitial  protein  concentration 

5.65e-07  moles/ml  (12) 

Interstitial  oncotic  pressure 

12  mmHg  (12) 

Capillary  protein  concentration 

7.97e-07  moles/ml  (12) 

Capillary  oncotic  pressure 

After  instillation  of  2.10  ml  of  fluid: 

20  mmHg  (12) 

Airspace  protein  concentration 

7.97e-07  moles/ml 

Airspace  oncotic  pressure 

20  mmHg 

Table  4.1  Parameters  used  in  a  mathematical  model  of  the  alveolar  fluid 
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Figure  4.1  The  time  courses  of  fluid  clearance  after  the  instillation  of  serum 
(black  circles)  and  SLS  (white  circles)  and  the  result  of  a  mathematical  model  of 
the  lung  fluid  clearance  (solid  line). 

In  order  to  simulate  the  observed  clearance  of  SLS,  and  the  corresponding 
time  courses  of  changes  in  the  lung's  H2O  and  D2O  contents,  it  was  necessary  to 
include  in  the  model  an  epithelial  fluid  absorption  rate  of  175  nl/sec  and  a  fluid 
secretion  rate  of  141  nl/sec.  For  the  purposes  of  this  model  it  was  assumed  that 
fluid  movements  across  the  endothelium  were  due  to  passive  forces  via 
paracellular  pathways  only,  because  the  endothelial  permeability  is  much  higher 
than  that  of  the  alveolar  epithelium  (13).  With  the  addition  of  active  fluid  fluxes 
across  the  epithelium,  we  were  able  to  generate  the  computed  time  course  of  lung 
fluid  clearance  shown  in  figure  4.1.  The  simulated  fluid  clearance  is  displayed  in 
comparison  to  the  experimentally  observed  fluid  clearance  of  serum  and  of  SLS. 
Also,  figure  4.2  illustrates  the  simulated  time  courses  of  the  lung's  H20  and  D20 
contents  after  the  instillation  of  SLS,  in  comparison  to  the  corresponding 
experimentally  observed  time  courses  of  the  lung's  H2O  and  D2O  contents. 
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Figure  4.2  The  time  courses  of  the  changes  of  the  H2O  and  D2O  contents  of  the 
lung  after  instillation  of  SLS.  The  corresponding  time  courses  simulated  with  a 
mathematical  model  of  the  lung  fluid  clearance  are  indicated  by  the  dotted  line 
and  the  dashed  line,  respectively.  Results  of  the  mathematical  model  are  plotted 
at  15  minute  intervals. 

The  active  fluid  absorption  and  secretion  rates  used  to  simulate  the 
alveolar  epithelium  in  the  model  were  initially  chosen  to  be  the  4  hour  average 
rates  observed  experimentally  for  the  entire  air-blood  barrier,  80  nl/sec  and  55 
nl/sec,  respectively.  These  values  were  then  adjusted  in  the  model  to  be  175 
nl/sec  and  141  nl/sec,  respectively,  to  minimize  the  sum  of  the  squares  of  the 
differences  between  our  experimental  measurements  and  the  simulated  amounts 
of  H2O  and  D2O  in  the  lung,  at  each  of  the  6  different  measurement  points  in  time 
after  fluid  was  instilled.  This  mathematical  model  is  thus  sufficiently  constrained 
that  the  estimated  rates  of  absorption  and  secretion  required  to  simulate  our 
experimental  observations  are  expected  to  be  unique  solutions.  Furthermore,  the 
simulated  fluid  clearance  depends  very  little  on  the  parameters  we  have  assumed 
to  describe  passive  fluid  movements  in  the  lung  (listed  in  Table  4.1),  and  supports 
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our  conclusion  that  alveolar  fluid  clearance  is  dominated  by  active  transport 
mechanisms.  As  a  result,  the  simulated  clearance  depends  primarily  on  the 
epithelial  absorption  and  secretion  rates  used  in  the  model  and  we  can  therefore 
assume  these  to  be  reasonable  estimates  of  those  in  the  cat  lung.  This  assumption 
is  supported  by  the  fact  that  Basset  et  al.  (1)  have  observed  an  alveolar  absorption 
rate  of  134  nl/sec  from  excised  rat  lungs  having  an  estimated  alveolar  surface  area 
of  5000  cm2,  in  comparison  to  our  estimated  absorption  rate  of  175  nl/sec  for  the 
same  alveolar  surface  area. 


Figure  4.3  The  time  courses  of  fluid  clearance  after  the  instillation  of  serum  plus 
PAF  (black  circles)  and  SLS  plus  PAF  (white  circles)  and  the  results  of  a 
mathematical  model  of  the  clearance  of  serum  plus  PAF  (solid  line)  and  the 
clearance  of  SLS  plus  PAF  (dashed  line)  from  the  airspaces  of  the  lung.  Results 
of  the  mathematical  model  are  plotted  at  15  minute  intervals. 


Simulations  of  the  observed  fluid  clearance  from  the  lung  after  the 
instillation  of  SLS  plus  PAF,  were  generated  by  introducing  supplemental 
epithelial  fluid  absorption  and  secretion  at  rates  of  1220  nl/sec  and  887  nl/sec. 


* 
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respectively.  Also,  it  was  necessary  that  these  supplemental  absorption  and 
secretion  rates  resulting  from  the  addition  of  PAF  to  the  instillate,  decay 
exponentially  with  a  half  time  of  8.3  minutes  measured  from  the  time  of  fluid 
instillation  into  the  lung.  Because  bi-directional  fluid  movements  could  not  be 
monitored  with  serum-based  instillates,  to  simulate  the  observed  clearance  of 
serum  plus  PAF  it  was  necessary  to  assume  the  same  supplemental  epithelial 
absorption  and  secretion  rates  as  those  introduced  to  model  the  clearance  of  SLS 
plus  PAF,  but  the  duration  of  PAF  action  was  assumed  to  be  greater  for  serum- 
based  instillates. 


Figure  4.4  The  time  courses  of  the  changes  of  the  H2O  and  D2O  contents  of  the 
lung  after  instillation  of  SLS  plus  PAF.  The  corresponding  time  courses 
simulated  with  a  mathematical  model  of  the  lung  fluid  clearance  are  indicated  by 
the  dotted  line  and  the  dashed  line,  respectively.  Results  ot  the  mathematical 
model  are  plotted  at  15  minute  intervals. 


To  model  the  clearance  of  serum  plus  PAF  then,  the  decay  half-time  of  the 
supplemental  active  fluid  fluxes  was  increased  to  33.3  minutes.  The  simulated 
clearances  of  serum  plus  PAF,  and  of  SLS  plus  PAF,  are  shown  in  figure  4.3  in 
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comparison  to  the  corresponding  experimentally  observed  clearance  time  courses 
of  these  solutions.  Also,  the  time  courses  of  the  experimentally  observed  H2O 
and  D2O  contents  of  the  lung  after  the  instillation  of  SLS  plus  PAP  are  shown  in 
figure  4.4  in  comparison  to  the  corresponding  time  courses  generated  with  the 
mathematical  model. 

The  simulated  clearances  of  either  serum  plus  PAF,  or  SLS  plus  PAP,  that 
we  have  generated  may  not  be  unique  because  we  cannot  be  certain  of  the  effects 
of  PAF  on  the  alveolar  epithelium.  However,  we  can  be  certain  that  the 
mechanisms  we  have  included  in  the  model  to  simulate  the  action  of  PAP  are  at 
least  plausible.  We  have  observed  that  the  absorption  of  fluid  from  the  alveolar 
space  is  enhanced  with  PAP,  and  it  has  been  shown  by  others  that  PAP  is  rapidly 
metabolized  by  alveolar  type  II  epithelial  cells  (7)  so  that  the  expected  duration  of 
PAF  action  in  the  airspaces  is  limited.  Surprisingly  though,  we  have  observed 
experimentally  that  PAF  only  enhances  the  rate  of  fluid  absorption  from  the 
airspaces,  and  yet  we  have  found  it  necessary  to  increase  the  rates  of  both 
epithelial  fluid  absorption  and  secretion  in  order  to  simulate  the  observed  fluid 
clearances  with  PAF.  The  difference  between  our  mathematical  model  and  our 
experimental  observations  lies  in  the  fact  that  experimentally  we  observed  the 
effects  of  fluid  movements  across  the  entire  air-blood  barrier,  and  experimentally 
we  were  unable  to  observe  separately  the  fluid  movements  across  the  alveolar 
epithelium  and  the  capillary  endothelium.  Nonetheless,  we  have  succeeded  in 
simulating  our  experimental  observations,  indicating  that  the  temporary 
enhancement,  or  augmentation,  of  both  epithelial  absorption  and  secretion 
mechanisms  is  a  plausible  explanation  of  the  effects  of  PAF  on  the  epithelial  fluid 
transport  mechanisms.  Also,  this  model  of  the  lung  fluid  balance  has 
demonstrated  that  the  clearance  of  fluid  from  the  alveolar  space  of  a  relatively 
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normal  lung  is  the  result  of  active  fluid  transport  mechanisms  with  the  effects  of 
hydrostatic  and  osmotic  forces  being  of  little  consequence,  as  was  first 
demonstrated  by  Basset  et  al.  (1). 

The  overall  picture  of  the  lung  fluid  balance  that  unfolds  from  our 
experimental  observations,  and  mathematical  modeling  of  the  lung,  is  as  follows. 
The  balance  of  fluid  in  the  normal  lung  is  dominated  by  active  fluid  transport 
mechanisms.  Moreover,  these  mechanisms  are  expected  to  be  responsible  for 
transporting  fluid  in  both  directions  across  the  alveolar  epithelium  of  the  normal 
lung  to  maintain  a  normal  amount  of  alveolar  surface  fluid  (5).  These 
mechanisms  must  of  course  work  in  concert  with  those  mechanisms  responsible 
for  maintaining  the  integrity  of  the  surfactant  layer  which  lines  the  alveolar 
surface.  The  interstitial  fluid  is  reasonably  expected  to  be  held  in  balance  by  the 
lymphatics  which  routinely  remove  any  excess  of  fluid  that  passes  across  the 
capillary  endothelium  as  a  result  of  hydrostatic  and  protein  osmotic  pressures. 
The  role  of  the  lymphatics  in  the  lung  fluid  balance  has  already  been  described  in 
considerable  detail  by  Matthay  (12)  and  Staub  (13).  In  the  event  of  a  considerable 
increase  in  the  capillary  blood  pressure,  however,  the  lymphatic  fluid  clearance 
cannot  match  the  flow  into  the  interstitial  space  and  interstitial  edema  quickly 
develops.  Also,  if  the  increase  in  capillary  pressure  is  sufficient,  the  active 
transport  mechanisms  in  the  epithelium  can  most  certainly  become  overpowered 
and  alveolar  edema  results.  Similarly,  in  the  event  of  increased  epithelial 
permeability,  the  mechanisms  of  active  fluid  transport  across  the  epithelium  likely 
become  ineffective  (2)  and  so,  once  again,  the  passive  forces  affecting  fluid 
movement  across  the  epithelium  can  become  dominant  and  result  in  alveolar 
edema.  Once  alveolar  flooding  has  occurred,  if  hydrostatic  pressures  and 
membrane  permeabilities  are  restored  to  normal,  it  appears  that  fluid  is  rapidly 


151 


taken  up  into  perivascular  fluid  cuffs  as  described  by  Gee  et  al.  (6)  and  bi¬ 
directional  fluid  transport  across  the  epithelium  continues  and  once  again 
dominates  the  lung  fluid  balance  mechanisms.  Once  the  excess  fluid  is  moved 
into  the  interstitial  space  it  can  be  cleared  fairly  rapidly  either  via  the  lymphatics 
or  by  crossing  the  endothelium  to  enter  the  blood.  Thus,  the  lung  is  returned  to  its 
normal  state  and  is  able  once  again  to  continue  in  its  primary  function  of  gas 
exchange. 
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5.0  Summary 

Light  microscopy  of  the  lung  tissue  after  either  serum  or  serum  plus  PAF 
instillation  showed  no  morphological  evidence  of  inflammation  and  that  the 
alveolar  epithelium  and  capillary  endothelium  were  apparently  intact  4  hours  after 
fluid  instillation.  Also,  the  125I-albumin  clearance  rate  was  ~3%/hour  regardless 
of  whether  the  instillate  was  serum,  serum  plus  PAF,  SLS,  or  SLS  plus  PAF. 
These  observations  support  that  our  belief  that  NMR  and  gravimetric 
measurements  were  done  in  normal  lungs  of  a  cat,  and  should  have  direct 
physiological  application. 

The  proton  NMR  images  that  we  have  obtained,  have  enabled  us  to 
observe,  non-invasively,  the  time  course  of  alveolar  fluid  clearance  in  living  cats. 
Moreover,  by  employing  both  and  2H  NMR  techniques  we  have  been  able  to 
quantify  the  absorption  and  secretion  components  of  this  net  fluid  clearance.  The 
clearance  of  a  serum  (or  serum-like)  solution  from  the  alveolar  spaces  was 
observed  to  proceed  at  a  rate  of  approximately  17%/4  hours  (ti/2  ~  11  hours),  and 
was  the  net  result  of  fluid  absorption  and  secretion  rates  of  55%/4  hours  and 
38%/4  hours,  respectively.  The  quoted  rates  of  fluid  movement  are  expressed  as 
percentages  of  the  initially  instilled  amount.  These  results  are  consistent  with  a 
model  in  which  the  net  clearance  of  alveolar  fluid  is  primarily  the  result  of  solute- 
coupled  ion  transport,  with  passive  fluid  movements  being  of  secondary 
importance. 

With  the  addition  of  PAF  to  the  serum  instillate,  the  clearance  was  bi- 
phasic  with  an  initial  rapid  clearance  with  a  1 1/2  of  only  ~30  minutes,  followed  by 
a  slower,  more  sustained  clearance  similar  to  that  observed  after  the  instillation  of 
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serum  alone.  Similarly,  in  the  first  hour  after  fluid  instillation,  the  rate  of 
absorption  of  SLS  was  increased  from  28%/hour  to  54%/hour  by  the  addition  of 
PAF,  but  the  fluid  secretion  rate  remained  unchanged  at  approximately  20%/hour. 
However,  the  effect  of  PAF  on  the  clearance  of  serum  was  considerably  more 
pronounced  than  that  observed  with  SLS,  as  demonstrated  by  their  respective  net 
4  hour  clearances  of  65%  and  30%,  respectively.  This  dependence  of  PAF  action 
on  the  instilled  fluid  (serum-based  or  SLS-based)  indicates  that  there  are  one  or 
more  components  of  serum,  not  included  in  the  SLS  preparation,  that  act  to 
enhance  the  effect  of  PAF  on  alveolar  fluid  clearance.  Furthermore,  we  have 
proposed  that  PAF  may  act  by  increasing  cytosolic  Ca2+  in  pulmonary  epithelial 
cells  to  enhance  the  absorptive  function  of  the  epithelium,  or  by  increasing  levels 
of  intracellular  cAMP,  thereby  stimulating  epithelial  ion  transport. 

Transverse  relaxation  curves  obtained  from  proton  spin-echo  images  were 
repeatedly  resolved  into  two  relaxation  components,  indicating  that  we  were 
detecting  NMR  signal  from  at  least  two  lung  water  compartments.  We  have 
demonstrated  the  faster  relaxing  of  these  two  components  (giving  rise  to  60%  to 
75%  of  the  NMR  signal)  to  be  associated  with  tissue-bound  fluid  in  the 
intracellular  and  interstitial  spaces  of  the  lung,  and  the  more  slowly  relaxing 
component  to  be  associated  with  the  instilled  alveolar  fluid.  Whereas  the  tissue- 
bound  fluid  had  a  relatively  constant  T2  value  ot  25  msec  to  32  msec,  the  airspace 
fluid  T2  value  varied  with  the  instilled  alveolar  fluid  (serum-based  or  SLS-based) 
and  with  the  mode  of  ventilation  (spontaneous  breathing  or  mechanical  ventilation 
with  PEEP),  and  ranged  from  150  msec  to  544  msec.  These  observations  are 
consistent  with  the  conclusion  that  lung  water  proton  relaxation  rates  are 
proportional  to  the  tissue  surface  area  in  contact  with  the  airspace  fluid,  because 
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relaxation  within  the  lung  arises  from  a  rapid  exchange  between  bulk  water  and 
water  protons  which  are  tightly  bound  to  biopolymer  segments. 

The  combination  of  and  2H  NMR  techniques  we  have  employed  has 
enabled  us  to  quantify  not  only  the  rate  of  alveolar  fluid  clearance,  but  also  the 
rates  of  alveolar  fluid  absorption  and  secretion  in  the  lungs  of  a  living  animal,  for 
the  first  time.  Also,  these  techniques  are  sufficiently  sensitive  that  we  were  able 
to  assess  the  effects  of  PAF  on  the  mechanisms  which  influence  alveolar  fluid 
clearance,  and  which  presumably  play  a  role  in  the  maintenance  of  the  lung  fluid 
balance.  Thus,  future  research  in  this  area  should  have  two  distinct  goals,  namely, 
the  clinical  application  of  the  NMR  imaging  techniques  that  we  have 
developed,  and  the  further  definition  of  alveolar  fluid  clearance  mechanisms  and 
their  modification  with  PAF,  in  animal  models  of  alveolar  edema.  Further  studies 
should  include  monitoring  the  time  course  of  changes  in  the  extravascular  lung 
water  content,  and  of  proton  transverse  relaxation  times  of  lung  fluid,  in  patients 
with  chronic  lung  diseases  as  well  as  during  the  resolution  phase  of  acute 
pulmonary  edema.  Future  research  with  animal  models  of  alveolar  edema  should 
also  include  evaluations  of  the  effects  of  agents  which  block  the  active  cellular 
transport  and/or  uptake  of  Na+,  CP  and  Ca2+  ions,  on  the  absorptive  and  secretory 
functions  of  the  pulmonary  air-blood  barrier.  In  addition,  monitoring  the 
clearance  of  Na+  from  the  airspaces,  and  the  effect  of  PAF  on  this  clearance, 
should  be  attempted  with  23Na  NMR  techniques. 

This  project  has  demonstrated  the  unique  capabilities  of  NMR 
measurement  techniques  for  monitoring  the  dynamics  of  extravascular  lung  fluid, 
non-invasively ,  in  a  living  animal.  Not  only  have  these  NMR  techniques  enabled 
us  to  observe  separately  the  absorption  and  secretion  of  airspace  fluid,  but  through 
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transverse  relaxation  time  measurements  have  also  enabled  us  to  distinguish 
between  tissue-bound  and  airspace  fluids.  Furthermore,  the  ability  to  monitor  bi¬ 
directional  fluid  movements  across  the  pulmonary  air-blood  barrier  may  prove  to 
be  a  powerful  research  tool.  As  a  result,  it  is  hoped  that  the  successful  completion 
of  this  project  will  provide  an  impetus  for  future  NMR  studies  of  the  lung  fluid 
balance,  in  both  research  and  clinical  settings. 
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Appendix  1.1:  Detailed  Calculations  of  the  Motions  of  Spins  in  the  Presence 

of  a  Rotating  Magnetic  Field 

In  order  to  describe  the  motion  of  the  net  magnetization  of  a  system  of 
non-interacting  spins  under  the  influence  of  a  rotating  magnetic  field,  it  is  useful 
to  work  in  the  rotating  frame  of  reference  and  to  assume  the  simplest  case  to  start. 
The  simplest  case  is  that  in  which  the  rotating  magnetization,  Bi,  and  the  rotating 
frame  of  reference  have  the  same  angular  velocity,  Q,  which  in  combination  with 
a  static  B0  field,  provides  an  effective  field  in  the  rotating  frame  of  reference,  Beff, 
which  is  defined  to  be  along  the  z'p  axis  as  shown: 


^eff  =  Y  Beff 


Figure  A  1.1.1  Precession  of  the  net  magnetization,  M0,  about  an  effective 
magnetic  field  in  the  rotating  frame  of  reference,  Beff,  defined  to  be  along  the  zp' 
axis. 

The  magnetization  components  in  the  Xp'yp'  plane  and  along  the  zp  axis  are 
described  by: 


Mx'  =  Mo  sin0  cos(COit)  ip' 


(Al.1.1) 
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(Al.1.2)  My- =  -  M0  sin9  sin(COit)  jp' 

(A1 .1 .3)  Mz<  =  M0  cos0  kp' 


where  ip',  jp1,  and  kp'  are  unit  vectors 
along  the  xp',  yp'  and  Zp'  axes 
respectively. 

To  expand  this  description  to  the  more  general  case  in  which  the  effective  field  is 
not  aligned  with  the  z  axis  of  the  coordinate  system,  we  can  rotate  the  coordinate 
system  through  an  angle  9.  This  coordinate  rotation  allows  us  to  describe  the  case 

in  which  the  net  magnetization  was  initially  parallel  to  B0  along  the  z  axis,  before 
the  application  of  Beff  at  the  angle  0  from  the  z  axis,  as  shown: 


Figure  Al.1.2  Precession  of  the  net  magnetization,  M0,  about  an  effective 
magnetic  field  in  the  rotating  frame  of  reference,  Beff. 
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Equations  which  describe  a  rotation  of  the  (xp',yp',Zp')  coordinate  system  to  the 
(xp,yp,zp)  system  are  as  follows 


(A1 .1 .4)  Xp  =  Xp'  cos0  -  Zp'  sin0 

(A  1.1.5)  Zp  =  zp' cos0  +  Xp' sinO 

(A  1.1.6)  yp  =  yp' 

Applying  this  coordinate  transformation  to  equations  A  1.1.1  to  A  1.1. 3  yields 


(Al.1.7) 

Mxp  =  M0  (cos(COit)  -  1)  sin0  cos0 

(Al.1.8) 

Myp  =  -M0  sin0  sin(COit) 

(Al.1.9) 

Mzp  =  M0  [cos20  +  sin20  cos(OOit)] 

As  a  result,  the  magnetization  component  in  the  xpyp  plane,  under  the  influence  of 
a  rotating  Bi  field,  can  be  described  by  the  expression 

(A1 .1 .10)  Mxpyp  =  M0  sin0  [  cos20  (cos(COit)  -  l)2  +  sin2(C0it)](1/2) 

where  the  values  of  cos0  and  sin0  are  defined  in  terms  of  the  applied  fields  as 


follows: 
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Figure  Al.1.3  Definition  of  the  effective  magnetic  field  in  the  rotating  frame  of 


reference. 

(Al.1.11) 

Beff=[(B0  +  0/Y)2  +  Bi2](l/2) 

(Al.1.12) 

sin0  =  Bi/Beff 

(Al.1.13) 

cos0  =  (B0  +  Q/yyBeff 
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Appendix  3.1: — Deuteron  NMR  Techniques  Employed  for  Monitoring 

Instilled  Alveolar  SLS.  or  ST  S  plus  PAF 

In  order  to  determine  the  optimal  2H  NMR  technique  for  monitoring  the 
spatial  distribution  of  instilled  alveolar  fluid,  we  assessed  the  suitability  of  gradient- 
recalled  echo  and  free-induction  decay  (FID)  imaging  schemes.  The  8  msec 
minimum  required  by  our  NMR  hardware  to  form  a  gradient-recalled  echo  gave  rise 
to  insufficient  sensitivity  for  gradient-recalled  echo  imaging.  This  time  is 
considerably  longer  than  the  transverse  relaxation  time,  T2*,  with  which  the  echo 
amplitude  decays,  and  which  for  SLS  deuterons  in  the  lung  is  roughly  3  msec. 

The  FID  imaging  scheme,  although  being  more  prone  to  image  artifacts  than 
gradient-recalled  echo  imaging  (4),  provided  a  factor  of  14  greater  sensitivity  than 
the  latter  method,  sufficient  sensitivity  to  obtain  images  of  instilled  SLS  in  the  lungs 
of  a  living  cat.  The  improved  sensitivity  resulted  from  the  ability  to  sample  the  peak 
signal  amplitude  only  30  jlls  after  the  RF  excitation  pulse  was  transmitted.  As  a 
result,  the  2H  NMR  signal  intensity  was  not  degraded  by  the  rapid  deuteron 
transverse  magnetization  decay  in  the  lung.  Images  were  obtained  by  applying  a 
reconstruction  from  projections  technique  to  eight  NMR  signal  projections  from  a 
transverse  imaging  slice.  The  fact  that  only  eight  projections  were  used  limited  the 
spatial  resolution  to  roughly  3  mm.  Because  the  NMR  sensitivity  of  deuterons  is 
only  1%  of  that  of  protons,  it  was  necessary  to  average  512  FIDs  for  each 
projection,  over  a  period  of  51  seconds,  to  achieve  a  projection  signal  amplitude  to 
noise  ratio  of  approximately  10:1  immediately  after  SLS  instillation.  Although  an 
increase  in  the  number  of  projections  and  signal  averages  would  have  improved  the 
image  quality,  the  total  number  of  acquisitions  was  limited  to  4096  (8  projections  x 
512  averages)  by  the  capacity  of  the  computer  memory  allocated  for  storing 
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magnetic  field  gradient  values  for  each  acquisition.  The  image  was  reconstructed  by 
means  of  a  back-projection  algorithm  (2,  3). 

To  achieve  the  maximum  sensitivity  possible  to  2H  NMR,  and  hence  the 
maximum  precision  in  the  measurements  with  the  hardware  employed  for  this 
study,  it  was  necessary  to  forgo  spatial  resolution  altogether  and  measure  the  global 
2H  NMR  signal  intensity  with  a  spectroscopic  technique  similar  to  that  used  by 
Kim  et  al.  (1).  To  maximize  the  weak  deuteron  signal  in  a  given  time,  the  total 
signal  intensity  was  measured  from  1000  spectral  acquisitions  with  a  repetition  time 
of  100  msec.  When  applied  to  detecting  SLS  immediately  after  instillation  into  a 
cat's  lung  lobe  this  spectral  acquisition  technique  was  able  to  provide  a  signal  to 
noise  ratio  as  high  as  50:1  in  100  seconds. 
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Appendix  3.2: — Proton  Transverse  Relaxation  Time  Measurements  of 
Lung  Water,  after  Instillation  of  either  SLS  or  SLS  plus  PAF 

Measurements  of  lung  water  proton  transverse  relaxation  times  were  carried 
out  after  the  instillation  of  SLS-based  instillates  to  supplement  those  obtained  with 
serum-based  instillates  (Chapter  2).  These  measurements  enabled  us  to  monitor 
changes  in  the  microscopic  environment  experienced  by  lung  water  protons  during 
the  clearance  of  instilled  alveolar  fluid. 

Methods 


Transverse  relaxation  time  measurements  of  lung  water  protons  were 
obtained  from  the  lungs  of  anesthetized,  ventilated  cats,  after  the  instillation  of  either 
SLS  or  SLS  plus  PAF,  as  described  in  Chapter  3.  For  this  study,  transverse 
magnetization  decay  curves  were  described  by  integrating  the  pixel  intensities  of  the 
lung  region  of  interest  in  each  of  16  spin-echo  images.  These  spin-echo  images 
were  acquired  using  a  Carr-Purcell-Meiboom-Gill  (CPMG)  sequence  having  an 
inter-echo  interval  of  18.1  msec.  Spatially  averaged  transverse  relaxation  times 
were  then  derived  for  that  lung  region  of  interest  using  a  non-linear  least  squares 
algorithm  (1). 

Results 


The  decay  of  the  transverse  proton  magnetization  from  the  region  of  interest 
in  the  fluid-instilled  lung,  could  be  reproducibly  resolved  into  two  exponential 
components,  with  characteristic  times  T2f  (fast  relaxing)  and  T2s  (slow  relaxing). 
Measurements  of  T2p  and  T2s  (again  time  averaged  over  the  40  minute  intervals 


used  for  the  integrated  signal  intensity  measurements),  are  shown  as  a  function  of 
time  following  instillation  in  Figure  A3 .2.1. 
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Figure  A3. 2.1  A  comparison  of  the  time  course  of  lung  water  proton  transverse 
relaxation  times  of  the  more  slowly  relaxing  component,  T2S>  and  of  the  faster 
relaxing  component,  T2F-  White  circles  are  used  to  indicate  relaxation  times 
measured  after  the  instillation  of  SLS  plus  PAF  whereas  the  black  circles  indicate 
measurements  obtained  after  the  instillation  of  SLS  alone.  The  two  time  courses  of 
T2F  values  illustrated  in  this  figure  overlap  so  that  one  obscures  the  other.  All 
measurements  are  averages  over  40  minute  data  acquisition  periods  and  over  4  cats, 
with  error  bars  indicating  the  standard  error  of  the  mean. 


These  results  indicate  that  the  measured  transverse  relaxation  properties  do  not  vary 
systematically  over  the  4  hours  of  observation.  Assuming  that  fluctuations  are  due 
to  measurement  uncertainty  we  can  derive  an  overall  average  (over  4  hours)  of  the 
relaxation  times.  All  of  the  average  relaxation  times  are  expressed  as  the  mean  of  m 
measurements,  plus  or  minus  the  standard  error  of  the  mean.  After  the  instillation 
of  SLS  the  overall  value  of  T2F  was  32  msec  ±  1  msec  (m  =  18),  and  was  similar 
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to  the  value  of  29  msec  ±  2  msec  (m  =  25)  measured  after  instilling  SLS  plus  PAF. 
The  overall  average  value  of  T2S>  on  the  other  hand,  was  544  msec  ±  55  msec  (m  = 
18)  after  the  instillation  of  SLS  alone,  and  337  msec  ±13  msec  (m  =  25)  for  the 
case  when  PAF  was  added  to  the  SLS  instillate.  Throughout  the  4  hour 
observation  period  the  proportion  of  the  total  signal  from  the  more  slowly  relaxing 
component,  ps,  was  25%  ±  2%  (m  =  43),  irrespective  of  whether  SLS  or  SLS  plus 
PAF  was  instilled  . 

Discussion 


The  ability  to  decompose  the  transverse  proton  magnetization  decay  into  two 
reproducible  components  indicates  that  we  were  detecting  signal  from  two  distinct 
fluid  compartments  in  the  lung.  These  two  compartments  can  be  identified  by  their 
transverse  relaxation  times,  T2F  and  T2S,  respectively.  Previously  (6)  we  have 
associated  the  faster  relaxing  water  protons  with  tissue  bound  fluid  in  the 
intracellular  and  interstitial  spaces  of  the  lung,  whereas  the  more  slowly  relaxing 
water  protons  were  associated  with  the  airspace  fluid. 

Comparison  of  the  present  results  with  those  from  our  previous  study 
demonstrated  that  spatially  and  time  averaged  values  of  T2F  were  similar  for  SLS- 
based  and  serum-based  instillates  (range  of  25  msec  to  32  msec).  Moreover,  the 
small  difference  in  T2F  observed  after  the  instillation  of  either  SLS  (32  msec  ±  1 
msec)  or  SLS  plus  PAF  (29  msec  ±  2  msec)  is  not  regarded  as  significant  and  does 
not  support  a  conclusion  of  a  PAF  dependence.  In  contrast,  the  similarly  averaged 
T2S  value  measured  after  the  instillation  of  serum-based  instillates  (range  of  150 
msec  to  170  msec)  was  considerably  shorter  than  that  with  SLS-based  instillates 
(range  of  337  msec  to  544  msec)  (6).  This  difference  is  possibly  a  consequence  of 
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the  different  modes  of  ventilation  employed  with  the  two  instillate  types  since  it  is 
expected  that  T2S  values  are  elevated  for  fluid  in  larger  airways  (6).  For  example, 
the  cat  was  allowed  to  breathe  spontaneously  after  instilling  a  serum-based  solution 
(6),  whereas  after  the  instillation  of  an  SLS-based  solution  (Chapter  3)  mechanical 
ventilation  was  employed  with  a  positive  end-expiratory  pressure  of  2  cmH20, 
resulting  in  an  increased  average  lung  volume.  Similarly,  because  PAP  has  been 
demonstrated  to  cause  bronchoconstriction  in  a  variety  of  animal  species  (5),  the  use 
of  mechanical  ventilation  is  suspected  of  enabling  us  to  distinguish  between  T2S 
values  measured  with  either  SLS  (544  msec  ±  55)  or  SLS  plus  PAF  (337  msec  ± 
13  msec)  in  the  airspaces. 

The  dependence  of  T2S  values  on  the  airway  diameter  is  consistent  with  the 
conclusion  that  T2S  arises  from  a  fast  exchange  between  the  bulk  airspace  water  and 
a  second  faster  relaxing  water  component,  possibly  that  adsorbed  to  cellular  surface 
proteins  as  suggested  by  Kveder  et  al.  (2).  This  conclusion  is  also  supported  by  the 
fact  that  T2S  does  not  appear  to  have  a  strong  dependence  on  the  T2  of  the  bulk 
airspace  fluid.  Our  results  indicate  that  T2S  is  not  sensitive  to  the  considerable 
increase  in  the  airspace  fluid  protein  concentration  (and  corresponding  decrease  in 
the  bulk  airspace  fluid  T2  (4))  that  occurs  over  4  hours  of  fluid  clearance  (3). 
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Appendix  3.3:  Derivation  of  the  Equations  used  to  Compute  the 

Time  Course  of  Alveolar  Fluid  Absorption  and  of  Alveolar  Fluid 

Secretion 
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The  derivation  of  the  equations  employed  in  Chapter  3  to  compute  the  time 
course  of  alveolar  fluid  absorption,  A(t),  and  of  alveolar  fluid  secretion,  S(t), 
depends  on  the  premise  that  the  observed  clearance  of  2H  from  the  airspaces  is  only 
the  result  of  alveolar  fluid  absorption.  Assuming  this  premise  to  be  true,  the 
quantity  of  2H,  Qabs>  that  is  absorbed  from  the  lung  over  a  period  of  duration,  t,  is 
given  by 

(A3 .3.1)  Qabs(t)  —  Qd(0)  A(t) 

where  Qd(0)  is  the  quantity  of  2H  initially  instilled  into  the  lung.  The  quantity  of 
residual  alveolar  2H,  Qd(0,  is  therefore 

(A3 .3 .2)  Qd(0  =  Qd(0)  -  Qd(0)  A(t) 

Equation  A4.2  can  be  rearranged  to  yield  an  expression  for  A(t) 

(A3. 3 .3)  AU  =  1'W0) 

In  order  to  derive  a  similar  expression  for  S(t)  however,  we  must  first  look 
at  the  net  clearance  time  course,  N(t).  The  net  amount  of  fluid,  QNet>  remaining  in 
the  lung  at  a  time,  t,  after  fluid  instillation,  is  given  by 


(A3. 3 .4) 


N(t)  =  Qo(t)  +  Qh(0 
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where  Qh(0  is  the  quantity  of  residual  alveolar  lH.  Moreover,  the  quantity  of  fluid 
instilled  into  the  lung  initially  is  given  by,  where 

(A3 .3 .5)  Qinst  =  Qd(0)  +  Qh(0) 


N(t)  is  simply  the  amount  of  residual  alveolar  fluid  relative  to  the  instilled  amount, 
and  can  therefore  be  expressed  as 

(A3 .3 .6)  N(t)  =  aDCO.QHCO 


(A3 .3 .7) 


_  Qd(0  Qd(0)  .  Qh(0  Qh(0) 
m)~Qn( 0)  QInst  +Qh(0) 


For  the  study  discussed  in  Chapter  3,  Qh(0)  and  Qd(0)  are  equal,  and  the  ratios 
QH(0)/Qinst  and  QD(0)/QinSt  are  both  equal  to  1/2.  Asa  result 

(A3.3.8)  N(')  =  2[§M  +  §f§] 

Now,  returning  to  the  derivation  of  an  expression  for  S(t),  we  know  from  the 
definitions  of  N(t),  A(t)  and  S(t),  that 

(A3 .3 .9)  N(t)  =  1  -  A(t)  +  S(t) 

which  can  be  rearranged  to  give  the  following  expression  for  S(t). 


(A3 .3. 10) 


S(t)  =  N(t)  + A(t)-1 
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Substituting  for  A(t)  and  N(t)  from  equations  A3 .3 .3  and  A3 .3.8  yields 


(A3 .3. 11) 

sm  _ 1  r  Qd(o  ,  qH(o  1  qd(o 

“  2  LQd(O)  +  Qh(0)J  •  Qd(0) 

(A3 .3. 12) 

_  i  r  Qh(o  qd(o  ~i 
b(t)  -  2  LQh(O)  -  Qd(0)  J 

The  *H  and  2H  NMR  signal  intensities,  In(t)  and  Io(t),  are  assumed  to  be 
proportional  to  the  quantities  of  *H  and  2H  in  the  lung  region  of  interest,  so  that 

(A3 .3. 13) 

Ih(0  Qh(0 

Ih(0)  -  0h(0) 

(A3 .3. 14) 

Id(0  Qd(0 

Id(0)  -  Qd(0) 

By  substituting  the  relative  NMR  signal  intensities  for  the  quantities  of  *H  and  2H 

in  the  lung  at  a  time  t  after  fluid  instillation,  the  expressions  for  N(t),  A(t),  and  S(t), 

become 

(A3 .3. 15) 

N(t)  -  2 [lD(0)  +  IH(0)] 

(A3 .3. 16) 

A«  = 

(A3 .3. 17) 

1  r  !h(0  ID(t)  "I 

S(t)-  2LIh(0)‘ID(0)J 
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Appendix  4.1:  Flow  chart  of  the  Mathematical  Model  of  the  Lung  Fluid  and 
Solute  Exchange 

Definition  of  variables: 

abs+  =  the  additional  epithelial  fluid  absorption  caused  by  PAF 
abs_norm  =  the  normal  epithelial  fluid  absorption  rate 
ap  =  the  pore  area 

asd  =  the  effective  area  for  solute  diffusion  through  a  single  pore 

asf  =  the  area  for  solute  filtration  through  a  single  membrane  pore 

awf  =  the  area  for  water  filtration  through  a  single  membrane  pore 

a_alb  =  the  diameter  of  an  albumin  molecule 

a_water  =  the  diameter  of  a  water  molecule 

C  =  the  mean  concentration  of  solutes  across  the  membrane 

chod_as  =  the  alveolar  HOD  concentration 

chod_is  =  the  interstitial  HOD  concentration 

decay_time  =  the  characteristic  decay  time  of  the  PAF  effect 

At  =  the  element  of  time  over  which  fluid  and  solute  fluxes  are  integrated 

F_endl  =  the  number  of  endothelial  pores  of  size  #1 
F_end2  =  the  number  of  endothelial  pores  of  size  #2 
F_ep  =  the  number  of  epithelial  pores 
Js  =  the  trans-membrane  solute  flux 
Js_as  =  the  solute  flux  across  the  epithelium 

Js  isl  =  the  solute  flux  across  the  endothelium  through  pore  size  #1 
Js_is2  =  the  solute  flux  across  the  endothelium  through  pore  size  #2 
Jv  =  the  trans-membrane  fluid  flux 
Jv_as  =  the  fluid  flux  across  the  epithelium 

Jv  isl  =  the  fluid  flux  across  the  endothelium  through  pore  size  #1 
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Jv_is2  =  the  fluid  flux  across  the  endothelium  through  pore  size  #2 

lp  =  the  membrane  fluid  permeability 

P_1  =  the  hydrostatic  pressure  on  side  1  of  the  membrane 

P_1  =  the  hydrostatic  pressure  on  side  1  of  the  membrane 

7i_l  =  the  protein  osmotic  pressure  on  side  1  of  the  membrane 

P_2  =  the  hydrostatic  pressure  on  side  2  of  the  membrane 

P_2  =  the  hydrostatic  pressure  on  side  2  of  the  membrane 

n_2  =  the  protein  osmotic  pressure  on  side  2  of  the  membrane 

rc_as  =  the  protein  osmotic  pressure  of  the  alveolar  space  fluid 

7t_is  =  the  protein  osmotic  pressure  of  the  interstitial  space  fluid 

• 

Qfas  =  the  rate  of  fluid  flux  across  the  epithelium 
Qfis  =  the  rate  of  fluid  flux  across  the  endothelium 
Qsas  =  the  rate  of  solute  flux  across  the  epithelium 

Qsis  =  the  rate  of  solute  flux  across  the  endothelium 
qhod_as  =  the  quantity  of  HOD  in  the  alveolar  space 
qhod_is  =  the  quantity  of  HOD  in  the  interstitial  space 
Qlymph  =  the  rate  of  lymph  flow  out  of  the  interstitial  space 
qs_as  =  the  quantity  of  solutes  in  the  alveolar  space 
rp  =  the  radius  of  a  membrane  pore 
a  =  the  solute  reflection  coefficient 

sec+  =  the  additional  epithelial  fluid  secretion  caused  by  PAF 
sec_norm  =  the  normal  epithelial  fluid  secretion  rate 
Sin_as  =  the  solute  flux  rate  into  the  alveolar  space 
Sin_is  =  the  solute  flux  rate  into  the  interstitial  space 
Sout_as  =  the  solute  flux  rate  out  of  the  alveolar  space 
Soutjs  =  the  solute  flux  rate  out  of  the  interstitial  space 
t  =  time  elapsed  since  the  instillation  of  fluid  into  the  lung 
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Vin_as  =  the  fluid  flux  rate  into  the  alveolar  space 
Vin_is  =  the  fluid  flux  rate  into  the  interstitial  space 
Vout_as  =  the  fluid  flux  rate  out  of  the  alveolar  space 
Vout_is  =  the  fluid  flux  rate  out  of  the  interstitial  space 
v_as  =  the  alveolar  fluid  volume 

v_instilled  =  the  volume  of  fluid  instilled  into  the  alveolar  space 
v_int  =  the  interstitial  fluid  volume 

v_int_initial  =  the  interstitial  fluid  volume  of  the  normal  lung 
0)  =  the  permeability  of  the  membrane  to  solutes 
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